Mathematical Modeling of Shear Stress and Direct Shear Test for Compressible Soil: Case of Soil Bordering the Wouri River ()
1. Introduction
Coastal regions have been always attractive strategic sites for construction projects. The great property of the soil in these regions is the compressibility of its structure [1] [2]. Moreover, because of their proximity to the sea, border areas are strategic implantation for both persons and companies [3] [4] [5]. Many borders constructions have been made possible after both earth embankment and compaction. In many cases, the water courses affect the surrounding soils and make them unfit to be used as the foundation stone and sometimes require a special treatment before use [6]. Many authors have focused their researches on the theoretical and experimental content analyses of soil properties. To name just a few, Kozicki et al. [7] proposed discrete modeling results of a direct shear test for granular materials versus finite element results, in which the comparison of results obtained for both discrete and finite elements calculations is made, while Xia and coworkers [8] developed a shear behavior of rock joints under constant normal stiffness (CNS) boundary condition. They limit their studies on 48 groups of shear tests on plaster specimens replicated from three rough rock joints, and as results the effects of joint roughness, normal stiffness, and initial normal stress on mechanical characteristics of joints were significantly investigated. Recently Liu and collaborators [9] proposed a new device designed to make direct shear test on sandstone using a constant normal load, while Larsson and Flansbjer [10] developed an approach to compensate the influence of the system normal stiffness in CNS direct shear tests. In their obvious works, they introduced the effective normal stiffness approach applicable to closed-loop control systems, and they nest established the effects of the system normal stiffness on the applied normal load. More recently, Ouagni et al. [11] proposed a direct oedometer test to show that the soil of coast area is highly compressible and out of the phenomenon of shrinkage-swelling that creates differential settlement on construction.
Despite the above mentioned works, there is not enough consideration on direct shear tests for compressible soil issues to the best of our knowledge which still remains challenging and thus deserves more dissemination. It's worth noting that the coastal region in general, and particularly the Douala town bordering the Wouri river of Cameroon, is congested and the population density is too high. This is why inhabitants and companies assault these areas in order to set buildings and houses. Among some of soil tests necessary to be made before any construction project, there is the bearing capacity test through the penetrometer, and in coastal area, it gives very low results (around 0.01 Mpa), leading these areas to be classified as non-constructible, highlighting then the interest of the present work, which aims to give more details about the results obtained for the compressible soils. One may wonder whether some of soils properties can be deduced through mathematical formulations.
In the present study, a laboratory test was conducted to determine the bearing capacity of the soil which is helpful to dimension structural elements such as footings. The results of the shear test at ten points collected at different depths in the area were then analyzed. Accurate determination of the shear strength parameters (angle of internal friction and cohesion) is a major interest in the design of different geotechnical structures. We recall that the shear strength is the main engineering property of soil, helping engineers to ensure the soil mass stability under structural loads [12] [13]. These parameters can be found either in laboratory or on the site. Then understand and master soil behavior is a great challenge in the civil engineering field. In fact, soil can be considered as an assembly of elements [14] [15] [16], with the presence of voids between particles, which contain water and/or air, making the soil a three-phase material [17] [18] [19]. When collecting samples, we observed multiple layers of soil with different colors and the succession of its structure in its profile, as illustrated in Figure 1. Granular classification enables engineers to predict the soil’s behavior, which has important implications for both the mechanical and erosion properties of sediment [20]. The physical properties of the soil of Douala River can be easily found in [21]. According to these obvious reasons, this study aims to determine the average value of the bearing capacity (a function of the depth h) of the soil through the direct shearing test
by applying the Terzaghi formula, which would be a hand use tool for engineers carrying out a project in those boarders and which could help the town planners in their land use plan in the region. This paper is organized as follows: in the next section we present the materials and method to determine the mechanical parameters of the soil; then follows the computational principles of the direct shear test. In Section 3 we present the
Figure 1. Sample of soil log position P1.
results and discussions. In Section 4 we introduce the mathematical modeling and end up in the final section by our conclusions.
2. Materials and Methods
2.1. Materials and Sampling
The data under considerations in the present work was based on ten samples of undisturbed soil collected at different depths of the coast, positioned at locations points
and distanced from each other by ten to fifteen meters, with their geographical location indicated in Figure 2. The inner dimensions of samples are 60 mm × 60 mm in plan which was also the inner dimension of shear box. The thickness of box was about 50 mm while the thickness of sample is 25 mm. It is important to mention that the increasing of the population density in the center region impulses population to build many stores, which impacts on the load acted on the soil. Therefore in addition to the penetrometer test necessary for constructions, civil engineers need to go deeper in soil investigations and carry on the direct shear test. Soil samples are then collected using a variety of approaches and tools subjecting on the depth of the desired sample and the soil type. Near-surface soils are easily sampled using a trowel, and scoop. At greater depths, continuous flight auger (screw) consisted of a trier and a “T” handle is used, After collecting the samples, we brought them to the laboratory to perform the test, via the phases and process which can be found in [22] [23], the US and UK standards which define how the test should be performed being ASTM D3080, AASHTO T236 and BS 1377-7:1990 respectively.
For each sample collected at each point as enumerated above, we took three samples and applied three different loads
, (1)
Which are normal stresses and next we plotted for each of the ten points three curves, corresponding each one to the corresponding normal stress. For each point, we found the breaking stress
and computed the mean value:
(2)
Therefore having three breaking stresses
corresponding to the three normal stresses at each point, we plotted the coulomb line and we extended to the y-axis. Then we obtained the cohesion
of the soil, while the frictional angle
is given by the slope of the coulomb line. This frictional angle as well as the cohesion of the compressible soil is given by their mean values as follows:
(3)
(4)
To model the function which can best represent the graph of the shear stress, the least squared method is used to adjust the results obtained by experiment.
2.2. Fieldwork and Localization of the Site of Study
The main concern for this ancient soil strength test is to keep intact the samples until the laboratory. Next, one uses the most open method to assess the resistance at the interface among two frames. The geographical coordinates of the area where soil samples are collected is shown in Table 1.
2.3. Computation Principles for Direct Shear Test
The principle of test is illustrated in Figure 3, where the soil sample, confined inside the upper and lower rigid boxes (60 mm × 60 mm), is subjected to the normal load
and is sheared by the shear force
, A being the surface where forces are applied. All the above samples are tested after setting the apparatus, under varying normal loads given by Equation (1), to determine
(a)(b)
Figure 2. (a) Location of the study area in the Littoral region of Cameroon. (b) Points of collect samples.
Table 1. Geodetic coordinates of the site location.
Figure 3. Direct shear test principle [24].
the effects upon shear resistance and displacement, and strength properties such as Mohr strength envelope. As results, the maximum shear stresses versus the vertical (normal) confining stresses are plotted for each sample. Next the maximum shear stress
known as critical shear stress or breaking stress is obtained for a specific vertical confining stress. From the plot, a straight-line approximation of the Mohr–Coulomb failure envelope curve is drawn. The cohesion C and the friction angle φ are then computed from the shear stress equation:
(5)
2.4. Experimental Overview
A civil engineer needs to watch out damages on a structure and master correctly all parameters that can influence the structure stability. When a portion of soil is loaded (as illustrated in Figure 4), the ruin can occur by settlement, punching or shearing [25]. Both frictional angle and the cohesion are parameters which are affected.
Knowing these parameters, we assess the bearing capacity of the soil, from where we deduce the ultimate strength called allowable bearing. The cost and safety of any structural building are thoroughly relevant for these values, which has a direct impact of the dimensioning of structural elements and the reinforcement.
As far as cohesion and internal friction are concerned, both simple and elaborated including laboratory methods of determining soil parameters have been developed by geotechnical researchers and engineers (e.g. [26] [27] ). Let us mention that the shear stress of soil can also be determined through triaxial compressing test [28]. The most commonly applied field tests are vane shear test, standard penetration test (SPT), and cone penetration test (CPT) and so on [29] [30]. The mechanical parameters
depends on the size of particles. In this work, the direct shear test was considered due to its reliability and simplicity [31], and the experimental setup is shown in Figure 5.
Figure 4. Failure pattern under a loaded footing [24].
3. Results and Discussions
3.1. Preliminary
This section presents results and discussions about shearing test for coastal areas, which is an experiment to determine the mechanical characteristics of the soil as far as cohesion and frictional angle are concerned, as described in Section 2. For illustrative purpose, Figure 6 depicts the shearing stress versus shearing displacement from where we clearly observe that the evolution is linear and reaches the critical value which is a shear stress of the soil. We then have
representing the critical stress for the above loads given by Equation (1). Following the same process, we have the shear stress for the other samples points (P3…P10) as shown in Figure 7, while the shear stress for the ten points location and their means values are given in Table 2. Therefore, from these values we deduce the values of cohesion and frictional angle
of the soil at these positions as shown in Figure 8 for illustrative purpose for the ten points. The same approaches were applied for the rest of points P3 to P10 and the results are presented in Figure 9.
The means values of mechanicals properties of the soil studied are shown in Table 3 for all the ten points. Statistically, we exploit those values to determine the mean values of cohesion and frictional angles of the particles of soil. Having the standard deviation in good range and the value of the correlation R tends to 1 (R2 in Figure 9) and considering the individual value, we then determine the frictional angle and the cohesion of the compressible soil of the border coast which for our case is Wouri river of Cameroon.
,
(6)
From the curves previously obtained (Figure 8, Figure 10) we can determine the mean values for the breaking stress for each normal stress:
,
,
(7)
By using again the above normal stresses for the entire direct shearing test, we
Figure 6. Shear stress at P1 and P2, tangential stress versus shearing displacement. at point (a) P1:
10−1 Mpa; (b) at point P2:
10−1 Mpa.
Figure 7. Shear stress (peak) at points P3 to P10, tangential stress versus shearing displacement.
Table 2. Shear stress at different points location according to the normal load and their means values.
Figure 8. Morh Coulomb straight line at two locations points, P1: (a): C = 0.60, φ1 = 29.65˚, R2 = 0.9924; P2: (b): C = 0.40, φ2 = 20.820˚, R2 = 0.9893
Table 3. Mean values of frictional angle and cohesion for ten points with relevant density and moisture contents.
have plotted the Mohr Coulomb straight line failure stress vs normal stress, as shown in Figure 10 representing the direct shearing test of the soil with
as cohesion and frictional angle that are the means values. From Figure 8, we deduce average values of the mechanical properties of the soil:
,
. (8)
It is obvious to observe the similarity between these results and those obtained above (5 and 6).
Figure 9. Morh Coulomb straight line at eigth locations points for P3 to P10 in which we can compute
and
.
3.2. Validation of Results
Soils can be visibly known by their sizes as shown in Table 4 [1]. As shown in subsection 3-1, soil bordering the Wouri coast is a sort of clay melt with a low percentage of sand, with clay very soft, since its frictional angle verifies
[24]. The sieve analysis is shown on Figure 11.
Figure 10. Plot of the direct shearing test of the soil for the ten points (P1…P10) considering the mean values.
Table 4. Classification of soil type based on particle sizes and shapes [1].
Figure 11. Particle size analysis of soil borders Wouri River.
Frictional angle and contact force play major role in both soil and foundations classification. Soil of Wouri border of Douala has been recently classified as highly compressible (illites), medium clay according to the range of the compression index Cc [11] and is amongst the most important type of mineral clay. This type of soil can be easily rolled and molded. Shear failure takes place when both cohesive resistance and friction (stress and friction) are not sufficient (see Figure 12). Nevertheless, the normal effective stress rather than total (conventional) normal stress is ideal to evaluate the shear resistance.
As one can see from Figure 8, the relationship between
and
is the affine function, meaning that Equation (5) is satisfied.
As an example, if a structure with a gravity-based foundation is built on this site (see Figure 13), with the applied normal effective stress (
) and shear stress (
) at location P5,
and
, respectively, as obtained experimentally. In order to evaluate whether the soil’s shear stress is sufficient at location P5, From Equation (5), one can have
. It is then obvious that,
and therefore, the shear stress of the soil is unsatisfactory, and the construction will collapse.
The failure criterion is often expressed by “Mohr–Coulomb failure criterion”, also called coulomb line which describes the boundary between soil’s linear-elas- tic and plastic behavior in terms of effective stress parameters, as shown in Figure 10 and Figure 11. It also shows the frontier between elastic and plastic domains (See Figure 8). We note that shear stress of clay is triggered even at zero normal effective stress when the past stress-strain history allows it. From the coulomb-straight-line shown in Figure 8, the extension of the straight-line on
Figure 12. Shear failure in Foundation... [26].
Figure 13. Illustration of shear stress of soils [24].
y-axis doesn’t pass by the origin; proving the lower sand percentage in this soil. Both the mean values
referred to the strength parameters of soils.
· Case of saturated clays,
, leading to the undrained strength
,
Here, based on the method illustrated by Prandtl [32] to study the penetration of hard bodies into softer materials, Terzaghi [33] proposed that the bearing capacity of shallow foundations at a depth h measured from the ground surface can be determined as an ultimate effective pressure
reflecting the general shear failure mode as shown in Figure 4.
, (9)
being constants representing the influence from unit weight, overburden pressure, and cohesion, respectively, and they are functions of friction angle of sands and can be checked from relevant geotechnical handbooks. (The reader can referred to the nomenclature presented in this document for details and can also have an excellent text in Ref. [34] [35] ).
· Illustrative example
For rectangular footings, we denote B and L the width and the length respectively,
,
,
. (10)
Due to the fact that the study was carried out at the boarder of the coast, when the underground water at a depth D is below or above the footing, the effective unit weight of the soil in the third term may be replaced by a weighted average unit weight and (10) becomes:
(11)
Equation (12) enables to calculate the bearing capacity of the soil and hence allowable bearing capacity at any depth h. Let us consider a classical case for the Equation (12) such as
. One obtains:
(12)
,
(13)
Let’s h = 2.00 m. The characteristics of soil obtained previously are
,
. The rest of the parameters values are founded by NAVFAC DM-7.2 [36] ) as follows:
,
,
a) When the soil is saturated, the underground water is below the footing;
,
(14)
b) When the underground water is above the footing,
,
(15)
Having the means values of mechanical properties through the direct shearing test and considering that the bearing strength depends only on the depth h; Equation (12) will be an important tool for engineers, who will do only the dynamic penetrometer test and confront them with these values at any depth in order to have an idea of the range of the allowable bearing strength. This will have a direct impact of the dimensioning the structural elements and then cost and safety. The low value of qallowable shows that shallow foundations and isolated footings are prohibited in the border of coast, since they will be subjected to differential settlement. Meaning that particular precautions are required to avoid these drawbacks.
4. Mathematical Model of Direct Shear Curve Test
4.1. Fitting of the Curves by the Least Square Method
After obtaining the results previously presented by experiments, it is natural to establish an accurate mathematical model that will allow the estimation of the shear stress. When observing different curves obtained through the Casagrande box shown in Figures 14(a)-(c), each of ten points (P1...P10) has been grouped considering three normal loads Equation (1). As one can see, these graphs seem nonlinear with one peak on the x-axis, the maximum representing the breaking stress. The goal in this section is to find the mean curve for the above mentioned normal loads and to determine the mean breaking stress for each one. And next make a comparison with the observational data and generalize the characteristics of the shearing test for the soil in the coast zone. We point out three methods of
Figure 14. Graphs of the shearing test at the ten sampling points in the area (a) 0.1 Mpa, (b) 0.2 Mpa; (c) 0.3 Mpa.
adjustments which highlight more details to our strategy: (a). Fourth order Fourier; (b). Second order exponential; (c): Nonlinear model, among many optimizations methods that are used. The best one being determined by these values: SSE sum of squares due to error, RMSE root mean squared error, which are called goodness of fitting
4.1.1. Fourth-Order Fourier Method (FOFM)
A Fourier series is a periodic function composed of harmonically related sinusoids, combined by a weighted summation. Let’s denote a0 a mean value of f(x). and ai, bi, the fundamental and the upper harmonics respectively, ωx denotes the corresponding phase for the variable x. To increase the accuracy of the approximation, the fourth-order equation has been chosen due to the fact that it is taking in account more information about the measure achieved. This interpolation is shown on Figure 15 corresponding to the three normal stresses with the function defined as follows:
(16)
where
,
and
are constant real parameters to be determined.
Figure 15. Pictures of the shearing test after fitting, (a) 0.1 Mpa, (b) 0.2 Mpa; (c) 0.3 Mpa with FOFM.
4.1.2. Exponential Function Model (EFM)
Exponential function is chosen as a combination of one or two exponential terms linearly independent.
This approximation takes into consideration all variations and curvatures in the parameters space. The graphs obtained are presented in Figure 16 by using the function defined as follows:
(17)
being also the constant real parameters.
4.1.3. Nonlinear Polynomial Function (NPFM)
This method is considered in order to point out some complex aspects input- output which aren’t linked by a constant. A nonlinear function is defined by terms having a variable of degree two or higher. The relevant function in view is presented as follows and the graphs obtained are shown in Figure 17.
(18)
With
the real constant parameters
4.2. Comparative Study
The main results from the methods presented above are recorded in Table 5. The comparative analysis is made in order to choose the fitting method which
Figure 16. Pictures of the shearing test after fitting, (a) 0.1 Mpa, (b) 0.2 Mpa; (c) 0.3 Mpa with EFM.
Figure 17. Pictures of the shearing test after fitting, (a) 0.1 Mpa, (b) 0.2 Mpa; (c) 0.3 MPa with NPFM.
Table 5. Results summary from different fitting methods (FOFM, EFM, NPF)
represents the best approximation of the direct measure of the shear stress presented in the first part of this study. We clearly observe that for the fourth Fourier order method the values of the breaking stress are too closed to the experimental values. For the EFM and NPFM, the results are obtained and the output for each normal stress is derived but with low accuracy compared to the FOFM. From this analysis, it’s obvious to reveal that the FOFM stands as the best fitting method to model the shear stress. Nevertheless, even though, the EFM and NPFM appear unsuitable for the coast, they may be exploited for other type of soils.
5. Concluding Remark
In this work, through a direct shear test, we have computed the characteristics of compressible soil (internal friction angle and the cohesion) and developed a corresponding mathematical model for the coastal area. Firstly, we provided by experiments the mechanical properties of the soils. It’s come out as far as cohesion and frictional angle are concerned that the soil is a mixture of clay and sand with 0 < C < 0.48 and cohesive (
) respectively. The case of Wouri borders in Cameroon is chosen to validate the results and the outcome has been satisfactory. These values will bring many conveniences to engineers of geotechnical field whenever they have identified that the soil is highly compressible. It will help to plan the construction of an appropriate structure with the soil conditions without having to make all tests for determining soil properties by doing only the penetrometer test and compare the values at different depths. Secondly, the comparative study has shown that the best fitting method to model the shear stress is the fourth order Fourier function. Some advantages can be listed through these results such as helping engineers to take decision faster and more accurately, typically offer convenience and cost advantages over other means of obtaining the required information on reality, helping civil engineers and architects to rapidly access the type and capacity of building as far as the elastic and plastic domains are known. It is relevant to point out that during this test the water pressure and drainage conditions are uncontrolled, additionally the non-uni- formity of the failure plan. We then have only an approximation of pure shear specimen, so it will be interesting to carry out the same test with another test as triaxial method and compared the results.
Nomenclatures
C: Cohesion (kPa);
R: Coefficient of correlation;
: Bearing capacity factors;
: Shape factors;
: Depth factors;
: Load inclination factors;
h: Depth of embedment of foundation (m);
Greek Symbols
: Dry unit of the soil (kN/m3);
: Moisture content, percentage (%);
: Normal stress (kN/m2);
: Critical (breaking) shearing stress (kN/m2);
: Angle of internal friction (degree);
: Effective width of a rectangular foundation (m);
: Effective unit weight of the soil (kN/m3).