# Effects of Flotage on Immersion Indentation Results of Bone Tissue: An Investigation by Finite Element Analysis

Authors: Mingxing Zhou; Zunqiang Fan; Zhichao Ma; Yue Guo; Liguo Yang; Long Qian; Xingdong Sun
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

In reality, nanoindentation test is an efficient technique for probing the mechanical properties of biological tissue that soaked in the liquid media to keep the bioactivity. However, the effects of flotage imposed on the indenter will lead to inaccuracy when calculating mechanical properties (for instance, elastic modulus and hardness) by using depth-sensing nanoindentation. In this paper, the effects of flotage on the nanoindentation results of cortical bone were investigated by finite element analysis (FEA) simulation. Comparisons of nanoindentation simulation results of bone samples with and without being soaked in the liquid media were carried out. Conclusions show that the difference of load-displacement curves in the case of soaking sample and without soaking sample conditions varies widely based on the change of indentation depth. In other words, the nanoindentation measurements in liquid media will cause significant error in the calculated Young’s modules and hardness due to the flotage. By taking into account the effect of flotage, these errors are particularly important to the accurate biomechanics characterization of biological samples.


---

## Body

## 1. Introduction

It is now appreciated that nanoindentation technique has been successfully applied to characterize elastic and plastic behaviors of biological tissue such as Young’s modulus, hardness, and viscoelasticity, when combined with the in situ observation of electron microscope [1–3]. Recently, Young et al. [4] probed the viscoelasticity of a spider’s biological vibration filter by utilizing an atomic force microscopy (AFM) technique, and they investigated the elastic modulus over a temperature range of 15–40°C at various loading frequencies. Yoo et al. [5] applied an indentation technique to characterize viscoelastic properties of small ocular and orbital tissue specimens, and they developed quantitative models to predict the wide range response of the tissue’s biomechanics based on the experimental measures. To maintain the freshness and prevent dehydration, all of experimental biological tissues should be soaked in the liquid media during the indentation. However, the method to keep the bioactivity of biological sample is seldom mentioned in the above-mentioned experiments. And this will lead to the inaccurate characterization of mechanical properties. Nevertheless, Selby et al. [6] conducted the indentation experiment of the contact lens soaked in phosphate buffered saline during the test. They investigated the influence of hydrogel film thickness on the measured elastic modulus. Li et al. [7] utilized the indentation method to characterize the elastic and viscoelastic mechanical properties of brain tissues, and they immersed the tissue sample in the artificial cerebrospinal fluid to maintain its bioactivity during the whole experimental process. However, when conducting the nanoindentation test in the liquid, the indenter will inevitably bear the flotage, and this will result in the inaccurate collection of reaction force about the indenter. Specifically, the effects of flotage in the experimental microenvironment were not investigated in detail in their work. Regarding the fact that precise evaluation of mechanical properties of biological tissue in immersion indentation is based on fitting unloading stage of load-displacement curves and that these curves are determined by accurate collection of contact force between the indenter and sample, the influence of flotage caused by the fluid media should be paid more attention to in biomechanical tests.This paper focuses on the quantitative influence of flotage on the elastic-plastic indentation mechanical properties of a cortical bone specimen. Three kinds of typical indenters were selected for comparative studies of the load-displacement(P-h) curves in traditional nanoindentation and immersion indentation. When several kinds of indenters penetrate into the liquid media, different magnitudes of flotage will be applied to the corresponding indenters resulting from their geometric diversities. Moreover, percent errors of the specimen’s elastic modulus and hardness calculated based on the obtained load-displacement (P-h) curves were listed in very intuitive terms to evaluate the effect of flotage. The research may have certain reference values on indentation of bone, brain, and other biological materials in wet environment.


## 2. Theory of Traditional Indentation and Immersion Indentation

### 2.1. Traditional Indentation

