Measurements of Viscoelastic Functions of Polymers in The Frequency-Domain Using Nanoindentation
Measurements of Viscoelastic Functions of Polymers in The Frequency-Domain Using Nanoindentation
Measurements of Viscoelastic Functions of Polymers in The Frequency-Domain Using Nanoindentation
C Springer 2005
1. Introduction
The nanoindentation technique developed in the last two decades (Pethica and
Oliver, 1987, 1989; Oliver and Pharr, 1992; Li and Bhushan 2002) has been demon-
strated to be an effective approach to measure the material properties of very small
amounts of materials such as thin films, microstructures, and inhomogeneous ma-
terials. Loads as small as a few µN and displacements of a fraction of a nm can
be measured by the commercially available nanoindentation systems (Pethica and
Oliver, 1987; Syed et al., 1999); material properties are extracted from the load-
displacement relation. In recent years, the study on material properties at nanometer
and micrometer scales has received increasing attention due to their perceived im-
portance in applications such as microelectronics, nano-structured materials, and
346 G. HUANG ET AL.
where E and E are the uniaxial storage modulus and the loss modulus, respec-
tively, S the contact stiffness, C the damping coefficient, and A the contact area
between the indenter and the workpiece. This method was used to measure the com-
plex modulus of polyisoprene. In an effort to examine the formulas in Equation (1),
we conducted measurements using an MTS nano indenter XP system with CSM,
and used the formulas in Equation (1) to compute the complex modulus of polycar-
bonate (PC) and polymethyl methacrylate (PMMA) at 75 Hz. Data are computed
at all times at this frequency, and plotted in Figure 1, but only steady state values
represent the complex viscoelastic function. Also shown in Figure 1 are the con-
ventional data measured from dynamic mechanical analysis (DMA) (for details,
please see Sections 3 and 4) for the same batch of PC and PMMA. The uniaxial
storage modulus of PC measured by DMA at 75 Hz is 2.29 GPa. However, the stor-
age modulus computed using Equation (1) is at least 40% higher than this value,
indicating the difficulty associated with the method described in Equation (1) for
measurement of storage modulus for PC. Similar problem is evident for PMMA.
Also shown in these figures are storage modulus data measured by the proposed
method that will be discussed in Section 4.
This study is intended to develop a method to measure the complex viscoelastic
functions of time-dependent materials in the frequency-domain using nanoinden-
tation with a spherical indenter. Based on solutions for the indentation of an ax-
isymmetric indenter into a linearly elastic material, viscoelastic indentation under
a time-harmonic loading condition is analyzed using a hereditary integral operator
as proposed by Lee and Radok (1960). Formulas are derived to process the am-
plitudes of load and displacement as well as the out-of-phase angle between load
and displacement to determine the storage and loss parts of the complex compli-
ance (or modulus) function using a spherical indenter. The Lee–Radok approach is
applicable to situations where the contact area between the nanoindenter and the
workpiece does not decrease. When the condition of non-decreasing contact area is
not satisfied, the Ting approach (1966) needs to be used; the difference between the
approximation from Lee–Radok approach and the solution obtained from the Ting
approach is estimated. Dynamic nanoindentation tests were conducted on PC and
PMMA to determine the complex compliance, and results are compared with data
obtained from DMA on the identical materials to validate the method presented.
2. Theoretical Background
In this section, we present the derivation of the formulas for the computation of com-
plex compliance for a linearly viscoelastic material. The formula for the complex
modulus is also presented, as it is simply the reciprocal of the complex compliance.
Formulations will be given for a spherical indenter that will be used in experimental
verification in this study.
Figure 2 shows the geometry of a spherical indenter. We first consider the prob-
lem of a spherical indenter indenting into a half space composed of a homogenous,
isotropic and linearly elastic material. The diamond indenter is assumed to be rigid
because its Young’s modulus is at least two orders of magnitude higher than that of
polymer samples; and the material occupies the half space (z ≥ 0). The spherical
indenter has a tip radius R. Based on the Hertzian solution, under the condition
that the ratio of indentation depth to radius of indenter is not higher than 0.16
(Giannakopoulos, 2000), the relation between the applied indentation load and the
indentation displacement can be expressed by Hertz (1881) and Ling et al. (2002)
√
8 R
P= Gh 3/2 , (2)
3(1 − ν)
where G is the shear modulus, ν the Poisson’s ratio, P the applied indentation load,
and h the indentation displacement.
For the nanoindentation of a spherical indenter into a viscoelastic material, we
assume a constant Poisson’s ratio. During nanoindentation experiments in which a
relatively short time (such as ∼250 s used in this study) is involved, the Poisson’s
ratio (Lu et al., 1997) does not change significantly for some polymers in the glassy
state, such as PMMA. Thus, the use of a constant Poisson’s ratio will not cause
much error in the complex compliance data.
For a half space composed of a linearly viscoelastic material, we consider first
the case in which the contact area between the indenter tip and the work material is
MEASUREMENTS OF VISCOELASTIC FUNCTIONS OF POLYMERS 349
Figure 2. Geometry of the spherical indenter (a) Schematic diagram of the indenter. (b) A
TEM image of the spherical indenter tip.
where J (t) is the creep compliance function in shear in the time-domain, P(t) = 0
for t < 0.
Consider a sinusoidal nanoindentation load superimposed on a step loading,
represented by
where H(t) is the Heaviside unit step function, Pm is the carrier load, or main load,
P0 is the amplitude of the harmonic load. Equation (4) implies P(t) = 0 for
t < 0.
Inserting Equation (4) into Equation (3), we have
t
3(1 − ν)
h (t) =
3/2
√ Pm J (t) + ωP0 J (t − θ ) cos ωθ dθ . (5)
8 R 0
√
The contact radius is a(t) = Rh(t) for h(t) R. Considering that the complex
compliance is defined after the harmonic response has reached a steady state (or
equivalently t → ∞), we have
3(1 − ν)
h 3/2 (t) = √ {Pm J (t) + P0 [J (ω) sin ωt − J (ω) cos ωt]} (6)
8 R
∞ ∞
where J (ω) = ω 0 J (t) sin ωt dt and J (ω) = −ω 0 J (t) cos ωt dt are stor-
age compliance and loss compliance in shear, respectively. Note that the complex
compliance in shear is J ∗ (ω) = J (ω) − i J (ω).
On the other hand, the total displacement as output from a nanoindenter is
expressed as
where o(h 0 ) represents high-order terms of h 0 , and they are negligible under
condition of h 0 h m (t). Comparing Equation (6) with Equation (8), we find
that
3(1 − ν)
m (t) =
h 3/2 √ Pm J (t), (9)
8 R
√ 1/2 √ 1/2
4 R h m (t)h 0 4 R h m (t)h 0
J (ω) = cos δ, and J (ω) = sin δ.
1−ν P0 1−ν P0
(10)
MEASUREMENTS OF VISCOELASTIC FUNCTIONS OF POLYMERS 351
Under the loading condition in which a small sinusoidal load is superimposed upon
a ramp loading, i.e.,
Comparing Equation (12) with Equation (8), the same formulas as in Equation (10)
for the complex compliance can be derived for a small oscillatory load superimposed
upon a ramp loading.
Equation (10) is used to determine the complex compliance in shear from nanoin-
dentation. The uniaxial complex compliance, D(ω), can be computed from
J (ω) − i J (ω)
D ∗ (ω) = D (ω) − i D (ω) = , (13)
2[1 + ν ∗ (ω)]
where D (ω) and D (ω) are uniaxial storage modulus and loss modulus, respec-
tively, ν ∗ (ω) the complex Poisson’s ratio. Assuming that ν ∗ (ω) is a constant during
short-time nanoindentation tests for glassy polymers, from Equations (10) and (13),
D (ω) and D (ω) can be computed by
√ 1/2 √ 1/2
2 R h m (t)h 0 2 R h m (t)h 0
D (ω) = cos δ, and D (ω) = sin δ.
1 − ν2 P0 1 − ν2 P0
(14)
approach when the contact area is non-decreasing. The Ting approach, however, is
necessary in the case where decreasing contact area occurs in nanoindentation.
We next provide a condition under which non-decreasing contact area is main-
tained so that the solution derived
√ from the Lee–Radok approach is valid. In terms
of Equation (8) and a(t) = Rh(t), for a small harmonic loading superimposed
on a ramp loading, if v0 ≥ P0 ω, as seen from Equation (11) for a harmonic
loading superimposed on a ramp loading, the non-decreasing load leads to the non-
decreasing contact area in the entire indentation history. For a harmonic loading
superimposed on a step loading, when ω ≤ ḣ m /h 0 , from Equation (7) the contact
area will be non-decreasing during the whole process; as the frequency exceeds the
critical value ωc = ḣ m /h 0 , the contact area increases and decreases with time as
a result of the applied harmonic load so that the Ting approach should be adopted.
Nevertheless, as will be discussed in Section 4.2, under certain condition when
ω > ωc , the solutions derived from the methods by Lee and Radok, and Ting are
very close, even though the condition for the Lee–Radok approach is not justified.
Since a closed-form solution derived from the Lee–Radok approach exists, while
only numerical solution can be obtained using the Ting approach, the formulas
derived for a harmonic load superimposed on a step loading from the Lee–Radok
approach could be used to estimate the complex viscoelastic functions in the regime
of linear viscoelasticity when ω > ωc .
3. Measurements
We conducted two independent measurements, namely nanoindentation and DMA
measurements, to find the complex viscoelastic functions of the same materials. The
results would be compared to examine the measurement technique by nanoindenta-
tion. We describe in this section first DMA measurements and then nanoindentation
measurements.
Figure 5. Nanoindentation output from oscillation on a ramp loading at 10 Hz. (a) Carrier load-
depth curve. (b) History of amplitude of harmonic load. (c) Response of harmonic displacement
amplitude. (d) Out-of-phase angle with correction.
where P(t) is given by Equation (4). For the second half of the cycle, in which the
contact area is decreasing, the solution derived from the Ting approach is
a(t) = a(t1 ), (16)
where t1 < T /2 and t > T /2, with T being the period of one cycle, and t1 , as
shown in Figure 6, is obtained by
t1
8
P(t) = G(t − τ ) d(a 3 (τ )), (17)
3(1 − ν)R 0
where G(t) is the shear relaxation modulus. Using the creep functions of PC and
PMMA, the contact radius in one cycle can be determined after the sinusoidal
response has reached the steady state. Figure 6 shows the results of contact radius
Figure 6. Comparison of contact radius results at 75 Hz using the Lee–Radok approach and the
Ting approach. (a) Contact radius computed for PMMA under a harmonic load superimposed
on a step loading. (b) Contact radius computed for PC under a harmonic load superimposed
on a step loading.
358 G. HUANG ET AL.
Figure 7. Nanoindentation output from oscillation on a step loading at 75 Hz. (a) Carrier
load-depth curve. (b) The history of harmonic load amplitude. (c) Response of harmonic
displacement amplitude. (d) Out-of-phase angle with correction.
MEASUREMENTS OF VISCOELASTIC FUNCTIONS OF POLYMERS 359
carrier load-depth curve (without the harmonic component) for both PC and PMMA.
Figures 7(b) and 7(c) illustrate the amplitudes of harmonic load and displacement
at 75 Hz for PC and PMMA. As shown in Figure 7(b), at initial stage of the contact
between the indenter and sample surface, the amplitude of harmonic load generally
increased until t = 50 s, and was maintained constant thereafter. Accordingly, as
shown in Figure 7(c), there was an increase in amplitude of harmonic displacement
before t = 50 s. So the condition of Equation (4) is precise after about 50 s. A
steady state was reached at about t = 150 s. As shown in Figure 7, after t = 150
s, h 0 / h m is much less than 1%, thus the condition of using Equation (10) is
satisfied. Figure 7(d) shows the out-of-phase angle between harmonic load and
harmonic displacement for both PC and PMMA.
In order to ensure that the deformation of the polymer samples is in the lin-
early viscoelastic regime, the indentation depth into the sample surface for PC and
PMMA materials was controlled to within the limit of linearity. It was found that
in indentation within the limit of linearity, the deformation of PC and PMMA is
linearly viscoelastic, indicated by the fact that complex compliance is independent
of the magnitude of indentation carrier load. It is also noted that the condition
h/R < 0.16 needs to be satisfied as well to justify the use of Equations (2) and (3).
This leads to h < 544 nm when a spherical indenter with a tip radius of 3.4 µm was
used in our tests. As a result, in this investigation experimental data associated with
h(t) < 544 nm were processed to determine the complex viscoelastic functions in
the frequency domain.
For time-dependent materials under dynamic loading, the out-of-phase angle
between the harmonic load and displacement plays an important role in the com-
putation of complex compliance. Consequently, we discuss next the correction on
out-of-phase angle.
As presented by Pethica and Oliver (1989), the out-of-phase angle between the
force and displacement φ can be computed by
ω(Ci + Cc )
tan φ = , (18)
k − mω2
ωCi
tan δ = . (19)
k − mω2
Figure 8. Complex compliance in shear from nanoindentation under a harmonic load super-
imposed on a ramp loading.
MEASUREMENTS OF VISCOELASTIC FUNCTIONS OF POLYMERS 361
Figure 9. Complex compliance in shear from nanoindentation under a harmonic load super-
imposed on a step loading.
storage and loss compliance for PC are only 0.8 and 3.8%, respectively, while for
PMMA, the storage compliance decreases by 16.6%, and loss compliance decreases
by 35.4% between 3 Hz and 260 Hz.
In the computation of both storage compliance and loss compliance, Poisson’s
ratio is usually time-dependent for polymeric materials. In these measurements,
the time-scale is not large; the maximum time duration is less than 300 seconds.
Within this time period, it has been demonstrated that the Poisson’s ratio does not
change significantly (Lu et al., 1997). Therefore, we assumed a constant Poisson’s
ratio (ν = 0.3) in the computation of the complex compliance for both PC and
PMMA materials. The natural frequency of the indentation system is 180 Hz, thus
we avoided testing near the resonant frequency.
We turn next to the comparison in results from Loubet et al. and the new method
proposed in this paper. The method by Loubet et al. (1995) on measurement of the
uniaxial complex modulus is based on the analogy between dynamic indentation
and uniaxial dynamic analysis. Nanoindentation on viscoelastic materials, however,
is a complex viscoelastic problem involving moving contact interface, and needs to
be analyzed to derive formulas for the computation of viscoelastic functions in the
frequency-domain. Based on their method, the uniaxial storage modulus at 75 Hz
for both PC and PMMA were computed using the data of stiffness, damping and
contact area from nanoindentation measurements, as shown in Figure 1. During the
tests, a ramp loading superimposed by a small harmonic loading was applied. Since
tests were under a single frequency (75 Hz) within a time frame less than 150 s,
the storage modulus computed using Equation (1) was plotted with respect to time.
Also plotted in Figure 1 are the DMA data of uniaxial storage modulus at 75 Hz,
362 G. HUANG ET AL.
5. Conclusions
A method to measure the complex compliance has been presented using nanoin-
dentation with a spherical indenter for linearly viscoelastic materials. Following
the Hertzian solutions for indentation in linear elasticity and the consideration of
the Lee–Radok approach for a moving boundary problem in linear viscoelasticity,
formulas for the components of the complex compliance function in the frequency-
domain have been derived based on the load-displacement relation for linearly
viscoelastic materials under harmonic loading. The formulas should be used under
the conditions that h/R < 0.16 and the nanoindentation depth is below the limit of
linearity. When a ramp loading is used as the carrier load, the formulas are exact
under the condition that the loading rate is high enough so that the condition of non-
decreasing contact area is satisfied. While a step loading is used as the carrier load,
the formulas are considered to be approximate when the frequency is higher than
a critical value such that decreasing contact area occurs; in the case of decreasing
contact area, the Ting approach was used to find the solution in the steady state at
75 Hz. The difference between solutions obtained using the Lee–Radok approach
and the Ting approach was estimated. Very close correlation was found for the test
conditions used in this study. Dynamic nanoindentation tests on PC and PMMA
materials were performed, using a spherical indenter, to determine the amplitudes
of oscillating load and displacement, as well as the out-of-phase angle of displace-
ment with respect to the harmonic force output. The complex compliance functions
in the frequency-domain were determined using the proposed method and results
MEASUREMENTS OF VISCOELASTIC FUNCTIONS OF POLYMERS 363
were compared with conventional data obtained from DMA tests for the same PC
and PMMA materials. The condition for non-decreasing contact area, v0 ≥ P0 ω,
is satisfied under a ramp carrier load so that the complex compliance formulas de-
rived from the Lee–Radok approach are appropriate. Under a step carrier load, the
condition of non-decreasing contact is satisfied up to the frequency limit, ḣ m /h 0 ;
at frequencies higher than the frequency limit, complex compliance data could be
considered to be an approximation. In both ramp and step carrier loading conditions,
a good agreement between nanoindentation results and conventional data has been
reached, indicating the validity of the proposed method for measuring the complex
compliance function in the frequency-domain using dynamic nanoindentation.
Acknowledgments
This work was initially supported by the National Science Foundation under CMS-
9872350 and CMS-9985060. Funds from NASA under grant NNL04AA4ZG with
Dr. Thomas S. Gates as the technical manager allowed the completion of this study.
The authors are grateful to Dr. Brian P. Grady at the University of Oklahoma for
conducting DMA tests.
References
Cheng, L., Xia, X., Yu, W., Scriven, L.E. and Gerberich, W.W., ‘Flat-punch indentation of viscoelastic
material’, J. Polym. Sci [B]: Polym. Phys. 38, 2000, 10–22.
Ferry, J.D., ‘Mechanical properties of substances of high molecular weight. VI. Dispersion in con-
centrated polymer solutions and its dependence on temperature and concentration’, J. Am. Chem.
Soc. 72, 1950, 3746–3752.
Giannakopoulos, A.E. ‘Strength analysis of spherical indentation of piezoelectric materials’, J. Appl.
Mech. 67, 2000, 409–416.
Hertz, H., ‘Über die beruhrung fester elastischer körper’, J. für die Reine und Angewandte Mathematik
92, 1881, 156–171.
Knauss, W.G. and Zhu, W., ‘Nonlinearly viscoelastic behavior of polycarbonate. I. Response under
pure shear’, Mech. Time-Depend. Mater. 6, 2002, 231–269.
Lee, E.H. and Radok, J.R.M., ‘The contact problem for viscoelastic bodies’, J. Appl. Mech. 27, 1960,
438–444.
Li, X. and Bhushan, B., ‘A review of nanoindentation continuous stiffness measurement technique
and its applications’, Mater. Char. 48, 2002, 11–36.
Ling, F.F., Lai, W.M. and Lucca, D.A., Fundamental of Surface Mechanics, Springer-Verlag, New
York, 2002.
Loubet, J.L., Lucas, B.N. and Oliver, W.C., ‘Some measurements of viscoelastic properties with the
help of nanoindentation’, in International Workshop on Instrumental Indentation, San Diego, CA,
April 1995, D.T. Smith (ed.), 1995, pp. 31–34.
Lu, H., Zhang, X. and Knauss, W.G., ‘Uniaxial, shear and Poisson relaxation and their conversion to
bulk relaxation: studies on poly (methyl methacrylate)’, Polymer Engineering Science 37, 1997,
1053–1064.
Lu, H., Wang, B., Ma, J., Huang, G. and Viswanathan, H., ‘Measurement of creep compliance of solid
polymers by nanoindentation’, Mech. Time-Depend. Mater. 7, 2003, 189–207.
364 G. HUANG ET AL.
Lucas, B.N., Oliver, W.C., and Swindeman, J.E., ‘The dynamics of frequency-specific, depth-sensing
indentation testing’, Mater. Res. Soc. Symp. Proc. 522, San Francisco, Moody, N.R. et al. (ed.),
1998, 3–14.
Oliver, W.C. and Pharr, G.M., ‘An improved technique for determining hardness and elastic modulus
using load and displacement sensing indentation experiments’, J. Mater. Res. 7, 1992, 1564–1583.
Oyen-Tiesma, M., Toivola, Y.A. and Cook, R.F., ‘Load-displacement behavior during sharp indenta-
tion of viscous-elastic-plastic materials’, in Fundamentals of Nanoindentation and Nanotribology
II, Baker, S.P., Corcoran, S., Moody, G.N.R. and Cook, R.F. (ed.), Warrendale, PA, 2001, Q.5.1.
Pethica, J.B. and Oliver, W.C., ‘Tip surface interactions in STM and AFM’, Phys. Scr. T19, 1987,
61–66.
Pethica, J.B. and Oliver, W.C., ‘Mechanical properties of nanometer volumes of material: Use of the
elastic response of small area indentations’, in Materials Research Society Symposium Proceed-
ings 130, Pittsburgh, Bravman, J.C. et al. (ed.), 1989, 13–23.
Sane, S.B. and Knauss, W.G., ‘The time-dependent bulk response of poly(methyl methacrylate)’,
Mech. Time-Depend. Mater. 5, 2001, 293–324.
Sneddon, I.N., ‘The relation between load and penetration in the axisymmetric boussinesq problem
for a punch of arbitrary punch’, Inter. J. Eng. Sci. 3, 1965, 47–56.
Syed, S.A., Wahl, K.J. and Colton, R.J., ‘Nanoindentation and contact stiffness measurement using
force modulation with a capacitive load-displacement transducer’, Rev. Sci. Instrum. 70(5), 1999,
2408–2413.
Ting, T.C.T., ‘The contact stresses between a rigid indenter and a viscoelastic half-space’, J. Appl.
Mech. 33, 1966, 845–854.