The present invention is related to the indentation test methods for characterizing the various types of mechanical properties as well as the surface adhesion energy of elastic/elastoplastic/viscoelastic materials.
The conventional indentation test is a testing technique that evaluates various mechanical properties such as the hardness and the elastic modulus through measuring the size of residual impression that formed by pressing an indenter onto a surface of test specimen.
When an indenter is pressed onto the surface of a perfectly elastic body, its elastic modulus can be evaluated from the indentation-induced elastic deformation, while due to the plastic flow in elastoplastic indentation contact, the residual contact impression is formed after unloading. The size of this residual impression combined with the indentation load leads to the contact hardness. The elastic modulus and the contact hardness are the representative examples of several mechanical characteristics that are determined in indentation testing.
In ductile metals, the size of residual impression nearly coincides with the impression formed at the maximum indentation load due to the significant plastic flow. In ceramics and/or organic polymers with their significant elastoplastic and/or viscoelastic deformation and flow, on the other hand, their residual impressions after unloading are always smaller than those formed at the maximum indentation load due to the elastic recovery of impression induced in its unloading process.
Therefore, in order to quantitatively characterize the mechanical properties of a test specimen, it is necessary to make in-situ measurements of indentation contact, including the measurements of indentation load and penetration depth along with the indentation-induced contact area both in the loading and in the subsequent unloading processes.
In the conventional instrumented nanoindentation test method, it is only capable of measuring the penetration depth (h) under the applied indentation load (P). Outside the indentation contact area, the free-surface of an elastic body always sinks-in, while it piles-up for elastoplastic body. Without the information on these sink-in/pile-up contact profiles, therefore, it is by no means capable of conducting the quantitative characterization of mechanical properties of the test specimen indented. In order to circumvent this fatal issue of sink-in/pile-up contact profiles, undesirable approximations and assumptions have always been made in the conventional instrumented nanoindentation testing to estimate the indentation contact area (A) formed under the applied load (P). In other words, without the information on the indentation contact area A under the applied load P, no one determines in a quantitative manner the mechanical characteristics of any engineering materials.
The instrumented indentation microscope is the test system to overcome the difficulties included in the conventional instrumented nanoindentation test. The instrumented indentation microscope is capable of not only measuring the penetration depth (h), but also optically determining the projected contact area (A) in an in-situ manner under the applied load of indentation, leading to a quantitative characterization of mechanical characteristics without any of undesirable approximations/assumptions that are required in the conventional nanoindentation testing. The instrumented indentation microscope can therefore measure the projected contact area both for “sink-in” and “pile-up” contact surface profiles (refer to Patent Literature 1, Patent Literature 2, Patent Literature 3, Patent Literature 4, Patent Literature 5, Non Patent Literature 1, Non Patent Literature 2, Non Patent Literature 3, and Non Patent Literature 4).
However, in characterizing the mechanical properties of a specimen by indentation test method, it is requisite to take into account the effect of surface adhesion of the specimen on the indentation contact behavior. In the cases of so-called the hard materials having rather large elastic modus (the plane strain Young's modulus, E′≥100 GPa) such as the engineering materials of metals and ceramics, due to their large elastic moduli, the surface adhesion, even if there exists, makes a rather minor effect on their contact behavior. In most cases, therefore, we can neglect the surface adhesion in their indentation contact testing.
On the other hand, the elastic moduli of soft materials including organic polymers, biomaterials, and microbiological mediums are rather small, and fall in the range of E′≈1 Pa-100 MPa, where the surface adhesion plays an essential role in their mechanical characteristics and functions, leading to critical difficulties in quantitatively characterizing the mechanical properties through using the conventional nanoindentation test systems.
The present invention has been made in view of such circumstances and aims to provide a test technique and analysis for simultaneously characterizing the adhesion energy as well as the various mechanical properties of soft material through taking into account the adhesion force at the contact interface between the material tested and the indenter.
The present invention provides the following technical method and procedure for solving the problems and issues described above.
A test method for characterizing the mechanical properties including the surface adhesion energy γ on the basis of the experimentally derived P-A relationship,
where P means the indentation load under the penetration depth h of an indenter pressed onto a test specimen with surface adhesion, and
A means the contact area of indentation at the contact radius a under the applied load of P.
Based the present invention, it is possible to provide an indentation test method for simultaneously characterizing in a quantitative manner the adhesion energy as well as the various mechanical properties of soft materials with surface adhesion.
Hereinafter, the embodiments of the present invention will be described.
An axisym metric indenter with an arbitrary shape is penetrated into an elastic body having the elastic modulus (the plane strain Young's modulus) of E′ under the indentation load P, resulting in the penetration depth of h, and the contact radius of a at this time. It is assumed that the elastic modulus of the indenter is sufficiently larger than that of the test specimen (the E′-value of a diamond indenter is about 1000 GPa, and about 410 GPa of a sapphire indenter, by way of example). In the case where the elastic body has a surface adhesion, the indenter is withdrawn to the elastic body. That is, the surface adhesion induces “negative” contact pressure.
This fact implies that the indentation load P at the contact radius a will be smaller than that of the elastic body without surface adhesion. The JKR theory models the surface adhesion as the negative contact pressure acting on the flat-ended cylindrical punch with the radius a. In the present theoretical considerations, the JKR theory that has only been applied to spherical indentation will be extended to the conventional pyramid/cone indentations including the Vickers, Berkovich as well as the Rockwell indentations.
The contact pressure distribution p(r) of a flat-ended cylindrical indenter having a radius a is given by the following Formula.
In the above equation, the suffix F indicates a flat-ended cylindrical indenter (Flat punch), and r indicates the radius from the indentation axis (z-axis).
In the case where the surface adhesion is expressed using Formula (1), the coefficient pF has a “negative value” (pF<0) because the indenter is withdrawn to the surface of elastic body. Similarly, the indentation pressure distribution ps(r) of the spherical indenter can be expressed using the following Formula.
Furthermore, the indentation pressure distribution PC(r) of the conical indenter can be expressed using the following Formula.
[Equation 3]
p
C(r)=C cos h−1(a/r); 0≤r<a (3)
In Formulas (2) and (3), the suffixes S and C of the coefficients ps and pc represent a spherical indenter (Sphere) and a conical indenter (Cone), respectively.
Therefore, in the indentation to a perfectly elastic body having a surface adhesion, the contact pressure distribution generated immediately beneath the indenter can be expressed by the following Formula in the case of a spherical indenter by superimposing the pressure distributions of Formulas (1) and (2).
In the case of a conical indenter, the distribution can be expressed by the following Formulas by superimposing the pressure distributions of Formulas (1) and (3).
Meanwhile, the contact surface profile uz(r) (0≤r≤a) beneath a spherical indenter having the radius of R can be expressed by the following geometric relation.
Similarly, the contact surface profile uz(r) (0≤r≤a) induced beneath a conical indenter having the inclined face-angle of β can be expressed by the following geometric expression.
[Equation 7]
u
Z(r)=h−r tan β (7)
Through combining the contact pressure distributions p(r) of Formulas (1) to (5) and the contact surface profiles uz(r) of Formulas (6) or (7), the following expression is finally obtained for spherical indentation:
Notice the algebraic identity of Formula (8) as to the variable r, then the following relations are obtained;
By further using Formula (4), the indentation load (P) of spherical indentation is expressed with the following formula:
Similar mathematical operations conducted for spherical indentation finally lead to the key expressions of the conical indentation as follows;
Through these mathematical procedures, the coefficients of contact pressure distribution, ps and pc are related to the elastic modulus E′ in Formulas (9) and (13).
Meanwhile, it is impossible to determine the coefficient of contact pressure distribution pF of flat-ended cylindrical indenter in terms of the surface adhesion via the preceding mathematical operation. To overcome this difficulty, therefore, noticing the fact that the adhesion force γ (N/m) is equivalent to the adhesion energy γ (J/m2), an energy-based consideration will be made in what follows for determining pF.
In the first step of indentation contact process, suppose an axisym metric indenter pressed onto an elastic body without surface adhesion to the indentation load P1 as depicted in
[Equation A1]
U1=∫0h
The application of this integral to the spherical indentation results in
and to the conical indentation leads to
The elastic strain energy U1 stored in the body is given by the area OABO depicted in
In the subsequent second step of indentation contact process, at the point A (P1, h1, a1) in
This mechanical process implies that the indenter is progressively pulled to the contact surface, resulting in unloading, as shown in
The indentation contact state at the point C (P2, h2, a1) is the equilibrium state of the elastic body having the surface energy of γ. This mechanical process along the line AC will be equivalent to the unloading process of a flat-punch with the radius a1. The total energy released through this unloading process along the line AC is denoted by U2(<0) and given by the area ABDCA (=−U2) in
As described above, the unloading process along the line AC is equivalent to the unloading process of the flat-ended cylindrical punch with the radius a1, The P-h unloading path along the line AC is expressed by the following Formulas.
where pF1 (<0) is the coefficient of contact pressure distribution of flat-ended cylindrical punch with radius a1, having a negative value due to the surface adhesion.
The load P1 at point A shown in
P
1=(2πa13/3)pS1 [Equation A2]
The load P1 for conical indentation can be expressed by the following Formula.
P1=πa12pC1 [Equation A3]
On the other hand, the released energy U2(<0) associated with the incremental surface adhesion is given by
Therefore, the following expressions are finally obtained as the release energy U2
for the spherical indentation;
and for conical indentation;
Therefore, the elastic strain energy UE when a spherical indenter is pressed onto the perfectly elastic body with “surface adhesion” until the contact radius becomes a, that is, the area ACDOA in
Alternatively, the energy UE can also be expressed with the following Formula by substituting Formulas (9) and (10) into the above Formula.
On the other hand, the energy UE of conical indentation can be expressed by the following Formula.
Alternatively, the energy UE can also be expressed by substituting Formulas (13) and (14) into the above Formula;
As mentioned in the preceding considerations, in the present indentation contact problem, the adhesive surface force introduces the surface energy Us which decreases when the surfaces come into intimately contact and increases when they separate. Therefore, we can write
[Equation 25]
U
S=−2γπa2 (25)
The total energy (the Gibbs free energy) UT of the present mechanical system, therefore, is given by
[Equation 26]
U
T
=U
E
+U
S (26)
In the mechanical equilibrium under a fixed depth of penetration that means none of external works applied to the system, the variation of total energy associated with incremental contact radius δ a results in
By substituting Formulas (21) to (26) into Formula (27), and by using Formulas (22) and (24), the following expression is obtained both for spherical and conical indentations.
(∂UE/∂a)h=(π2a2/E′)F2 [Equation A5]
Accordingly, the pressure distribution coefficient pF of the flat-ended cylindrical indenter is finally correlated to the adhesion energy γ by the following Formula.
Substituting the coefficients pS (Formula (9)), pC (Formula (13)), and pF (Formula (28)) obtained in the preceding considerations into Formulas (10) and (11), or to Formulas (14) and (15) leads to the key expressions of the h vs. a and the P vs. a relations (h: penetration depth, P: indentation load, a: contact radius).
The key expressions for spherical indentation are;
The key expressions for conical indentation are;
Here, A (=πa2) represents the contact area of indentation.
The coefficient of surface adhesion in Formula (32) is defined by the following Formula.
λE is referred to as the adhesion toughness that stands for the fracture toughness of interfacial delamination between the tip-of-indenter and the material indented (the suffix E indicates elastic) The physical dimension of the adhesion toughness λE is [Pa·m1/2], being the same as the mode-I fracture toughness KIc;
[Equation 34]
K
Ic(≡√{square root over (2γE′)}) (34)
By substituting the adhesion energy γ=0 into Formulas (29) to (32), those Formulas are naturally reduced to the well-known indentation contact mechanics relations of the perfectly elastic body without surface adhesion.
The elastoplastic body with surface adhesion will be examined in what follows.
Unlike the perfectly elastic body examined in the preceding considerations, an elastoplastic body leads to a mechanical process in which the plastic deformation (plastic flow) in the vicinity of the surface reduces the surface adhesion.
In consideration of the effect of plastic deformation on the adhesion toughness, therefore, the JKR-based elastic theory (Formula (32)) can be extended to the elastoplastic region by the following Formula.
[Equation 35]
P=H
M
A−λ
EP
A
3/4 (35)
Here, HM represents the Meyer hardness, and the elastoplastic adhesion toughness λEP is defined by the following Formula.
The value γEP in Formula (36) represents the surface energy (surface force) under plastic flow, that is, represents the elastoplastic surface force (elastoplastic surface energy). There are no analytical solutions for the correlation between the elastoplastic adhesion toughness λEP or the elastoplastic adhesion energy γEP and the yield stress Y. There is no choice, therefore, but to derive these correlations as empirical rules through the FEA-based numerical analysis.
The preceding considerations have been made for pyramid/cone indentation with arbitrary inclined-face-angle of β, implying that all the experimental procedures combined with the analytical formulas given above are applicable to the conventional Vickers/Berkovich indentation as well as conical indentation.
The instrumented indentation microscope is designed for quantitatively determining the contact area A and the penetration depth h in an in-situ manner under the indentation load P applied to the test specimen. Accordingly, the elastic modulus E′, yield stress Y, as well as the adhesion energy γ can readily be determined through applying the experimental data to the Formulae given in the preceding considerations.
The effect of surface adhesion on the indentation contact mechanics of viscoelastic bodies will be examined in what follows.
By applying the “elastic-viscoelastic corresponding principle” to the JKR theory (Formula (32)), the following Formula is obtained in the Laplace space for the constitutive equation of a viscoelastic body with surface adhesion.
In Formula (37),
are, respectively, defined by the following equations:
Therefore, the inverse Laplace transform of Formula (37) directly results in “the constitutive equation of viscoelastic body having surface adhesion” in real space as follows;
Examine a stepwise indentation to a constant contact area A0, as an example of viscoelastic indentation test:
[Equation 43]
A(t)=u(t)A0 (43)
The function u(t) in Formula (43) stands for the Heaviside step function. The effect of surface adhesion on the indentation load relaxation will be discussed in what follows.
Formula (43) along with the relational of du(t)/dt=δ(t) (Dirac delta function) applied to Formula (42) results in;
Suppose a Maxwell viscoelastic liquid for simplicity in the following numerical procedures; the relaxation modulus E′relax(t) is given by
The numerical results of the indentation load relaxation curve (P(t) vs. t) obtained through substituting Formula (45) into Formula (44) is plotted in
In
Next, the present invention will be described in more detail with the following Examples.
In order to deepen the understanding of the effect of surface adhesion on the elastic/elastoplastic indentation contact problem, the finite element method was selected as a numerical analysis to examine the problem. In the present finite element analyses, the commercially available finite element software package of ANSYS was selected; having been recognized well in numerically analyzing the contact problems including elastic/elastoplastic/viscoelastic deformation and flow.
As an example of the finite element analysis applied to a perfectly elastic body having the elastic modulus E′=20 kPa, the numerical results of the loading-unloading P-A relations for the maximum penetration depth hmax=30 μm are shown in
The P-A loading-unloading relationship (the closed circle; ●) of the elastic body without surface adhesion (γ=0.0 J/m2) is linear and none of hysteresis is observed in its loading/unloading paths. Note that the broken line (the analytical solution of P=(E′ tan β/2)A) in
Note the fact in
The dotted line along the symbols ◯ in
The effect of plastic flow on the adhesion toughness is examined through the FEA-based numerical study.
The effect of the plastic flow on the elastoplastic adhesion energy γEP is shown in
The correlations between the elastoplastic adhesion toughness (the following formula) and the adhesion energy γ are plotted in
[Equation 46]
λEP2/E′(≡(16/√{square root over (π)})γEP) (46)
The following conclusions are obtained from
Both the elastoplastic adhesion toughness λEP and the elastoplastic adhesion energy γEP decrease with the reduction in the yield stress Y, that is, with the enhancement of plastic flow. In other words, the effect of surface adhesion on the elastoplastic indentation contact diminishes with enhancing plastic flow.
Conversely, both of γEP and λEP increase such that γEP→γ and λEP→λE together with the increase in the yield stress Y, The indentation contact behavior of these elastoplastic bodies, therefore, realizes the perfectly elastic body that is well described with the elastic JKR theory.
The aγ-value on the horizontal axis in
[Equation 47]
a
γ=γEP/γ (47)
From the above considerations and Formula (47), the shift factor aγ has a strong correlation with the plastic deformation and flow of elastoplastic body, that is, with the plastic index (PI (≡εIE′/cY)); aγ→1 for γEP→γ in the extreme of perfectly elastic deformation (PI↓0), while aγ→0 for γEP→0 in the extreme of fully plastic deformation (PI↑0).
To verify these considerations, the quantitative correlation (FEA-based numerical results) between the shift factor aγ (≡γEP/γ) and the plastic index PI (≡εIE′/cY) is shown in
The best-fitting empirical formula for duplicating the FEA-derived correlation shown in
[Equation 48]
a
γ=1/(3.5PI); PI≥0.286 (48)
Silicone rubber was selected as a perfectly elastic model specimen; both the loading and unloading P-A paths are linear and none of hysteresis is observed. A glue (3M, PN: 55) was coated on the silicone rubber to realize the surface adhesion.
All the indentation tests were conducted by the use of instrumented indentation microscope with a Berkovich indenter (diamond trihedral pyramid indenter with the inclined face-angle of β=24.75°).
Formula (32) is applied to the test results shown in
Aloe-gel was selected as an example of soft materials having an extremely low elastic modulus E′. Aloe leaves were sliced to make the test specimens with their thickness of about 3.5 mm. The translucent gel of the mesophyll was indented on the instrument indentation microscope.
Formulas (35) and (36) are applied to the loading P-A relation of aloe-gel shown in
The Meyer hardness HM was determined to be 3.0 kPa from the slope of the straight line passing through the origin of the graph shown in
The elastic modulus E′ obtained from the unloading modulus M given as an initial slope of the unloading P-A line might be significantly overestimated due to the effect of surface adhesion. In order to circumvent this issue, E′ must be determined using the unloading stiffness S given as the initial slope of the P-h unloading curve as shown in
The elastic modulus E′ of the aloe-gel was, therefore, successfully evaluated to be 19.0 kPa through applying the observed S-value in
The Meyer hardness (HM=3.0 kPa) as an elastoplastic parameter, the elastic modulus (E′=19.0 kPa) as an elastic parameter, and the yield value Y as a plastic parameter are correlated through the following formula on the basis of the “principle of the excluded volume of indentation”;
Using Formula (50), therefore, the yield stress as the plastic measure of the aloe-gel was determined to be Y=1.94 kPa through assuming the constraint fact c=2.65.
The elastoplastic adhesion toughness λEP thus obtained in the analysis made in
[Equation 51]
λEP=4√{square root over (γEPE′/π1/2)}=55 N/m3/2 (51)
finally leads to the adhesion energy (i.e. the surface energy or the force of surface tension) of the aloe-gel; γEP=17.4 mJ/m2. It will be worthy of note that this γEP-value is less than the surface energy of pure water (73 mJ/m2).
Preferred embodiments of the present invention have been described in detail through the preceding context. The present invention, however, is not limited to the specific embodiments, and various modifications and alterations may be made within the scope of the present invention described in the appended claims.
Number | Date | Country | Kind |
---|---|---|---|
2017-219816 | Nov 2017 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2018/042343 | 11/15/2018 | WO | 00 |