It is well known that the method of nanoindentation analysis to determine Young’s modulusE and hardness H proposed by Oliver and Pharr is applicable for bulk specimens [8]. The formula of Oliver and Pharr is established based on the following assumptions [9]:(i)
The material of specimen is homogenous during the deformation process.(ii)
The material of specimen undergoes elastic-plastic deformation on the loading process and only elastic recovery on the unloading process.(iii)
The creep and viscoelastic deformation of the material should be negligible.Obviously, the method of calculating elastic modulus and hardness is different from traditional method, which can be written as follows [8, 9]: (1)S=dpdh(2)hc=hmax-ε×pmaxS(3)Ac=fhc(4)H=PmaxAc(5)Er=sqrtπ×S2β×sqrtAc(6)1Er=1-vs2Es+1-vi2Ei.In (1) to (6),S is the contact stiffness by fitting the unloading force-displacement curves [8, 9]. Pmax is the maximum indentation load; Ac is the projected contact area under the peak indentation depth. hc is contact depth that is used to calculate the real contact area and is determined by Sneddon’s equation [8]. H is the hardness of specimen, β is a correction factor that depends on the geometry of indenter, and ε is a geometry constant that equals 0.75 for a Berkovich indenter [10].Obviously, the relation betweenE and H can be determined by the projected contact area Ac, contact stiffness S, Young’s modulus and Poisson’s ratio of indenter (Ei and νi), contact depth hc, reduced modulus Er, and indentation load-displacement (P-h) curve. Regarding the traditional indentation testing on a bone specimen, the unloading stage of load-displacement (P-h) curve directly responds to the elastic-plastic deformation behavior and determines the mechanical parameters of materials [11]. According to (1) and (4), H is directly determined by the ratio of Pmax to Ac, and S is defined by the slope of the initial unloading stage of P-h curve (relevant with hmax and residual depth hf) [8]. Furthermore, on the basis of (2) and (3), hc is determined by maximum indentation depth hmax andS[12], and Ac is proportional to the square of hc [9]. According to (5) and (6), the elastic modulus and hardness of the specimen will be easily obtained on the condition that the contact stiffness S and the projected contact area Ac are known.


### 2.2. Immersion Indentation

