STRESSES IN THIN FILMS AND THEIR RELAXATION, KRAFT, GAO 387
After the stress reaches s
y
on loading, plastic deformation occurs with a
linear increase in flow stress. Due to this increase, a back stress s
b
builds up. On reverse loading, the yield strength in compression is
lowered by this back stress and the film yields at s
reverse
s
y
s
b
.
Figure 8.16 also shows that this predicted behavior is in agreement with
the cyclic stress-strain behavior of the Cu films on polyimide substrates.
Shen et al. have shown in their original work that their model can be
applied fairly well to model thermal stress cycles of Cu films when a
temperature-dependent yield strength is assumed. The incremental
increase in stress, that is, the hardening rate, in the elastic-plastic regime
is given by:
∂s
,
(16)
where M is the biaxial modulus of the film, e
*
is the hardening parame-
ter, and s
y
(T) is the temperature-dependent yield strength. The yield
strength was assumed to decrease linearly with temperature, that is,
s
y
(T) s
o
(1 TT
*
), where s
o
is the yield strength at T 0 K and T
*
is the temperature at which the yield strength becomes zero. Hence, e
*
,
s
o
, and T
*
can be used to fit the model empirically to experimental data.
Nix and Leung
[50]
have adapted this approach and tried to reduce the
empirical fit parameters. They suggested that the yield strength is
related to the dislocation motion in the film, as given by Eq. (11).
Therefore, its temperature dependence is related to the temperature
dependence of the shear modulus, which can be assumed to be G
G
o
(1 c
*
(T 300)T
m
), where G
o
is the shear modulus at room tem-
perature, c
*
is a constant close to 0.5 for most FCC metals, and T
m
is the
melting temperature.
Figure 8.17 shows a quantitative comparison of this approach to
the data from Fig. 8.15. Again, a reasonable agreement between
experiment and model can be seen. In particular, the plateau in com-
pression is well reproduced. Note that the hardening rate, which is
adjusted to fit the experimental data, is of the order of the shear
modulus of Cu, which is about two orders of magnitude larger than
the hardening rate in bulk Cu. Furthermore, this approach does not
allow us to account for any time-dependent deformation; therefore, it
is not possible to account for strain rate effects or stress relaxation.
More recently, Kraft et al.
[51]
suggested that both thermally activated
dislocation glide and strain hardening be included by modifying
a∂T
M
1
s
e
y
(
*
T)