Abstract
Because indentation tests which are used in several engineering fields allow obtaining material strength parameters in a fast, reliable and reproducible manner, the present study was initiated with the aim of finding out whether such tests could also be used in geotechnical engineering. Experimental results are obtained from undrained indentation tests performed with a rigid flat-tipped, cylindrical indenter on Champlain clay specimens in order to deduce values of elastic moduli and yield stresses. These values are compared to those deduced from unconfined compression tests. Results of an Arbitrary Lagrangian–Eulerian (ALE) based finite element analysis that simulates such test are also presented. It is assumed that the clay behaves as a linear-elastic, perfectly plastic material obeying von-Mises yield criterion. A frictionless contact interface is chosen to realistically model interactions on the indenter-clay, clay-platen and clay-ring interfaces. The paper presents distributions of contact pressures beneath the indenter and along lateral and lower boundaries. Typical contours of von Mises deviator stress and equivalent plastic strain corresponding to different indentation depths in the clay specimen are presented. It is also shown that the presence of hairline cracks on the clay along the perimeter of the indenter and the limited thickness of the sample affect the deduced values of Young’s modulus and yield stress.
Similar content being viewed by others
Introduction
Indentation tests have been used for more than a century to characterize the mechanical properties of materials. Tests involve the application of load and the measurement of the depth of penetration or the size of the residual imprint after unloading. The value of the mean contact pressure \(p_{m}\) at which there is no additional increase with indenter load is related to the hardness number H [18, 29, 48]. The mean contact pressure is equal to the hardness for indentation methods that employ the projected area of contact. Experiments have repeatedly shown that the mean contact pressure \(p_{m}\) or the hardness H is directly proportional to the yield stress \(\sigma_{y}\) of the indented material,
where C is a constraint factor, the value of which depends on the type of specimen and the geometry of the indenter [48]. For instance, \(C \approx 3\) for materials with a large ratio \(E/\sigma_{y}\), and \(C \approx 1\) for low values of \(E/\sigma_{y}\), with E = Young’s modulus of indented specimen [48]. Various models which combine the elastic and plastic responses of materials have been proposed for the determination of the hardness, and for the estimation of their mechanical properties, on the basis of indentation test results obtained with indenters of different geometries (i.e., conical, pyramidal, spherical, or flat- ended cylindrical). These models variously describe the response of indented specimens in terms of elastic deformations, radial compressions, and slip lines.
The indentation problem was also successfully approached using the finite element method of analysis in which indentations are simulated by means of elastic or rigid indenters on elastic, elastic–plastic, and viscoelastic materials, with or without adhesion or friction [5, 6, 10, 21, 22, 27, 40].
Although indentations with spherical and sharp indenters have been widely used for the characterization of material properties [15, 17, 28, 31], application of flat-ended cylindrical indenters can be also found in the literature (See, for example, [13, 40]). While the contact surface area increases with increasing load for spherical and sharp indenters, it remains constant for flat-ended cylindrical indenters; as a result, the mean contact pressure increases linearly with applied load.
A type of indentation test has also been used in Scandinavia since 1915 to estimate the consistency of clays in routine geotechnical studies. Indeed, Olsson [38] described the fall-cone test in which a cone is allowed to fall freely under its own weight from a position at rest. Olsson [38] devised a strength number proportional to the mass of a cone of 60o tip angle which would produce a penetration of 10 mm. Hansbo [23] showed that for a cone of mass m and tip angle α, the undrained shear strength is directly proportional to the square of the penetration. However, as the fall-cone test is a dynamic test in which the cone initially accelerates during its fall and then decelerates to come to rest, the constant of proportionality varies from one clay to another. As a result, the fall-cone test is mainly used as an index strength test, rather than an absolute test, like the field vane and cone penetration tests [54]. Consequently, it is more appropriate to use a static indentation test in which boundary conditions are better controlled.
In this paper a flat-ended cylindrical indentation technique is used for the characterization of the mechanical properties of a natural sensitive clay of eastern Canada. The technique involves indentation of intact specimens of clay by means of a flat-ended steel punch and allows determination of load P versus depth of penetration h curves from which elastic moduli and undrained shear strengths are deduced. These parameters are compared to data obtained from conventional geotechnical laboratory tests, such as unconfined or uniaxial compression tests, miniature vane, and Swedish fall-cone tests.
The paper also presents a brief review of the most common theories used to interpret indentation test results obtained by means of flat-ended cylindrical indenters, and as well as a numerical finite element simulation of indentation tests carried out on the clay.
Indentation theories
Elastic response
The problem of determining the stresses within an elastic half-space indented by a rigid cylindrical punch seems to have been considered first by Boussinesq [9], using the methods of potential theory. Partial numerical results based on Boussinesq’s solution were also derived in the case of the flat-ended punch [36] and a conical punch [35]. However, as Boussinesq’s method of analysis does not lend itself to practical analyses, Harding and Sneddon [24] using integral transforms obtained the general solution for the stresses and displacements generated by an axisymmetric punch on an elastic half-space. The solutions appropriate to the flat-ended cylindrical punch and the conical punch have been discussed in detail by Sneddon [45, 46], respectively. Further, Sneddon [44] gives the complete solution for the stresses and displacements resulting from the indentation of an elastic half-space by a rigid flat-ended cylindrical punch. In addition, Sneddon [43] derived the relation between load and penetration in the axisymmetric Boussinesq problem for a punch of arbitrary shape, using a simpler method than that employed by Harding and Sneddon [24]. The results of Sneddon’s solution for the flat-ended cylindrical punch may be summarized with the help of Fig. 1 (see, also, [5, 19, 40]). The contact area is assumed to remain constant and equal to that of the indenter with radius a. The boundary conditions are:
The first condition indicates that the horizontal boundary outside the loaded area is free of the vertical stress. The second condition specifies that the interface between the indenter and the half-space is smooth. The third condition forces the half-space in contact with the indenter to undergo a uniform downward displacement \(u_{z}\) equal to h. In addition, the vertical stress \(\sigma_{z} \to\)\(\infty\) at \(r = a\) because of the shape of the sharp edge along the perimeter of the indenter. Although this condition causes the development of a localized plastic zone beneath the circular edge of the indenter, it does not invalidate the elastic solution, as discussed by Timoshenko and Goodier [52], Sneddon [43, 44] and Johnson [28].
The mean applied pressure \(p_{m}\), following a depth of penetration h, is given by:
where \(P =\) applied load, and \(\nu =\) Poisson’s ratio. The distribution of the contact stress \(\sigma_{z} \left( {r \le a,\,z = 0} \right)\) is found from:
which shows that there is a square root singularity at \(r = a\). The vertical and horizontal displacements of the horizontal boundary, \(u_{z}\) and \(u_{r}\), are given respectively by:
The remaining elastic relations in terms of cylindrical coordinates \(\left( {r,\theta ,z} \right)\) are the following [17, 19, 44, 46]:
and
with \(\rho = \frac{r}{a}\), \(\zeta = \frac{z}{a}\) and
where \(R = \left[ {\left( {\rho^{2} + \xi^{2} - 1} \right)^{2} + 4\zeta^{2} } \right]^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0pt} 2}}}\), \(r^{ * } = \left( {1 + \zeta^{2} } \right)^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0pt} 2}}}\), \(\xi \tan \theta = 1\), and \(\left( {\rho^{2} + \zeta^{2} - 1} \right)\tan \phi = 2\zeta\).
Elastic–plastic response
An indentation model which also takes into account plastic response is extremely complex to develop due to the nonlinearity of the constitutive relations (see, for instance, [55]). As a result, only approximate models have been proposed for the description of the elastic–plastic response of materials subjected to indentations by punches of different profiles, namely, the spherical cavity expansion model and the slip line theory.
For blunt indenters, including spherical and flat-ended cylindrical punches, the specimen responds in an elastic–plastic manner, and plastic flow is usually described in terms of the elastic constraint provided by the surrounding elastic material. Marsh [37], inspired by the work of Samuels and Mulhearn [41], compared the plastic deformation in the neighborhood of a blunt indenter to that which occurs during the radial expansion of a spherical cavity under internal pressure, an analysis of which was first given by Bishop et al. [7] and Hill [25]. Johnson et al. [30] replaced the expansion of the cavity with that of an incompressible hemispherical core of material subjected to internal pressure, as indicated in Fig. 2. In this model, which has been since referred to as the “expanding cavity model” or ECM, the core pressure is directly related to the mean contact pressure \(p_{m}\), that is,
where \(\alpha\) is the inclination of the indenter with the specimen surface and thus leads to a value of
for the constraint factor. Although the ECM cannot be directly applied to a flat-ended cylindrical punch indentation because the shape of the plastic zone evolves with the depth of penetration and tends to be similar to that around a spherical cavity only at great depth [6, 32, 55], it nonetheless represents a class of popular models that provide the capability to consider the effects of both elastic and plastic response on the evolution of the deformation fields during the indentation process [53]. For a flat-ended cylindrical punch, the relationship in Eq. (10) is usually replaced by the expression for the expansion of a spherical cavity, that is,
leading to
For deeper penetrations, the plastic zone is no longer constrained by the surrounding elastic material, and the volume of the material displaced by the indenter is accommodated by the upward flow around the indenter. The same phenomenon takes place with sharp conical indenters when they first touch the specimen [48]. A cutting mechanism is involved, and new surfaces are formed beneath the indenter as the volume displaced is accommodated by the upward and lateral flow of the plastically deformed material [17]. The constraint factor in this case arises due to flow and velocity considerations [25]. In such situation, slip line theory can be applied to the modelling of the indentation because the response is similar to that which takes place with rigid plastic materials which flow at constant shear stress. Shield [42] showed that the axisymmetric flow of a rigid plastic material could be described by a slip-line field if the material obeyed the Tresca criterion, and assuming the validity of the Haar-von Karman hypothesis which states that the circumferential stress \(\sigma_{\theta }\) is equal to one of the principal stresses in the meridional plane. Shield [42] obtained the solution in the case of the rigid plastic half-space indented by a smooth flat-ended cylindrical punch. The average contact pressure \(p_{m}\) was found to be equal to \(2.845\,\sigma_{y}\) or \(5.69\,k\), where \(k\) is the maximum shear stress (or \(S_{u}\)) in the terminology used in geotechnical engineering). The vertical pressure is equal to \(5.14\,S_{u}\) beneath the perimeter of the punch and the maximum value which equals \(7.2\,S_{u}\) occurs below the axis of symmetry. It should be recalled that Lockett [34], Houlsby and Wroth [26], and Chitkara and Butt [11] also obtained \(p_{m} = 5.69\,S_{u}\) in the case of the smooth flat-ended cylindrical punch. Eason and Shield [14] considered the case of a perfectly rough punch on a rigid plastic half-space and found \(p_{m} = 6.05\,S_{u}\), indicating that the material becomes fully plastic when \(p_{m}\) reaches about \(3\,\sigma_{y}\).
Finite element analysis
As mentioned above, the indentation of an elastic-perfectly plastic material constitutes a very difficult problem to be solved analytically. However, difficulties which arise due to the unknown shape and extension of the plastic field can be minimized by means of a finite element analysis (FEA), provided: (a) an adaptive refined mesh is used to model the indented material, (b) a realistic material response is considered, and (c) the geometry and boundary conditions are properly modelled. In the present paper, an Arbitrary Lagrangian–Eulerian (ALE) based finite element analysis using Abaqus/Explicit is conducted to simulate the experimental investigations performed on natural clay specimens.
All calculations were carried out by means of an axisymmetric model; it consists of a rigid flat-ended cylindrical punch (radius of 9.975 mm and height of 9.95 mm) that penetrates a cylindrical specimen of clay (radius of 31.5 mm and height of 19.0 mm), which is encased in a steel ring (thickness of 2 mm and height of 19.0 mm). The soil specimen and the steel ring rest on a rigid horizontal platen. The analysis was carried out by means of the general-purpose code ABAQUS [1]. Because of the severe distortion that occurs mostly beneath the perimeter of the punch in contact with the clay, modelling was performed using the Arbitrary Lagrangial Eulerian (ALE) auto-adaptive remeshing technique integrated in ABAQUS/Explicit. This technique maintained a high-quality mesh throughout the indentation process.
Instead of using a finite element approach, it is also possible to perform a finite difference analysis for the study of penetration interactions. Indeed, Ghazavi and Tovaioli [20] adopted a finite difference scheme to simulate a three-dimensional model for pile driving analysis using FLAC3D. It is an explicit finite difference program for simulating the response of three-dimensional soil structures built on rock or other soil materials undergoing plastic flow. Further applications of the use of FLAC3D to other pile driving analyses are reported by Tovaioli and Ghazavi [49,50,51].
Geometry of the assembled parts
As just mentioned, the flat-ended cylindrical punch (OFGH in Fig. 3) was modelled as a rigid body. A 0.1-mm radius circular boundary was created along the perimeter of the indenter in contact with the soil, to avoid numerical singularities. The cylindrical specimen (OBCD in Fig. 3) is bounded by a rigid circular ring CM and a rigid bottom platen BC.
Mechanical properties of materials
It is assumed that the clay behaves as a rate-independent isotropic linearly elastic perfectly plastic material obeying the von Mises flow rule. Values of Young’s modulus and yield stress were obtained from a series of unconfined compression tests performed on intact soil specimens. The unconfined compression tests were carried out according to ASTM [2] —Test Method for Unconfined Compressive Strength of Cohesive Soil which insures constant volume conditions, from which a value of 0.5 is deduced for Poisson’s ratio. The average value of Young’s modulus E was found equal to 15 MPa, and Poisson’s ratio \(\nu\) was assumed to be equal to 0.49 in the Finite Element Analyses to simulate constant volume conditions. The yield stress \(\sigma_{y}\) is equal to 113.4 kPa, as also found from the unconfined compression tests. Finally, the densities of the clay and of the steel indenter are considered to be: \(\rho_{clay} = 1.6 \times 10^{ - 9} \,{\text{t/mm}}^{ 3}\) and \(7.83 \times 10^{ - 9} \,{\text{t/mm}}^{ 3}\), respectively.
Interactions between assembled parts
Three surface-to-surface interactions are employed in the simulation. This type of interaction is integrated in ABAQUS/Explicit; it uses a kinematic contact method as a mechanical constraint formulation and a finite sliding formulation in which only the master surface is allowed to penetrate the slave surface. The first interaction occurs between the lateral surface (OFG in Fig. 3) of the rigid indenter which acts a master surface and the free upper horizontal surface of the specimen (OE in Fig. 3) which acts as a slave surface. The second interaction occurs on the interface represented by line BC in Fig. 3, i.e., between the horizontal bottom edge of the specimen which acts as a slave surface and the inner surface of the rigid low platen which acts as a master surface. Finally, a third interaction occurs on the interface represented by vertex CDM in Fig. 3, i.e., between the vertical cylindrical surface of the specimen (CD in Fig. 3) which acts as a slave surface and the inner vertical cylindrical surface of the container (CM in Fig. 3) which acts as a master surface. These aforementioned three interactions use the same frictionless formulation.
Experimental and numerical results
Clay
The clay tested in this study is a sensitive clay from a soil deposit in the town of Beloeil (Quebec) in eastern Canada. The site is located at a short distance from the shore of the Richelieu River, about 40 km to the south-east of Montreal. Undisturbed blocks of clay, measuring 0.3 m in width, were recovered at a depth of 3.5 m in a test excavation.
The most relevant properties of the clay are reported in Table 1, as given by Ewane et al. [16].
The soil is plastic, sensitive, and lightly overconsolidated. The undrained shear strength \(S_{u}\) of the clay was determined by means of unconfined uniaxial compression (UU) tests [2], miniature vane tests (VST) [3]—Test Methods for Laboratory Miniature Vane Shear Test for Saturated Fine-Grained Clayey Soil), and Swedish fall-cone tests [47]—Determination of Undrained Shear Strength and Sensitivity of Cohesive Soils using a Fall Cone Penetrometer).
The average value of \(S_{u}\) (Tresca) equals 56.7 kPa for the UU tests, 33.5 kPa for the VSTs, and 73.4 kPa for the Swedish fall-cone tests. The lowest value of 33.5 kPa is linked to the development of cracks around the blades of the miniature vanes upon insertion of the instrument. Because the average preconsolidation pressure, \(\sigma^{\prime}_{p}\), equals 143.2 kPa, the strength ratio \(S_{u} /\sigma^{\prime}_{p}\) equals 0.396 for the UU tests, 0.233 for the miniature vane tests and 0.512 for the Swedish fall-cone tests. While the value of 0.233 is slightly low for a typical sensitive clay of eastern Canada, the remaining ratios of 0.396 and 0.512 are rather high for an average overconsolidation ratio (OCR) of 6.42. Indeed, a more realistic value would range from 0.25 to 0.30 for such case [33].
A typical undrained uniaxial (UU) compression test result is reported in Fig. 4. Values of Young’s modulus E range between 15 and 20 MPa, and the yield stress \(\sigma_{y}\) (\(\sigma_{y} = 2S_{u}\)) for Tresca criterion varies between 104 and 108 kPa.
Indentation tests
Indentations were carried out using a computer-controlled universal testing machine. The 19.95 mm diameter punch was fixed to the top platen of the loading frame and was coated with graphite powder to minimize the adhesion with the clay specimen. Each clay specimen which measured 63 mm in diameter and 19 mm in height was encased in a polished Teflon coated oedometer steel ring of 2 mm in thickness and 19 mm in height to provide lateral confinement during indentation. The oedometer ring containing the clay specimen was placed on the lower platen of the loading frame. Indentation tests were performed by raising the lower platen at a constant rate of penetration of 0.5 mm/s. The preparation of the soil specimens and their installation in the oedometer rings were made following [4]—Test Method for One-Dimensional Consolidation Properties of Soils using Incremental Loading. The soil surface in contact with the cylindrical punch was trimmed flush with the plane ends of the rings. A wire saw was used for trimming the top and bottom of the specimens to minimize smearing. The load on the indenter and the depth of penetration were continuously recorded during each test. A LVDT with a precision of 0.0075 mm and a total displacement of 5 mm was used to monitor the depth of penetration of the punch. Typical results are presented in Fig. 5.
Examination of the experimental relationships indicates that the curves become practically linear following an initial nonlinear phase, probably caused by surface defects which induce an ill-defined initial contact with the clay surface, as also found in the indentation testing of metals [8, 56].
The slopes \(\Delta P/\Delta h\) of the linear segments of the curves shown in Fig. 5 range between 395 and 485 N/mm, with an average value of 455 N/mm. By assuming, for the moment, that the clay specimen can be considered to represent an elastic half-space, application of Eq. (3) yields \(E = 17.1\,{\text{MPa}}\) for \(P/h = \Delta P/\Delta h = 455\,{\text{N/mm}}\), \(\nu = 0.5\), and \(a = 9.975\,{\text{mm}}\). Such close agreement with the experimentally deduced values of E of 15–20 MPa as found from the uniaxial compression tests (see Fig. 4 and Table 1) is, however, purely coincidental. Indeed, results of FEA computations by Christian and Carrier [12] and Poulos [39] showed that Eq. (3) had to be modified when the thickness of the elastic medium is of limited extent. In such case, Eq. (3) should read
where \(I_{\rho }\) is a geometrical factor that is function of the ratio of the thickness of the medium to the diameter of the punch. For \(H/2a = 19/19.95 \approx 0.95\), \(I_{\rho }\) equals 0.4 for Christian and Carrier [12], and 0.6 for Poulos [39]. Therefore, if \(I_{\rho }\) is taken equal to an average value of 0.5, the value of Young’s modulus E deduced from Eq. (14) reduces to 8.55 MPa, which is about half of the experimental results obtained from the uniaxial compression tests. Such reduced value of E is undoubtedly related to partial damage induced by penetration of the indenter, as evidenced by the progressive development of hairline radial and concentric cracks on the clay surface, around the perimeter of the indenter. Concerning the determination of the hardness H of the clay, for which there is no additional increase in the applied load during penetration, unfortunately the indentation tests were not pursued to sufficiently large depths. This notwithstanding, if the mean contact pressure of 266 kPa which corresponds to the highest applied load of \(83\,{\text{N}}\) (see Fig. 5), is considered to be equal to the hardness H, then application of Shield’s [42] solution, \(H = 2.845\,\sigma_{y}\), allows to determine the yield stress \(\sigma_{y}\) in the case of Tresca criterion. Thus, \(\sigma_{y} = 93.3\,{\text{kPa}}\) for \(H = p_{m} = 266\,{\text{kPa}}\). Such value of \(\sigma_{y}\) yields an undrained shear strength \(S_{u} = \sigma_{y} /2\) or 46.7 kPa for Tresca, and \(S_{u} = \sigma_{y} /\sqrt 3 = 53.9\,{\text{kPa}}\) for von Mises. However, according to Table 1, a value of \(S_{u} = 46.7\,{\text{kPa}}\) is about 82% of that obtained from the uniaxial compression tests.
Results of finite element analyses
FEA results were obtained for particular values (i.e. h = 0.124, 0.243, 0.349, and 0.57 mm) displacements based upon experimental test results. Figures 6 and 7 present contours of the equivalent von Mises stress \(\sigma\) which is defined as
and in terms of principal stresses \(\sigma_{1}\), \(\sigma_{2}\), \(\sigma_{3}\), this definition can be written in a simpler form
While Fig. 8 presents contours of the total amount of equivalent plastic strain \(\bar{\varepsilon }^{pl}\) expressed by:
where \(\left. {\bar{\varepsilon }^{pl} } \right|_{0}\) is the initial equivalent plastic strain which is zero in this present analysis; and \(\dot{\bar{\varepsilon }}^{pl}\) is the equivalent plastic strain rate.
Although the results reported in Figs. 6, 7 and 8 indicate that most of the clay did not reach a fully plastic state for depths of indentation less than or equal to 0.57 mm, the yield stress and plastic strain contours of Figs. 7 and 8 respectively nevertheless point out that severe distortions occur immediately below the rounded perimeter of the indenter. The latter figure also shows the progression of the plastic zone toward the axis of symmetry of the clay with increasing depth of penetration.
The relationship between the computed load P versus the depth of indentation h up to a depth of penetration of 0.57 mm is presented in Fig. 9. It shows that the slope \(\Delta P/\Delta h\) is roughly constant, equal to 709.4 N/mm. Then, the adoption of Eq. (14) for the representation of an elastic response results in \(E/I_{\rho } = 26.7\,{\text{MPa}}\) for \(P/h = \Delta P/\Delta h = 709.4\,{\text{N/mm}}\), \(\nu = 0.49\), and \(a = 9.975\,{\text{mm}}\). Because the value of E used in the FEA equals 15 MPa, it follows that \(I_{\rho } = 0.56\), compared to 0.4 for Christian and Carrier [12] and 0.6 for Poulos [39], as mentioned above. In addition, the value of the constraint factor C equals 3.65 on the basis of the spherical cavity approach of Eq. (13), for \(E = 15\,{\text{MPa}}\), \(\nu = 0.49\) and \(\sigma_{y} = 113.4\,{\text{kPa}}\).
Furthermore, because the P versus h relationship shown in Fig. 9 shows that the initial “elastic” or linear phase terminates at h = 0.05 mm for P = 35.47 N, the corresponding mean contact pressure \(p_{m}\) is 113.5 kPa (Fig. 10). Thus, the value of the constraint factor \(C = p_{m} /\sigma_{y}\)\(= 1.0008\)\(= 113.5/113.4\) at the start of the non-linear response in the P versus h relationship. This observation indicates that the plastic response begins to play a role when the mean contact pressure \(p_{m} \approx \sigma_{y}\), as also found by Riccardi and Montanari [40], and as pointed out several years ago by Tabor [48] and Johnson [28] on the basis of sharp and spherical indenter test results. It should be noted that Fig. 10 also presents the central and the mean contact pressures.
Contact pressure distribution for the characteristic depths of penetration of 0.124, 0.243, 0.349, and 0.57 mm are presented in Fig. 11. Examination of the curves reported in this figure indicates that the contact pressures which generally increase from the axis of symmetry toward the perimeter of the indenter, reach a maximum value for \(r = 8.2\) to 9.7 mm, and thereafter decrease sharply as the rounded perimeter of the indenter is approached.
Contact pressures computed along the bottom platen are presented in Fig. 12. All the relationships reported in this figure are characterized by a classical bell-shaped curve.
Finally, Fig. 13 presents the horizontal contact pressure computed along the external vertical boundary at different vertical nodal elevations. The origin of these elevations is point C in Fig. 3. Examination of all the relationships shows that the contact pressure decreases from the top to the bottom of the medium. In addition, the nodal elevations of the clay specimen along the external vertical boundary remained practically unchanged.
Discussion
Indentation techniques and associated mathematical or numerical models used for the physical characterization of various materials require experimental validation in order to obtain accurate and precise measurements. In this study two simple mathematical models (i.e., the elastic and the elastic–plastic models) were tested in terms of their accuracy in producing values for strength parameters that could match those obtained from conventional geotechnical tests. To simplify analysis, investigators commonly assume large sample thickness and limit indentation depth to less than 10% of the sample thickness. In this study, the sample thickness being of the same order of magnitude of the diameter of the indenter, it was found necessary to take into account the limited thickness of the clay. Although such approach resulted in a better representation of the deformation mechanism involved, it nevertheless showed that the deduced values of Young’s modulus were much lower than those obtained from unconfined compression tests. The probable cause of the reduction was related to the development of concentric and radial cracks surrounding the indenter during the indentation process.
The numerical model used with the auto-adaptive remeshing technique allowed to simulate a deep penetration process. In the initial elastic phase of the simulation, the constant slope of the P versus h curve allowed to determine the value of geometrical factor \(I_{\rho }\). It was found that \(I_{\rho } = 0.56\) compared to 0.4 for Christian and Carrier [12], and 0.6 for Poulos [39]. It was also observed that while the plastic response started to become significant at a pressure \({{p_{m} } \mathord{\left/ {\vphantom {{p_{m} } {\sigma_{y} }}} \right. \kern-0pt} {\sigma_{y} }} \approx 1\), the maximum load P corresponded to \({{p_{\text{max} } } \mathord{\left/ {\vphantom {{p_{\text{max} } } {\sigma_{y} }}} \right. \kern-0pt} {\sigma_{y} }} = 3.3\). Such values which compare well with those proposed by Johnson [28] are also quite similar to the results reported by Riccardi and Montanari [40].
Conclusion
Considering the tests that employ sharp or spherical punches, the flat-ended cylindrical indentation is characterized by a constant indenter-specimen contact area; thus, elastic moduli and yield stress can be directly deduced from load or pressure-penetration curves.
Concerning the indentation tests, these were not taken to sufficiently high loads for the material to reach a full plasticity condition. This notwithstanding, the tests allowed to determine: (a) a minimum value of 266 kPa for the yield stress which permitted to obtain a Tresca shear strength of 46.7 kPa compared to the average value of 56.7 kPa obtained from the unconfined compression test; and (b) a value of 8.55 MPa for Young’s modulus, compared to an average value of 17.5 MPa, also obtained from the unconfined compression tests. The cause for the reduced value of Young’s modulus is attributed to the development of hairline concentric and radial cracks during the indentation process.
The yield stress \(\sigma_{y}\) was evidenced in the P versus h curve obtained from the finite element analysis. A value of 2.86 was derived for the constraint factor C at maximum load, based upon von Mises criterion, compared to 2.845 for Shield’s Tresca slip-line analysis.
References
ABAQUS (2018) Johnston. Dassault Système Simulia Corp, RI
American Society for Testing and Materials, ASTM D2166/D2166M—16, Standard Test Method for Unconfined Compressive Strength of Cohesive Soil. 2016 Annual Book of ASTM Standards, Vol.04.08. ASTM International 100 Barr Drive, West Conshohocken, PA, USA
American Society for Testing and Materials, ASTM D4648/D4648M—16, Standard Test Methods for Laboratory Miniature Vane Shear Test for Saturated Fine-Grained Clayey Soil. 2016 Annual Book of ASTM Standards, Vol.04.08. ASTM International 100 Barr Drive, West Conshohocken, PA, USA
American Society for Testing and Materials, ASTM D2435/D2435M—11, Standard Test Methods for One-Dimensional Consolidation Properties of Soils Using Incremental Loading. 2016 Annual Book of ASTM Standards, Vol.04.08. ASTM International 100 Barr Drive, West Conshohocken, PA, USA
Barquins M, Maugis D (1982) Adhesive contact of axisymmetric punches on an elastic half space: the modified Hertz-Huber’s stress tensor for contacting spheres. Journal de Mécanique Théorique et Appliquée 1(2):331–357
Battacharya AK, Nix WD (1988) Finite element simulation of indentation experiments. Int J Solids Struct 24(12):141–147
Bishop RF, Hill R, Moh NF (1945) The theory of indentation and hardness tests. Proc Phys Soc 57(3):147–159
Briscoe BJ, Fiori L, Pelillo E (1998) Nano-indentation of polymeric surfaces. J Phys D Appl Phys 31(19):2395–2405. https://doi.org/10.1088/0022-3727/31/19/006
Boussinesq J (1985) Application des potentiels à l’étude de l’équilibre et du mouvement des solides élastiques: des notes étendues sur divers points de physique mathématique et d’analyse. Gauthier, Villars, p 736
Cheng Y-T, Li Z (2000) Hardness obtained from conical indentations with various cone angles. J Mater Res 15(12):2830–2855
Chitkara NR, Butt MA (1992) Numerical construction of axisymmetric slip-line fields for indentation of thick blocks by rigid conical indenters and friction at the tool-metal interface. Int J Mech Sci 34(11):849–862
Christian JT, Carrier WD III (1978) Janbu, Bjerrum and Kjaemsli’s chart reinterpreted. Can Geotech J 15:123–128
Delaine-Smith RM, Burney S, Balkwill FR, Knight MM (2016) Experimental validation of flat punch indentation methodology calibrated against unconfined compression tests for determination of soft tissue biomechanics. J Mech Behav Biomed Mater 60:401–415
Eason G, Shield RT (1960) The plastic indentation of a semi-infinite solid circular punch. Zeitschrift für angewandte Mathematik und Physik ZAMP 11(1):33–43
Eswar Prasad K, Chollacoop N, Ramamurty V (2011) Role of indenter angle on the plastic deformation underneath a sharp indenter on representative strains: an experimental and numerical study. Acta Materalia 59(11):4343–4355
Ewane MS, Silvestri V, James M (2018) Indentation of a sensitive clay by a flat-ended circular punch. J Geotech Geoenviron Eng. https://doi.org/10.1007/s10706-018-0561-4
Fischer-Cripps AC (2007) Introduction to contact mechanics, 2d edn. Springer-Verlag, New York
Fischer-Cripps AC (2002) Nanoindentation. Springer-Verlag, New York
Gazonas GA, Wildman RA. The mechanics of axisymmetric indentation revisited. Army Research Laboratory, Aberdeen Proving Ground, MD 21005-5066, ARL-TR-6528, 2013.
Ghazavi M, Tavasoli O (2012) Characteristics of non-uniform cross-section piles in drivability. Soil Dyn Earthquake Eng 43:287–299
Giannakopoulos AE, Larsson PL, Vestergaard R (1994) Analysis of Vickers indentation. Int J Solids Struct 31(19):2679–2708
Guha S, Sangal S, Basu S (2014) Numerical investigations of flat punch molding using a higher order strain gradient plasticity theory. Int J Mater Form 7(4):459–467
Hansbo S. A new approach to the determination of the shear strength of clay by the fall-cone test. In: Royal Swedish Geotechnical Institute, Stockholm, Proceedings No. 14, 1957, p. 47.
Harding JW, Sneddon IN (1945) The elastic stresses produced by the indentation of the plane surface of a semi-infinite elastic solid by a rigid punch. Mathematical Proc. Cambridge Phil. Society, Cambridge University Press, Cambridge, pp 16–26
Hill R (1950) The mathematical theory of plasticity. Oxford University Press, London
Houlsby GT, Wroth CP.Direct solution of plasticity problems in soils by the method of characteristics. In: Proc. 4th Int. Conf. on Numerical Meth. in Geomech., Edmonton, 1982, vol. 3, pp. 1059–71
Hu Z, Lynne K, Delfanian F (2015) Characterization of materials’ elasticity and yield strength through micro-/nano-indentation testing with a cylindrical flat-tip indenter. J Mater Res 30(04):578–591
Johnson KL (1985) Contact mechanics. Cambridge University Press, Cambridge
Johnson KL (1970) The correlation of indentation experiments. J Mech Phys Solids 18(2):115–126
Johnson W, Mahtab F, Williams A (1965) Experiments concerning geometric similarity in indentation. Int J Mech Sci 7(6):389–398
Karimzaleh A, Ayatollahi MR, Alizadeh M (2014) Finite element simulation of nano-indentation experiment on aluminum 1100. Comput Mater Sci 81:595–600
Lange K (1980) Handbook of metal forming. McGraw-Hill Book Company, New York
Leroueil S, Tavenas F, Le Bihan J-PL (1983) Propriétés caractéristiques des argiles de l’est du Canada. Can Geotech J 20(4):681–705
Lockett FJ (1963) Indentation of a rigid/plastic material by a conical indenter. J Mech Phys Solids 11(5):345–355
Love AEH (1939) Boussinesq’s problem for a rigid cone. Quart J Math 10(1):161–175
Love AEH (1929) The stress produced in a semi-infinite solid by pressure applied on part of the boundary. Phil Trans Royal Soc London A 228:377–420
Marsh D (1964) Plastic flow in glass. Proc Royal Society of London A 279(1378):420–435
Olsson J (1921) Method for undersokning av lerors hallfasthetsegnskaper, tillampad vid de geotekniska undersokningarma vid Statens Jarnvagar. Geologiska Forening Stockholm, Forhandlingar 43(5):502–507
Poulos H (1988) Cyclic stability diagram for axially loaded piles. J Geotech Eng. https://doi.org/10.1061/(ASCE)0733-9410(1988)114:8(877)
Riccardi B, Montanari R (2004) Indentation of metals by a flat-ended cylindrical punch. Mater Sci Eng A 381:281–291
Samuels L, Mulhearn T (1957) An experimental investigation of the deformed zone associated with indentation hardness impressions. J Mech Phys Solids 5(2):125–134
Shield RT (1955) On the plastic flow of metals under conditions of axial symmetry. Proc Royal Society London A 233(1193):267–287
Sneddon IN (1965) The relation between load and penetration in the axisymmetric Boussinesq problem for a punch of arbitrary profile. Int J Eng Sci 3(1):47–57
Sneddon IN (1951) Fourier transforms. McGraw-Hill Book Company, New York
Sneddon IN (1948) Boussinesq problem for rigid cone. Math Proc Cambridge Phil Soc 44(4):492–507
Sneddon IN (1946) Boussinesq’s problem for a flat-ended cylinder. Math Proc Cambridge Phil Soc 42(1):29–39
Standard CAN/BNQ 2501-110. Soils-Determination of Undrained Shear Strength and Sensitivity of Cohesive Soils Using a Fall Cone Penetrometer. Bureau de Normalisation et Certification du Québec. 333 Rue Franquet, Quebec, QC, G1P 4C7, Canada, 2014.
Tabor D (1951) The hardness of metals. Oxford University Press, London
Tavasoli O, Ghazavi M (2019) Driving behavior of stepped and tapered offshore piles due to hammer blows. Mar Georesour Geotechnol. https://doi.org/10.1080/1064119X.2019.1609631
Tavasoli O, Ghazavi M (2019) Performance evaluation of the lightweight concrete tapered piles under hammer impacts. Geomech Eng 18(6):615–626
Tavasoli O, Ghazavi M (2018) Wave propagation and ground vibrations due to non-uniform cross-sections piles driving. Comput Geotech 104:13–21
Timoshenko S, Goodier J (1951) Theory of elasticity. McGraw-Hill Book Company, New York
Wang Z, Basu S, Murthy TG, Saldana C (2016) Modified expansion theory formulation for circular indentation and experimental validation. Int J Solids Struct 97–98:129–136
Wood DM (2007) Soil behaviour and critical state soil mechanics. Cambridge University Press, Cambridge
Yu W, Blanchard JP (1996) An elastic-plastic indentation model and its solutions. J Mater Res 11(9):2358–2367
Zeng K, Chiu CH (2001) An analysis of load-penetration curves from instrumented indentation. Acta Mater 49(17):3539–3551
Author information
Authors and Affiliations
Contributions
The corresponding author GAS performed all numerical simulations, drew all the curves of the graphs and drafted the manuscript. MSE provided numerical values of the parameters of the clay’s mechanical model, which have been deduced from experimental tests. VS participated in the theoretical analyses of the spherical cavity expansion model and the expanding cavity model. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Abou-Samra, G., Ewane, M.S. & Silvestri, V. Physical modeling and arbitrary Lagrangian–Eulerian finite element analysis of indentation of a sensitive clay by a flat-ended axisymmetrical indenter. Geo-Engineering 11, 6 (2020). https://doi.org/10.1186/s40703-020-00113-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40703-020-00113-4