A general sketch map of immersion indentation experiment is shown in Figure1(a). Cortical bone specimen is rigidly fixed on the glass slide that soaked in the liquid media. Assuming that the bone specimen itself will not be influenced by the flotage, only the indenter will bear flotage during the whole experiment process. When conducting the immersion indentation, the indenter moves downward bit by bit, and then the indenter will bear flotage exactly at the time when it touches the liquid level. And this flotage will increase along with the increases of immersion depth h during the whole immersion state. Therefore, the maximum load Pmax in immersion indentation might be higher than that in traditional indentation due to the inaccurate collection of reaction force on the indenter. Accordingly, the loading and unloading stages of P-h curve in immersion indentation will both exhibit different extent of deviation from the traditional indentation P-h curve. Thus, it will be of great significance to be certain of the effect of flotage on these two kinds P-h curves by a numerical method. Figure 1(b) shows a two-dimensional axisymmetric immersion indentation model, whose Berkovich indenter tip is exactly in touch with the specimen. And h is the measured immersion depth, hflo is the submergence depth of specimen, and hind is the tip height of indenter. The clamping end of the indenter is equivalent to a cylinder and its radius is Ra. The relationship between the submergence depth of specimen hflo and the tip height of indenter hind will be clearly clarified here.Figure 1
(a) Sketch map of immersion indentation experiment and (b) two-dimensional axisymmetric experiment model of nanoindentation.
(a)
(b)Obviously, the measured immersion depthh can be considered as the sum of the maximum indentation depth hmax and the submergence depth hflo: namely, h=hmax+hflo. Moreover, the flotage applied to the indenter will increase along with the increase of immersion depth h; thus, the measured immersion depth h can be classified into two different conditions; namely, hflo≤hind and hflo>hind.Assuming that the fluid keeps relatively static in the whole indentation process and the liquid media’s surface tension adhered to the indenter can be neglected, thus the indenter will only bear the flotage in liquid media. It is clear that the magnitude of this flotage depends on the geometry parameters of indenters, and Table1 shows some geometry parameters about three indenters used in this paper [13, 14]. Here, φ is the equivalent semicone angle, α is the included angel between indenter’s geometry center line and faced plane, and it will be utilized to calculate the indenters’ geometric volume that determines the flotage.Table 1
Geometric parameters of different indenters.
Berkovich indenter
Cubic indenter
Spherical indenter
A
c
24.56hc2
2.5981hc2
π2Rhc-hc2
φ
70.32°
42.28°
α
65.3°
35.2644°Regarding the Berkovich indenter, when the submergence depthhflo is smaller than the tip height of indenter hind, namely, hflo≤hind, the formula of flotage can be written as follows:(7)Pflo=3ρgh3tan2⁡αh≤hindPflo=3ρghind3tan2⁡α+ρgπRa2h-hindh>hind.Here, Pflo is the flotage applied to the indenter, ρ is the density of liquid media, g is the ratio of gravity and mass, and h is the measured immersion depth. On the other hand, the formula of flotage can be written as follows when the condition is hflo>hind: (8)Pflo=3ρghind3tan2⁡α+ρgπRa2h-hind.Similarly, the flotage applied to other spherical indenter and cubic indenter can be deduced by using trigonometric functions and Archimedes flotage principle. Regarding the spherical indenter, when the submergence depthhflo is smaller than the tip height of indenter hind, namely, hflo≤hind, the formula of flotage can be written as (9)Pflo=ρgπh2Ra-h3h≤hind.Pflo=23ρgπRa3+ρgπRa2h-hindh>hindOn the other hand, when the condition is hflo>hind, the formula of flotage can be written as (10)Pflo=3ρghind3tan2⁡α+ρgπRa2h-hind.Regarding the cubic indenter, when the submergence depthhflo is smaller than the tip height of indenter hind, namely, hflo≤hind, the formula of flotage can be written as (11)Pflo=3ρgh3tan2⁡αh≤hindPflo=3ρghind3tan2⁡α+ρgπRa2h-hindh>hind.On the other hand, when the condition is hflo>hind, the formula of flotage can be written as(12)Pflo=3ρghind3tan2⁡α+ρgπRa2h-hind.


## 2.1. Traditional Indentation

It is well known that the method of nanoindentation analysis to determine Young’s modulusE and hardness H proposed by Oliver and Pharr is applicable for bulk specimens [8]. The formula of Oliver and Pharr is established based on the following assumptions [9]:(i)
The material of specimen is homogenous during the deformation process.(ii)
The material of specimen undergoes elastic-plastic deformation on the loading process and only elastic recovery on the unloading process.(iii)
The creep and viscoelastic deformation of the material should be negligible.Obviously, the method of calculating elastic modulus and hardness is different from traditional method, which can be written as follows [8, 9]: (1)S=dpdh(2)hc=hmax-ε×pmaxS(3)Ac=fhc(4)H=PmaxAc(5)Er=sqrtπ×S2β×sqrtAc(6)1Er=1-vs2Es+1-vi2Ei.In (1) to (6),S is the contact stiffness by fitting the unloading force-displacement curves [8, 9]. Pmax is the maximum indentation load; Ac is the projected contact area under the peak indentation depth. hc is contact depth that is used to calculate the real contact area and is determined by Sneddon’s equation [8]. H is the hardness of specimen, β is a correction factor that depends on the geometry of indenter, and ε is a geometry constant that equals 0.75 for a Berkovich indenter [10].Obviously, the relation betweenE and H can be determined by the projected contact area Ac, contact stiffness S, Young’s modulus and Poisson’s ratio of indenter (Ei and νi), contact depth hc, reduced modulus Er, and indentation load-displacement (P-h) curve. Regarding the traditional indentation testing on a bone specimen, the unloading stage of load-displacement (P-h) curve directly responds to the elastic-plastic deformation behavior and determines the mechanical parameters of materials [11]. According to (1) and (4), H is directly determined by the ratio of Pmax to Ac, and S is defined by the slope of the initial unloading stage of P-h curve (relevant with hmax and residual depth hf) [8]. Furthermore, on the basis of (2) and (3), hc is determined by maximum indentation depth hmax andS[12], and Ac is proportional to the square of hc [9]. According to (5) and (6), the elastic modulus and hardness of the specimen will be easily obtained on the condition that the contact stiffness S and the projected contact area Ac are known.


## 2.2. Immersion Indentation

A general sketch map of immersion indentation experiment is shown in Figure1(a). Cortical bone specimen is rigidly fixed on the glass slide that soaked in the liquid media. Assuming that the bone specimen itself will not be influenced by the flotage, only the indenter will bear flotage during the whole experiment process. When conducting the immersion indentation, the indenter moves downward bit by bit, and then the indenter will bear flotage exactly at the time when it touches the liquid level. And this flotage will increase along with the increases of immersion depth h during the whole immersion state. Therefore, the maximum load Pmax in immersion indentation might be higher than that in traditional indentation due to the inaccurate collection of reaction force on the indenter. Accordingly, the loading and unloading stages of P-h curve in immersion indentation will both exhibit different extent of deviation from the traditional indentation P-h curve. Thus, it will be of great significance to be certain of the effect of flotage on these two kinds P-h curves by a numerical method. Figure 1(b) shows a two-dimensional axisymmetric immersion indentation model, whose Berkovich indenter tip is exactly in touch with the specimen. And h is the measured immersion depth, hflo is the submergence depth of specimen, and hind is the tip height of indenter. The clamping end of the indenter is equivalent to a cylinder and its radius is Ra. The relationship between the submergence depth of specimen hflo and the tip height of indenter hind will be clearly clarified here.Figure 1
(a) Sketch map of immersion indentation experiment and (b) two-dimensional axisymmetric experiment model of nanoindentation.
(a)
(b)Obviously, the measured immersion depthh can be considered as the sum of the maximum indentation depth hmax and the submergence depth hflo: namely, h=hmax+hflo. Moreover, the flotage applied to the indenter will increase along with the increase of immersion depth h; thus, the measured immersion depth h can be classified into two different conditions; namely, hflo≤hind and hflo>hind.Assuming that the fluid keeps relatively static in the whole indentation process and the liquid media’s surface tension adhered to the indenter can be neglected, thus the indenter will only bear the flotage in liquid media. It is clear that the magnitude of this flotage depends on the geometry parameters of indenters, and Table1 shows some geometry parameters about three indenters used in this paper [13, 14]. Here, φ is the equivalent semicone angle, α is the included angel between indenter’s geometry center line and faced plane, and it will be utilized to calculate the indenters’ geometric volume that determines the flotage.Table 1
Geometric parameters of different indenters.
Berkovich indenter
Cubic indenter
Spherical indenter
A
c
24.56hc2
2.5981hc2
π2Rhc-hc2
φ
70.32°
42.28°
α
65.3°
35.2644°Regarding the Berkovich indenter, when the submergence depthhflo is smaller than the tip height of indenter hind, namely, hflo≤hind, the formula of flotage can be written as follows:(7)Pflo=3ρgh3tan2⁡αh≤hindPflo=3ρghind3tan2⁡α+ρgπRa2h-hindh>hind.Here, Pflo is the flotage applied to the indenter, ρ is the density of liquid media, g is the ratio of gravity and mass, and h is the measured immersion depth. On the other hand, the formula of flotage can be written as follows when the condition is hflo>hind: (8)Pflo=3ρghind3tan2⁡α+ρgπRa2h-hind.Similarly, the flotage applied to other spherical indenter and cubic indenter can be deduced by using trigonometric functions and Archimedes flotage principle. Regarding the spherical indenter, when the submergence depthhflo is smaller than the tip height of indenter hind, namely, hflo≤hind, the formula of flotage can be written as (9)Pflo=ρgπh2Ra-h3h≤hind.Pflo=23ρgπRa3+ρgπRa2h-hindh>hindOn the other hand, when the condition is hflo>hind, the formula of flotage can be written as (10)Pflo=3ρghind3tan2⁡α+ρgπRa2h-hind.Regarding the cubic indenter, when the submergence depthhflo is smaller than the tip height of indenter hind, namely, hflo≤hind, the formula of flotage can be written as (11)Pflo=3ρgh3tan2⁡αh≤hindPflo=3ρghind3tan2⁡α+ρgπRa2h-hindh>hind.On the other hand, when the condition is hflo>hind, the formula of flotage can be written as(12)Pflo=3ρghind3tan2⁡α+ρgπRa2h-hind.


## 3. FEA Simulation Results and Discussion

### 3.1. Simulation Model

In order to investigate the effects of flotage on the results of immersion indentation, finite element analysis (FEA) was carried out by using the commercial software ABAQUS 6.11 version [15]. Figure 2 showed a two-dimensional simulation model of nanoindentation, whose bottom of the specimen was frozen and its left side was applied with axisymmetric boundary conditions. It can be confirmed from previous FEA results that the tip of the indenter can be equivalent to a half cone angle, and the equivalent semicone angle of Berkovich indenter tip was 70.32° [13, 14]. The specimen was meshed with four node bilinear axisymmetric reduced integration elements (CAX4R element type), and the whole mesh size was 0.3 μm. No friction was considered between the indenter and specimen. Due to the large deformations involved in the simulation, adaptive meshing was used. The indenter had only one degree of freedom and was applied with displacement boundary to accomplish the loading and unloading steps. Compared with the sample, the hardness of diamond indenter was extremely large, so the indenter was treated as a rigid body in the simulation. The total size of 2D axisymmetric meshed simulation model was 15 × 15 μm and the indenter height was set to be 1.5 μm here. A homogenous isotropic cortical bone was chosen as the simulation sample and it was treated as ideal elastic-plastic body. The detailed mechanical properties were given in Table 2 [16].Table 2
Elastic and plastic properties of cortical bone.
Material
Young’s modulus,E (GPa)
Poisson’s, ratio, ν
Yieldstrength
σ
s (GPa)
Cortical bone
13.56
0.3
0.301Figure 2
The two-dimensional meshed axisymmetric simulation model of nanoindentation.The FEA indentation process of immersion indentation and traditional indentation was divided into two steps of loading and unloading. Both loading and unloading steps were all controlled by the displacement method with the maximum indentation depth of indenterhmax being 200, 400, 600, 800, and 1000 nm. Regarding the immersion indentation, indenters will bear some flotage increasing with the immersion depth h, and therefore the indenter will suffer an equivalent flotage load when it immerses into the liquid medium in the simulation. Assuming that the homogeneous liquid media are 3.5% normal saline and the density of this normal saline is 1.03 g/cm3, when the submergence depth of specimen hflo is 2.28 mm, it can be calculated from (7) that the flotage of Berkovich indenter Pflo is 1 mN. Similarly, regarding the cubic indenter whose geometry standard is identical with Berkovich indenter, the flotage of it is 0.106 mN. According to (9), assuming that the radius of spherical indenter Ra is 3 mm, the flotage of this indenter will be 0.377 mN.


### 3.2. Simulation Results

Figure3 illustrates corresponding P-h curves of three kinds of indenters in the immersion indentation and traditional indentation. It can be clearly seen that the loads at the maximum indentation depth hmax in immersion indentation are higher than that in traditional indentation, and the envelope area of P-h curves in immersion indentation is apparently larger than that in traditional indentation. In addition, the loading and unloading stages of P-h curves both exhibit different extent of deviations from traditional indentation P-h curves. Specifically, when the maximum indentation depth hmax equals 200 nm, the deviations of P-h curves between immersion indentation and traditional indentation become most obvious no matter what kind of indenter is employed in the simulation. That is to say, these deviations will cause inaccuracy of sample’s elastic modulus and hardness to a certain extent.Figure 3
Load-displacement curves of immersion indentation and traditional indentation with different indenters, (a) cubic indenter, (b) Berkovich indenter, and (c) spherical indenter.
(a)
(b)
(c)To have a quantitative explanation for the effects of flotage on the results of immersion indentation, the percent errors of elastic modulus and hardness are calculated based on the following equation:(13)ErrorE%=Eflo-EsEs×100%ErrorH%=Hflo-HsHs×100%.Here,Eflo is the calculated elastic modulus of immersion indentation and Es is the calculated elastic modulus of traditional indentation. Similarly, Hflo is the calculated hardness of immersion indentation and Hs is the calculated hardness of traditional indentation.As drawn in Figure4, the percent errors of elastic modulus and hardness for three kinds of indenters in immersion indentation are compared with that of traditional indentation. It is obvious that the percent errors of elastic modulus ErrorE and percent errors of hardness ErrorH will decrease with the increase of the maximum indentation depth hmax. Specifically, when the maximum indentation depth hmax equals 200 nm, the percent errors of elastic modulus ErrorE exceed 100%, regardless of the kinds of indenters. Meanwhile, the percent errors of hardness ErrorH are about 120.6%, 82.7%, and 56.7%, which are obtained from Berkovich indenter, cubic indenter, and spherical indenter, respectively.Figure 4
(a) Percent errors of elastic modulus for different indenters with immersion indentation compared with traditional indentation. (b) Percent errors of hardness for different indenters with immersion indentation compared with traditional indentation.
(a)
(b)Similarly, when the maximum indentation depthhmax increases to 1000 nm, the ErrorE values decrease to 17.5%, 4.4%, and 10.9%, while the ErrorH values decrease to 5.8%, 5.1%, and 7.2% when using Berkovich, cubic, and spherical indenters, respectively. Finally, the percent errors of elastic modulus for cubic indenter have always been the smallest one when compared with other indenters at the same indentation depth hmax.


## 3.1. Simulation Model

In order to investigate the effects of flotage on the results of immersion indentation, finite element analysis (FEA) was carried out by using the commercial software ABAQUS 6.11 version [15]. Figure 2 showed a two-dimensional simulation model of nanoindentation, whose bottom of the specimen was frozen and its left side was applied with axisymmetric boundary conditions. It can be confirmed from previous FEA results that the tip of the indenter can be equivalent to a half cone angle, and the equivalent semicone angle of Berkovich indenter tip was 70.32° [13, 14]. The specimen was meshed with four node bilinear axisymmetric reduced integration elements (CAX4R element type), and the whole mesh size was 0.3 μm. No friction was considered between the indenter and specimen. Due to the large deformations involved in the simulation, adaptive meshing was used. The indenter had only one degree of freedom and was applied with displacement boundary to accomplish the loading and unloading steps. Compared with the sample, the hardness of diamond indenter was extremely large, so the indenter was treated as a rigid body in the simulation. The total size of 2D axisymmetric meshed simulation model was 15 × 15 μm and the indenter height was set to be 1.5 μm here. A homogenous isotropic cortical bone was chosen as the simulation sample and it was treated as ideal elastic-plastic body. The detailed mechanical properties were given in Table 2 [16].Table 2
Elastic and plastic properties of cortical bone.
Material
Young’s modulus,E (GPa)
Poisson’s, ratio, ν
Yieldstrength
σ
s (GPa)
Cortical bone
13.56
0.3
0.301Figure 2
The two-dimensional meshed axisymmetric simulation model of nanoindentation.The FEA indentation process of immersion indentation and traditional indentation was divided into two steps of loading and unloading. Both loading and unloading steps were all controlled by the displacement method with the maximum indentation depth of indenterhmax being 200, 400, 600, 800, and 1000 nm. Regarding the immersion indentation, indenters will bear some flotage increasing with the immersion depth h, and therefore the indenter will suffer an equivalent flotage load when it immerses into the liquid medium in the simulation. Assuming that the homogeneous liquid media are 3.5% normal saline and the density of this normal saline is 1.03 g/cm3, when the submergence depth of specimen hflo is 2.28 mm, it can be calculated from (7) that the flotage of Berkovich indenter Pflo is 1 mN. Similarly, regarding the cubic indenter whose geometry standard is identical with Berkovich indenter, the flotage of it is 0.106 mN. According to (9), assuming that the radius of spherical indenter Ra is 3 mm, the flotage of this indenter will be 0.377 mN.


## 3.2. Simulation Results

Figure3 illustrates corresponding P-h curves of three kinds of indenters in the immersion indentation and traditional indentation. It can be clearly seen that the loads at the maximum indentation depth hmax in immersion indentation are higher than that in traditional indentation, and the envelope area of P-h curves in immersion indentation is apparently larger than that in traditional indentation. In addition, the loading and unloading stages of P-h curves both exhibit different extent of deviations from traditional indentation P-h curves. Specifically, when the maximum indentation depth hmax equals 200 nm, the deviations of P-h curves between immersion indentation and traditional indentation become most obvious no matter what kind of indenter is employed in the simulation. That is to say, these deviations will cause inaccuracy of sample’s elastic modulus and hardness to a certain extent.Figure 3
Load-displacement curves of immersion indentation and traditional indentation with different indenters, (a) cubic indenter, (b) Berkovich indenter, and (c) spherical indenter.
(a)
(b)
(c)To have a quantitative explanation for the effects of flotage on the results of immersion indentation, the percent errors of elastic modulus and hardness are calculated based on the following equation:(13)ErrorE%=Eflo-EsEs×100%ErrorH%=Hflo-HsHs×100%.Here,Eflo is the calculated elastic modulus of immersion indentation and Es is the calculated elastic modulus of traditional indentation. Similarly, Hflo is the calculated hardness of immersion indentation and Hs is the calculated hardness of traditional indentation.As drawn in Figure4, the percent errors of elastic modulus and hardness for three kinds of indenters in immersion indentation are compared with that of traditional indentation. It is obvious that the percent errors of elastic modulus ErrorE and percent errors of hardness ErrorH will decrease with the increase of the maximum indentation depth hmax. Specifically, when the maximum indentation depth hmax equals 200 nm, the percent errors of elastic modulus ErrorE exceed 100%, regardless of the kinds of indenters. Meanwhile, the percent errors of hardness ErrorH are about 120.6%, 82.7%, and 56.7%, which are obtained from Berkovich indenter, cubic indenter, and spherical indenter, respectively.Figure 4
(a) Percent errors of elastic modulus for different indenters with immersion indentation compared with traditional indentation. (b) Percent errors of hardness for different indenters with immersion indentation compared with traditional indentation.
(a)
(b)Similarly, when the maximum indentation depthhmax increases to 1000 nm, the ErrorE values decrease to 17.5%, 4.4%, and 10.9%, while the ErrorH values decrease to 5.8%, 5.1%, and 7.2% when using Berkovich, cubic, and spherical indenters, respectively. Finally, the percent errors of elastic modulus for cubic indenter have always been the smallest one when compared with other indenters at the same indentation depth hmax.


## 4. Conclusions

Regarding the immersion indentation whose specimen is soaked in the liquid media, the effects of flotage on the indentation results used to be neglected, and this will cause the properties such as elastic modulus and hardness of the sample to be somehow inaccurate. In this paper, the effects of flotage on the results of cortical bone nanoindentation were investigated by finite element analysis simulation. Conclusions are drawn as follows:(1)
The comparison studies inP-h curves of immersion indentation and traditional indentation indicate that the flotage could to some extent affect the envelope area of the P-h curves. When under the same indentation depth, the envelope area of immersion indentation will be larger than that of traditional indentation.(2)
From the calculated percent errors of elastic modulus and hardness, it is known that when the maximum indentation depth ranges from 200 nm to 1000 nm, the percent errors of elastic modulus for cubic indenter is the smallest and the Berkovich indenter is the largest one all the time.(3)
With the decrease of maximum indentation depth, the percent errors values of elastic modulus and hardness become larger on the contrary.



---

*Source: 4569351-2017-09-11.xml*