Models for thin films ====================== .. include:: includes.rst The following parts give a short overview of models existing in the literature used for the extraction of mechanical properties of thin films deposited on a substrate from indentation experiments with conical indenters. Before everything, it is is worth to mention the work performed by Jennett N.M. and Bushby A.J. [#Jennett_2001]_, about nanoindentation test on coatings, during the European project INDICOAT (SMT4-CT98-2249). Progress from this project help in the development of the ISO standard (ISO 14577 - 1-4). The ISO 14577 - 4 is dedicated to nanoindentation on coatings. * `ISO 14577 - 1 , "Metallic materials -- Instrumented indentation test for hardness and materials parameters -- Part 1: Test method", (2002). `_ * `ISO 14577 - 2 , "Metallic materials -- Instrumented indentation test for hardness and materials parameters -- Part 2: Verification and calibration of testing machines", (2002). `_ * `ISO 14577 - 3 , "Metallic materials -- Instrumented indentation test for hardness and materials parameters -- Part 3: Calibration of reference blocks", (2002). `_ * `ISO 14577 - 4 , "Metallic materials -- Instrumented indentation test for hardness and materials parameters -- Part 4: Test method for metallic and non-metallic coatings", (2007). `_ Some authors overviewed/reviewed already the nanoindentation technique applied to coatings : - `Pharr G.M. and Oliver W.C., "Measurement of Thin Film Mechanical Properties Using Nanoindentation" (1992). `_ - `Menčík J. et al., "Determination of elastic modulus of thin layers using nanoindentation" (1997). `_ - `Nix W.D., "Elastic and plastic properties of thin films on substrates: nanoindentation techniques" (1997). `_ - `Van Vliet K.J. and Gouldstone A., "Mechanical Properties of Thin Films Quantified Via Instrumented Indentation" (2001). `_ - `Šimůnková Š., "Mechanical properties of thin film–substrate systems" (2003). `_ - `Bull S.J., "Nanoindentation of coatings" (2005). `_ - `Fischer-Cripps, A.C., "Nanoindentation 3rd Ed." (2011) `_ - `Wen, W. et al., "Determination of material properties of thin films and coatings using indentation tests: a review" (2017) `_ Nanoindentation tests on thin films ###################################### Composite reduced Young's modulus and composite hardness ---------------------------------------------------------- For indentation test on a coated specimen or on a multilayer sample (e.g.: thin films deposited on a substrate), the evolution of the Young's modulus or the hardness calculated with models used for bulk materials, is function of the material properties and the thickness :math:`t` (in :math:`\text{m}`) of each underlying film (substrate included), and the properties and the geometry of the indenter. .. figure:: ./_pictures/multilayerSEM.png :scale: 60 % :align: center *SEM cross-sectional observation of a multilayer sample.* .. figure:: ./_pictures/evolution_EandH_thinfilm.png :scale: 60 % :align: center *Typical evolution of Young's modulus and hardness for a coated specimen as a function of the normalized indentation depth.* Thus, the composite reduced Young's modulus :math:`E^{'}` and the composite hardness :math:`H` calculated with the models used for bulk materials, can generally be expressed as a combination of respectively the reduced Young's moduli :math:`(E^{'}_\text{f})` or the hardness :math:`(H_\text{f})` of each underlayer and respectively the reduced Young's modulus :math:`(E^{'}_\text{s})` or the hardness :math:`(H_\text{s})` of the substrate. The reduced Young's moduli and the hardness are in :math:`\text{GPa}`. .. math:: E^{'} = f(E^{'}_{\text{f},{i \rightarrow N}}, t_{i \rightarrow N}, E^{'}_\text{s}) :label: youngs_modulus_evolution .. math:: H = f(H_{\text{f},{i \rightarrow N}}, t_{i \rightarrow N}, H_\text{s}) :label: hardness_evolution With :math:`i` the indice of the layer and :math:`N` the total number of layers. Finally, it is important to avoid confusion between `indentation size effect `_. and a decrease of hardness as a function of indentation depth observed in the case of a hard thin film deposited on a soft substrate. It might happen that both effects occur during indentation experiment, which complicates the analysis [#Zhang2007]_. Indentation contact topography --------------------------------- For nanoindentation tests on thin films, the contact topography is function of both thin film and substrate properties. .. figure:: ./_pictures/contact_topography_thinfilm.png :name: contact_topography_thinfilm :scale: 60 % :align: center *Schematic depiction of a) "pile-up" and b) "sink-in" observed during thin film indentation.* The :numref:`contact_topography_thinfilm` a ("pile-up") is typical of the case of a soft film on a hard substrate and the :numref:`contact_topography_thinfilm` b ("sink-in") of a hard film on a soft substrate [#ChenVlassak2001]_. The pile-up can be emphasized in case ot thin films, because of the material confinement by the substrate. To determine the depth of contact, `the same models `_ described for bulk material indentation are used, especially when the ratio of the indentation depth over the coating thickness remains low (<10%, see the rule proposed by `Bückle `_). But, some quantitative approaches are proposed in the literature to take into account both pile up and substrate effects in the material properties estimation of the coating [#Zhou2008]_. Corrections to apply for thin film indentation ------------------------------------------------ During nanoindentation tests of thin film on substrate, the thickness of the film beneath the indenter is smaller than its original value, because of plastic flow during loading. The use of the original film thickness :math:`t` in the regression model cause a systematic shift or distortion of the Young's modulus curve. A correction proposed by Menčík et al. can be applied, assuming a rigid substrate and determining the effective thickness :math:`t_\text{eff}` (in :math:`\text{m}`) [#Menčík_1997]_, [#Saha_2002]_, [#Chen_2004]_, [#Bec_2006]_ and [#Li_2010]_. .. math:: \pi a^2 t_\text{eff} = \pi a^2 t - V :label: general_thickness_correction With :math:`V` the volume displaced by the indenter and approximated by :math:`\pi a^2 h_\text{c} / 3`, for a conical indenter and contact depths :math:`h_\text{c}` smaller than the film thickness. .. math:: t_\text{eff} = t - \frac{h_\text{c}}{3} :label: thickness_correction .. figure:: ./_pictures/thickness_correction.png :scale: 60 % :align: center *Indentation penetration of a thin film on a sample.* Recently, Li et al. proposed to express the local thinning effect as [#Li_2010]_: .. math:: t_\text{eff} = t - \eta h_\text{c} :label: thickness_correction_Li With :math:`\eta` a parameter depending on the mechanical properties of the film and the substrate, and on the geometry of the indenter. Preliminary finite element calculations show that :math:`\eta` should be independent of indentation depth and that its value ranges from :math:`0.3` to :math:`0.7` for materials that do not work harden. Method of the tangent ----------------------- First of all, the easiest way to assess approximately the mechanical properties of the film is to use the tangential model. If one plots a tangent from the diagram's zero point to the curves, then the tangential point is the Young's modulus or the hardness of the film (see :numref:`tangente-method`). .. figure:: ./_pictures/tangente-method.png :name: tangente-method :scale: 60 % :align: center *Schematic depiction of the tangent method to assess film's properties.* Elastic properties of a thin film on a substrate ################################################## Several models to analyse indentation tests on bilayer sample and multilayer sample and to extract intrinsic material properties of the upper film (or the top coating) are detailed in the following part. Some of the following models developed are detailed in the chapter 8 ("Nanoindentation of Thin Films") of the book "Nanoindentation" written by A.C. Fischer-Cripps [#FischerCripps_2004]_. Bückle (1961) -------------- Bückle proposed an empirical law ("rule of thumb") to characterize thick coatings on a substrate [#Bückle_1961]_. It is possible to estimate the Young's modulus of the coating for indentation depth lower than ":math:`10\text{%}`" of the film thickness. But, because of the imperfections of the indenters, the roughness and the surface pollution, it is more meaningful to use this rule of the ":math:`10\text{%}`", for film thicker than :math:`500\text{nm}`. Doerner and Nix (1986) ------------------------- The model of Doerner and Nix is detailed in many papers [#DoernerNix_1986]_, [#King_1987]_, [#Pharr_1992]_ and [#Saha_2002]_, and is described by the following equation : .. math:: \frac{1}{E^{'}} = \frac{1}{E_\text{f}^{'}} + \left(\frac{1}{E_\text{s}^{'}} - \frac{1}{E_\text{f}^{'}}\right) e^{-\alpha\left(x\right)} :label: doerner_nix With :math:`x=t/h_\text{c}` and :math:`\alpha` an empirically constant determined using the method of least squares. The equation was modified by King [#King_1987]_ with the replacement of :math:`t/h_\text{c}` by :math:`t/a_\text{c}`, and then by Saha and Nix [#Saha_2002]_ with the replacement of :math:`t/h_\text{c}` by :math:`(t-h_\text{c})/a_\text{c}`. An empirical formulae based on the model of Doerner and Nix was proposed by Chen et al. in 2004 [#Chen_2004]_: .. math:: \frac{1}{E^{'}} = \frac{1}{E_\text{f}^{'}}\left[1-e^{-\alpha\left(x\right)^{(2/3)}}\right] + \frac{1}{E_\text{s}^{'}}e^{-\alpha\left(x\right)^{(2/3)}} :label: chen With :math:`x=t/h_\text{c}` and :math:`\alpha` an empirically constant determined using the method of least squares. Find here the |matlab| function for the Doerner and Nix model [#DoernerNix_1986]_: `model_doerner_nix.m `_. Find here the |matlab| function for the Doerner and Nix model modified by King [#King_1987]_: `model_doerner_nix_king.m `_. Find here the |matlab| function for the Doerner and Nix model modified by Saha [#Saha_2002]_: `model_doerner_nix_saha.m `_. Find here the |matlab| function function for the Doerner and Nix model modified by Chen [#Chen_2004]_: `model_chen.m `_. Gao et al. (1992) ------------------ The model of Gao is described by the following equation [#Gao_1992]_ : .. math:: E^{'} = E^{'}_\text{s} + \left(E^{'}_\text{f} - E^{'}_\text{s}\right) \phi_{Gao_0}\left(x\right) :label: gao .. math:: \phi_{Gao_0} = \frac{2}{\pi} arctan \frac{1}{x} + \frac{1}{{2\pi\left(1-\nu_\text{c}\right)}} \left[\left(1-2\nu_\text{c}\right)\frac{1}{x}ln \left(1 + x^2 \right) - \frac{x}{{1+x^2}}\right] :label: phi_gao_0 .. math:: \nu_{c} = 1 + \left[{\frac{\left(1-\nu_\text{s}\right) \left(1-\nu_\text{f}\right)} {1-\left(1-\phi_{Gao_1}\right)\nu_\text{f} - \phi_{Gao_1}\nu_\text{s}}}\right] :label: poisson_ratio_composite With :math:`\nu_{c}` the composite Poisson's ratio, :math:`\nu_{s}` the Poisson's ratio of the substrate and :math:`\nu_{f}` the Poisson's ratio of the thin film. .. math:: \phi_{Gao_1} = \frac{2}{\pi} arctan \frac{1}{x} + \frac{1}{{x\pi}}ln \left(1+x^2 \right) :label: phi_gao_1 With :math:`x=a_\text{c}/t`. Find here the |matlab| function for the weighting function :math:`\phi_{Gao_0}` : `phi_gao_0.m `_. Find here the |matlab| function for the weighting function :math:`\phi_{Gao_1}` : `phi_gao_1.m `_. Find here the |matlab| function for :math:`\nu_{c}` the composite Poisson's ratio : `composite_poissons_ratio.m `_. Find here the |matlab| function for the Gao et al. model : `model_gao.m `_. Menčík et al. (1997) ----------------------- Menčík et al. proposed the following structures to express the combination of :math:`E^{'}_\text{f}` and :math:`E^{'}_\text{s}` [#Menčík_1997]_. .. math:: E^{'} = E^{'}_\text{s} + \left(E^{'}_\text{f} - E^{'}_\text{s}\right) \phi\left(x\right) :label: youngs_modulus_evolution_phi .. math:: E^{'} = E^{'}_\text{s} + \left(E^{'}_\text{f} - E^{'}_\text{s}\right) \psi\left(x\right) :label: youngs_modulus_evolution_psi Where :math:`x` is the ratio of the contact radius (:math:`a_\text{c}`) or the contact depth (:math:`h_\text{c}`), to the film thickness (:math:`t`), and :math:`\phi` and :math:`\psi` are weight functions of the relative penetration :math:`x`. :math:`\phi` is equal to :math:`1` when :math:`x` is equal to :math:`0` and :math:`0` when :math:`x` is infinite. .. note:: If the difference between Poisson's ratio of the thin film and substrate is small, the values for uniaxial loading Young's moduli, :math:`E`, :math:`E_\text{f}`, :math:`E_\text{s}` can be used in previous equation. Menčík et al. (linear model) (1997) ------------------------------------ Menčík described too the linear model by the following expression [#Menčík_1997]_ : .. math:: E^{'} = E^{'}_\text{f} + \left(E^{'}_\text{s} - E^{'}_\text{f}\right)\left(x\right) :label: mencik_linear With :math:`x=a_\text{c}/t`. Find here the |matlab| function for the Menčík et al. linear function : `model_menick_linear.m `_. Menčík et al. (exponential model) (1997) ------------------------------------------- Menčík described the exponential model by the following expression [#Menčík_1997]_ : .. math:: E^{'} = E^{'}_\text{s} + \left(E^{'}_\text{f} - E^{'}_\text{s}\right) e^{-\alpha\left(x\right)} :label: mencik_exponential With :math:`x=a_\text{c}/t` and :math:`\alpha` is an empirically constant determined using the method of least squares. Find here the |matlab| function for the Menčík et al. exponential function : `model_menick_exponential.m `_. Menčík et al. (reciprocal exponential model) (1997) ----------------------------------------------------- Menčík described the reciprocal exponential model by the following expression [#Menčík_1997]_ : .. math:: \frac{1}{E^{'}} = \frac{1}{E_\text{s}^{'}} + \left(\frac{1}{E_\text{f}^{'}} - \frac{1}{E_\text{s}^{'}}\right) e^{-\alpha\left(x\right)} :label: mencik_reciprocal_exponential With :math:`x=a_\text{c}/t` and :math:`\alpha` is an empirically constant determined using the method of least squares. Find here the |matlab| function for the Menčík et al. reciprocal exponential function : `model_menick_reciprocal_exponential.m `_. Perriot et al. (2003) ----------------------- The following system of equation describes the model developed by Perriot et al. [#Perriot_2004]_ : .. math:: E^{'} = E^{'}_\text{f} + \frac{E^{'}_\text{s}-E^{'}_\text{f}} {1+\left(\frac{k_0}{x}\right)^n} :label: perriot_barthel .. math:: log(k_0) = -0.093+0.792log\left(\frac{E^{'}_\text{s}} {E^{'}_\text{f}}\right) + 0.05\left[log\frac{E^{'}_\text{s}}{E^{'}_\text{f}}\right]^2 :label: perriot_barthel_k0 With :math:`x=a_\text{c}/t`, and :math:`k_0` and :math:`n` are adjustable constants determined using the method of least squares. Find here the |matlab| function for the Perriot et al. model : `model_perriot_barthel.m `_. Jung et al. (2004) --------------------- Jung et al. [#Jung_2004]_ have adapted for conical indentation of thin films, the simple empirical approach of Hu and Lawn [#Hu_1998]_ developed initially for spherical indentation on bilayer structures. The following power-law relationship allows the evaluation of the Young's modulus of a thin film deposited on a substrate from nanoindentation experiments : .. math:: E = E_\text{s} {\left(\frac{E_\text{f}}{E_\text{s}}\right)}^L :label: jung with :math:`L` is the exponent term described by a dimensionless function : .. math:: L = \frac{1}{\left[1+A{x}^B\right]} :label: jung_sigmoidal With :math:`x=h_\text{c}/t` and where :math:`A` and :math:`B` are adjustable coefficients. Jung et al. founded :math:`A=3.76` and :math:`B=1.38` after regression fits of :eq:`jung` to different data sets. These coefficients are not universal and need to be "calibrated" with experimental data or with finite element data for specified material systems. Finally, to be more consistent with other analytical models implemented in this toolbox, the model of Jung is modified by using the reduced form of the Young's moduli : .. math:: E^{'} = E^{'}_\text{s} {\left(\frac{E^{'}_\text{f}}{E^{'}_\text{s}}\right)}^L :label: jung_reduced Find here the |matlab| function for the sigmoidal function used in the Jung's model : `sigmoidal_jung.m `_. Find here the |matlab| function for the Jung et al. model : `model_jung.m `_. Bec et al. (2006) ------------------- The elastic model of Bec et al. is based on indentation by a rigid cylindrical punch (radius :math:`a_\text{c}`) of a homogeneous film deposited on a semi-infinite half space [#Bec_2006]_. This system is modelled by two springs connected in series : .. figure:: ./_pictures/Bec_bilayer_model.png :scale: 60 % :align: center *Schematic description of the bilayer model of Bec et al.* .. math:: K_\text{f} = \pi a_\text{c}^2 \frac{E^{'}_\text{f}}{t} :label: bec_kf .. math:: K_\text{s} = 2 a_\text{c} E^{'}_\text{s} :label: bec_ks .. math:: K_\text{z} = 2 a_\text{c} E^{'} :label: bec_kz .. math:: \frac{1}{K_\text{z}} = \frac{1}{{f_\text{f}\left(a_\text{c}\right) K_\text{f}}} + \frac{1}{{f_\text{s}\left(a_\text{c}\right)K_\text{s}}} :label: bec_sum_kf_ks .. math:: f_\text{f}\left(a_\text{c}\right) = f_\text{s}\left(a_\text{c}\right) = 1 + {\frac{2t}{\pi a_\text{c}}} :label: bec_f_ac .. math:: \frac{1}{2 a_\text{c} E^{'}} = \frac{t}{\left(\pi a_\text{c}^2 + 2ta_\text{c}\right)E^{'}_\text{f}} + \frac{1}{2\left(a_\text{c} + \frac{2t}{\pi}\right)E^{'}_\text{s}} :label: bec_Eeq Find here the |matlab| function for the Bec et al. model : `model_bec.m `_. Korsunsky and Constantinescu (2009) ---------------------------------------------- Korsunsky and Constantinescu proposed a simple model response function to the analysis of indentation of elastic coated systems [#Korsunsky_1998]_ and [#Korsunsky_2009]_. They expressed the reduced Young's modulus in terms of a linear law of mixtures of the form: .. math:: E^{'} = E_\text{1}^{'} + \frac{E_\text{2}^{'} - E_\text{1}^{'}}{1 + \left(\frac{h_\text{t}}{\beta_{0}t}\right)^{\eta}} :label: korsunsky_model Here :math:`E_\text{1}^{'}`, :math:`E_\text{2}^{'}`, :math:`\eta` and :math:`\beta_{0}` are positive constants to be determined from fitting. It may be expected that for very shallow indentation (:math:`\beta_{0} \ll 1`), the corresponding parameter :math:`E_\text{1}^{'}` ought to approach the Young's modulus of the coating :math:`E_\text{f}^{'}`. Similarly, one might also expect that for very deep indentation (:math:`\beta_{0} \gg 1`), the corresponding parameter :math:`E_\text{2}^{'}` ought to approach the Young's modulus of the substrate :math:`E_\text{s}^{'}`. Values of :math:`\eta` and :math:`\beta_{0}` depend on the indentation boundary condition, if the coating is a defined as a freely sliding layer or as a perfectly bonded layer. Hay et al. (2011) ------------------- The present model of Hay et al. [#Hay_2011]_ is a development of the Song–Pharr model [#Rar_2002]_ and [#Xu_2006]_, which is already inspired by the Gao model [#Gao_1992]_. .. math:: \frac{1}{\mu_\text{c}} = \left(1-\phi_{Gao_0}\right) \frac{1}{\mu_\text{s} + F\phi_{Gao_0}\mu_\text{f}} + \phi_{Gao_0}\frac{1}{\mu_\text{f}} :label: hay_model Where :math:`\mu_\text{c}` (in :math:`\text{GPa} = \text{N/m}^2`) is the composite shear modulus calculated from the composite Young's modulus as : .. math:: \mu_\text{c} = \frac{E}{2\left(1+\nu_\text{c}\right)} :label: hay_muc Where :math:`\nu_\text{c}` is the composite Poisson's ratio given in Gao's model. .. math:: E = \left(1 - \nu_\text{c}^2\right)E^{'} :label: hay_Ec Knowing :math:`\mu_\text{c}`, it is possible to calculate :math:`\mu_\text{f}` : .. math:: \mu_\text{f} = \frac{-B + \sqrt{B^2-4AC}}{2A} :label: hay_muf Where :math:`A = F \phi_{Gao_0}` .. math:: B = \mu_\text{s} - \left(F \phi_{Gao_0}^2 - \phi_{Gao_0} + 1\right)\mu_\text{c} :label: hay_B With :math:`F = 0.0626`, a constant obtained from finite element simulations. .. math:: C = -\phi_{Gao_0}\mu_\text{c}\mu_\text{s} :label: hay_C .. math:: \mu_\text{s} = \frac{E_\text{s}}{2\left(1+\nu_\text{s}\right)} :label: hay_mus Finally, the Young's modulus of the film is calculated from the shear modulus and Poisson’s ratio of the film : .. math:: E_\text{f} = 2\mu_\text{f}\left(1 + \nu_\text{f}\right) :label: hay_Ef Find here the |matlab| function for the Hay et al. model : `model_hay.m `_. Bull (2014) ------------------- This model was originally proposed by Bull S.J. in 2011, in a paper about the mechanical characterization of ALD Alumina coatings [#Bull_2011]_. Then, this simple method to determine the elastic modulus of a coating on a substrate using nanoindentation was developed in another paper in 2014 [#Bull_2014]_. This model is based on the load support of a truncated cone of material beneath the indenter. .. math:: E^{'} = \frac{F_\text{c}}{2a_\text{c}\left(h_\text{c} + h_\text{s} \right)} :label: bull_contact_modulus .. math:: h_\text{c} = \frac{F_\text{c}}{\pi E_\text{f}} \left[ \frac{1}{a_\text{c} tan\alpha} - \frac{1}{a_\text{c} tan\alpha + t_\text{f} tan^2\alpha} \right] :label: bull_disp_coating .. math:: h_\text{s} = \frac{F_\text{c}}{\pi E_\text{s}} \left[ \frac{1}{a_\text{c} tan\alpha + t_\text{f} tan^2\alpha} - \frac{1}{a_\text{c} tan\alpha + (t_\text{f} + t_\text{s}) tan^2\alpha} \right] :label: bull_disp_substrate Where :math:`E_\text{f}` and :math:`E_\text{s}` are the Young's Modulus of the coating and substrate, :math:`t_\text{f}` and :math:`t_\text{s}` are the coating and substrate thickness, and :math:`\alpha` is the semi-angle of the cone material which supports the load. In fact, by assuming that the material thickness is very much greater than the contact radius, it is possible to replace in the previous equation :math:`tan\alpha` by :math:`2\pi = 32.48°`. Finally, by assuming that the substrate is very much thicker than the coating :math:`(t_\text{s} >> t_\text{f})`, the equation :eq:`bull_contact_modulus` can be rewritten : .. math:: E^{'} = \frac{1}{\frac{1}{E_\text{f}} \left[\frac{2t_\text{f}}{\pi a_\text{c} + 2t_\text{f}} \right] + \frac{1}{E_\text{s}} \left[\frac{\pi a_\text{c}}{\pi a_\text{c} + 2t_\text{f}} \right]} :label: bec_bull_contact_modulus This last equation from Bull S.J. is exactly the same as equation :eq:`bec_Eeq`, proposed by Bec et al. in 2006, with a unique difference which is the use by Bull S.J. of the non reduced form of the Young's moduli of the coating and the substrate. Argatov I.I. and Sabina F.J. (2016) ------------------------------------- Argatov I.I. and Sabina F.J. proposed an analytical model, based on the first-order asymptotic solution to the unilateral frictionless indentation problem. This model takes into account the effect of a compliant substrate in the small-scale indentation of thin coatings [#Argatov_2016]_. Elastic properties of a thin film on a multilayer system ############################################################### In 2008, Pailler-Mattei et al. proposed an extension of the Bec's model to a bilayer system deposited on a substrate [#PaillerMattei_2008]_. But more recently, Mercier et al. established a generalization of the Bec's model to :math:`N+1` layers sample. Mercier et al. (2010) ---------------------- The elastic model of Mercier et al. for a multilayer sample on :math:`N+1` layers is an extension of the Bec et al. model [#Mercier_2010]_ and [#Mercier_2013]_. .. math:: \frac{1}{2 a_{\text{c},0} E^{'}} = \sum_{i=0}^{N} \frac{t_i}{\left(\pi a_{\text{c},i}^2 + 2t_ia_{\text{c},i}\right)E^{'}_{\text{f},i}} + \frac{1}{2\left(a_{\text{c},N} + \frac{2t_N}{\pi}\right)E^{'}_\text{s}} :label: mercier_Eeq .. math:: a_{\text{c},i+1} = a_{\text{c},i} + \frac{2t}{\pi} :label: mercier_ac With :math:`a_{\text{c},0}` equal to :math:`a_{\text{c}}`. Thus, the Young's modulus of the film can be calculated as : .. math:: E^{'}_{\text{f},0} = \left[\frac{\pi a_{\text{c},0}^2 + 2t_0a_{\text{c},0}}{t_0} \left[\frac{1}{2 a_{\text{c},0} E^{'}} - {\left(\sum_{i=1}^{N} {\frac{t_i}{\left(\pi a_{\text{c},i}^2 + 2t_ia_{\text{c},i}\right)E^{'}_{\text{f},i}} + \frac{1}{2\left(a_{\text{c},N} + \frac{2t_N}{\pi}\right)E^{'}_\text{s}}}\right)}\right]\right]^{-1} :label: mercier_Ef .. figure:: ./_pictures/multilayer_sample_elastic_model.png :scale: 30 % :align: center *Schematic of elastic multilayer model.* It is advised to perform nanoindentation tests on each layer of the multilayer sample, from the substrate up to the final stack of layers (see :numref:`multilayer_sample_methodology`). y successive iterations using the model of Mercier et al., values of Young's modulus of each layer are extracted from the contact stiffness. .. figure:: ./_pictures/multilayer_sample_methodology.png :name: multilayer_sample_methodology :scale: 30 % :align: center *Experimental process to apply for elastic multilayer model.* Find here the |matlab| function for the Mercier et al. model : `model_multilayer_elastic.m `_. Puchi-Cabrera et al. (2015) ----------------------------- Puchi-Cabrera et al. proposed in 2015, a description of the composite elastic modulus of multilayer coated systems [#PuchiCabrera_2015]_, based on the physically-based concept advanced by Rahmoun et al. [#Rahmoun_2009]_ : .. math:: \frac{1}{E} = \sum_{i=1}^{N} {\frac{x^\text{i}_\text{v}}{E^\text{i}_\text{f}} + \frac{x^\text{s}_\text{v}}{E_\text{s}}} :label: PuchiCabrera_model With :math:`x^\text{i}_\text{v}` and :math:`x^\text{s}_\text{v}`, respectively the volume fraction of each layer and the corresponding volume fraction of the substrate material. In his paper, Puchi-Cabrera modified and extended from the bilayer to the multilayer specimen, the models of Doerner and Nix, Gao, Bec, Menčík, Perriot and Barthel, Antunes, Korsunsky and Constantinescu and Bull. Plastic properties of a thin film on a substrate ################################################### Bückle (1961) ----------------- Bückle proposed an empirical law ("rule of thumb") to characterize thick coatings on a substrate [#Bückle_1961]_. It is possible to estimate empirically the hardness of the coating for indentation depth lower than ":math:`40\text{%}`" of the film thickness. Like the rule of the ":math:`10\text{%}`" for the extraction of elastic properties, it is a "rule of thumb". Because of the imperfections of the indenters, the roughness and the surface pollution, it is more meaningful to use this rule of the ":math:`40\text{%}`", for film thicker than :math:`500\text{nm}`. In the paper of Rother and Jehn [#Rother_1996]_, a good agreement with Bückle's "rule of thumb" is found for a range of indentation depths over coating thickness ratios between 1/10 (10%) and 1/7 (14.3%). According to Veprek-Heijman M.G.J. and Veprek S., the Bückle's rule is not valid in the case of superhard coatings the deformation of a softer substrate, such as steel [#Veprek_2016]_. To obtain correct measurements of the hardness of superhard coatings (H ≥ 40 GPa) their thickness must be sufficiently large (typically 5–7 μm or more on steel and about 3 μm on Si or cemented carbide). In the case of ultrahard coatings (H ≥ 80 GPa), the thickness should be even larger. In parallel, Bückle proposed an expression of the composite hardness :math:`H` in the case of a two-layer material with a weighted sum of the different layer hardnesses during indentation process [#Bückle_1961]_. .. math:: H = aH_\text{f} + bH_\text{s} :label: buckle_hardness With :math:`H_\text{f}` the hardness of the film, :math:`H_\text{s}` the hardness of the substrate and :math:`a + b = 1`. :math:`a` varies from :math:`1` when the hardness is not affected by the substrate, to :math:`0` when the indentation depth is approaching the film thickness. Kao and Byrne (1981) --------------------------- Kao and Byrne proposed the following model to describe the evolution of the composite hardness as a function of the reciprocal indentation depth [#Kao_1981]_, based on Bückle's model [#Bückle_1961]_ : .. math:: H \simeq H_\text{s} + 2 k_1 t_\text{f}\left(H_\text{f} - H_\text{s} \right)\frac{1}{h} :label: kao_hardness With :math:`k_1` a weighting factor of about 9%, independent of material characteristics. Find here the |matlab| functions for the model of Kao : - `model_kao.m `_. - `model_kao_film.m `_. Jönsson and Hogmark (1984) --------------------------- Jönsson and Hogmark used a simple geometrical approach based on a area "law of mixtures" to separate the substrate and film contributions to the measured hardness from Vickers indentation [#Jönsson_1984]_. .. math:: H = \frac{A_\text{f}}{A}H_\text{f} + \frac{A_\text{s}}{A}H_\text{s} :label: jonsson_hardness With :math:`A_\text{f}` the area on which the mean pressure :math:`H_\text{f}` acts and :math:`A_\text{s}` the area on which the mean pressure :math:`H_\text{s}` acts. The total area :math:`A` is the sum of :math:`A_\text{f}` and :math:`A_\text{s}` and the following expressions for the area ratios are given by Jönsson and Hogmark: .. math:: \frac{A_\text{f}}{A} = 2C\frac{t}{d} - C^2{\frac{t}{d}}^2 :label: jonsson_area_ratio_1 .. math:: \frac{A_\text{s}}{A} = 1 - \frac{A_\text{f}}{A} :label: jonsson_area_ratio_2 With :math:`d` the diagonal of the indent, :math:`t` the film thickness and :math:`C` a constant equal to :math:`0.5` for hard coatings on very soft substrates (:math:`6.3 < \frac{H_\text{f}}{H_\text{s}} < 12.9`) or to :math:`1` when the coatings and substrate hardnesses are more similar (:math:`1.8 < \frac{H_\text{f}}{H_\text{s}} < 2.3`). Burnett and Rickerby (1984) --------------------------- Burnett and Rickerby proposed afterwards a model based on a “volume law of mixtures” similar to Jönsson's relation, considering the volumes of the plastic zones, :math:`V_\text{f}` and :math:`V_\text{s}` respectively in the film and in the substrate [#Burnett1_1987]_ [#Burnett2_1987]_. .. math:: H = \frac{V_\text{f}}{V}H_\text{f} + \frac{V_\text{s}}{V}H_\text{s} :label: burnett_hardness With :math:`V = V_\text{f} + V_\text{s}`. Later, Iost and Bigot proposed to compare the predictions obtained using Jönsson-Hogmark and Burnett-Rickerby models and to understand why it is often reported that the Jönsson and Hogmark model does not hold for indentation prints less than the coating thickness. Methods were reviewed for calculating the composite hardness, and it was found that the simplifications made by the authors are not always valid. By taking in account all the terms of the equations of Jönsson and Hogmark, it was found that the relation between the hardness and the reciprocal length of the indentation print is not linear and depends on the ratio between the film thickness and the indentation print, as well as the variation of the hardness of the substrate and the film with the applied load. Comparison of the Burnett and Rickerby experimental data with the modified model led to very good agreement [#Iost_1996]_. Bhattacharya and Nix (1988) --------------------------- Bhattacharya and Nix proposed the following model from numerical simulations to extract the hardness of a coating [#Bhattacharya_1988]_ : .. math:: H = H_\text{s} + \left(H_\text{f} - H_\text{s} \right) e^{-\alpha\left(x\right)^n} :label: bhattacharya_hardness With :math:`x=h/t` and :math:`\alpha` and is an adjustable constants determined using the method of least squares. The variable :math:`n` is equal to 2 when the coating is harder than the substrate and to 1 when the substrate is harder than the coating. In the original paper, :math:`\alpha` is equal to :math:`\frac{H_\text{f}/H_\text{s}}{\sigma_\text{f}/\sigma_\text{s}\sqrt{E_\text{f}/E_\text{s}}}`, when the coating is harder than the substrate and to :math:`\frac{\sigma_\text{f}/\sigma_\text{s}}{E_\text{f}/E_\text{s}}`, when the substrate is harder than the coating. Find here the |matlab| functions for the model of Bhattacharya and Nix : `model_bhattacharya.m `_. “Volume law of mixtures” --------------------------- Several authors developed composite hardness model to calculate hardness of a coating using a volume law of mixture: - Sargent (1986) [#Sargent_1986]_ - Chicot and Lesage (1995) [#ChicotLesage_1995]_ - He et al. (1996) [#He_1996]_ - Korsunsky et al. (1996) [#Korsunsky_1998]_ - Fernandes et al. (2000) [#Fernandes_2000]_ - Tsui et al. (2003) [#Tsui_2003]_ Saha and Nix (2002) --------------------------- Based on the methodology proposed by Joslin and Oliver (1990) [#Joslin_1990]_ for a bulk material, extended to the coated system by Page et al. [#Page_1998]_, Saha and Nix proposed to use the following equation, giving the evolution of the hardness as a function of indentation depth, even when pile-up occurs [#Saha_2002]_: .. math:: H = {\beta}^2 \frac{4}{\pi} \frac{F_\text{c,max}}{S^2} (E*)^{2} :label: hardness_saha .. math:: \frac{1}{E^{*}} = \frac{1}{E_\text{i}^{'}} + \frac{1}{E_\text{f}^{'}} + \left(\frac{1}{E_\text{s}^{'}} - \frac{1}{E_\text{f}^{'}}\right) e^{-\alpha\left(x\right)} :label: doerner_nix_saha With :math:`x=(t-h_\text{c})/a_\text{c}` and :math:`E_\text{i}^{'}` the reduced Young's modulus of the indenter. This model was reused later by Han et al. [#Han_2005]_ and [#Han_2006]_. .. note:: This model is valid only for the case of elastically inhomogeneous film/substrate systems. Chen (2004) --------------------------- In this paper, Chen et al. proposed a modified volume fraction (CZ) model, in which a conical deformation volume shape is assumed, but the conical tip angle is undetermined. The conical tip angle contains the influences of pile-up or sink-in since the total depth is used in the modified CZ model. Due to the introduction of the conical angle, the modified CZ model could describe both a soft film on a hard substrate system and a hard film on a soft substrate system [#Chen_2004]_. Iost (2005) --------------------------- More recently, Iost et al. compared 8 different models in term of robustness, taking into account the indentation size effect on film and substrate proposed in the literature to estimate the hardness of coatings from indentations where both coating and substrate interact [#Iost_2005]_. Plastic properties of a thin film on a multilayer system ########################################################## Engel et al. (1992) --------------------------- Engel et al. proposed a simple method of interpreting the superficial (Vickers) hardness of multilayered specimen [#Engel1992]_, by expanding the concept of Jönsson and Hogmark [#Jönsson_1984]_ : .. math:: H = \frac{A_\text{s}}{A}H_\text{s} + \sum_{i=1}^{N}\frac{A_\text{f,i}}{A}H_\text{f,i} :label: engel_model .. math:: A = A_\text{s} + \sum_{i=1}^{N} A_\text{f,i} :label: engel_model_area With :math:`N` the number of layers deposited on the substrate, :math:`A_\text{f,i}` and :math:`H_\text{f,i}` respectively the flow pressure area and the hardness of an intermediate layer. According to the author, this model is applicable as long as the ratio of the film thickness :math:`t_\text{f}` over the indentation imprint size (i.e. the diagonal :math:`d` of square imprint in case of Vickers indentation) is less than 0.2, resulting in a parallel displacement of indenter, layers, and substrate. Rahmoun et al. (2009) --------------------------- The paper written by Ramhoun et al. is about the description of a multilayer model, developed to estimate coating hardness evolution in case of a multilayer specimen [#Rahmoun_2009]_. After some simplifcation, the following equations are established to take into account the hardness variation to the applied load and the crumbling of the coating. This model, based on the Jönsson and Hogmark's area law of mixtures, makes it possible to calculate the minimum in hardness value which depends on the hardness of the different phases, the thickness of the coatings and the compaction of the layers under the load transmitted by the Vickers indenter. From this model (which has some features in common with the Kao and Byrne approach [#Kao_1981]_), it can be shown that the substrate begins to influence the coating hardness at a critical penetration depth equal to the sum of interlayer and top coating thickness. .. math:: H = \frac{A_\text{f}}{A}H_\text{f} + \frac{A_\text{i}}{A}H_\text{i} + \frac{A_\text{s}}{A}H_\text{s} :label: rahmoun_model .. math:: \frac{A_\text{f}}{A} = \frac{2t_\text{f}}{h} - \frac{t_\text{f}^2}{h^2} :label: rahmoun_model_area_film .. math:: \frac{A_\text{i}}{A} = \frac{2t_\text{f}}{h} - \frac{t_\text{i}^2}{h^2} - 2\frac{t_\text{i}t_\text{f}}{h^2} :label: rahmoun_model_area_intermediate .. math:: \frac{A_\text{s}}{A} = 1 - \frac{2(t_\text{f}+t_\text{i})}{h} + \frac{(t_\text{f}+t_\text{i})^2}{h^2} :label: rahmoun_model_area_substrate With :math:`t_\text{f}` the coating thickness, :math:`t_\text{i}` the interlayer thickness, :math:`h` the indentation depth, :math:`A` the total contact area, :math:`A_\text{f}`, :math:`A_\text{i}` and :math:`A_\text{s}` respectively the coating, the interlayer and substrate area transmitting the mean contact pressure, and :math:`H_\text{f}`, :math:`H_\text{i}` and :math:`H\text{s}` respectively the coating, the interlayer and substrate hardness. Arrazat et al. (2010) --------------------------- For complex multilayer samples, Arrazat et al. proposed a model based on the volume law of mixture proposed by Sargent et al. [#Sargent_1986]_. In principle, the composite hardness is equal to the sum of intrinsic hardness of all materials plastically deformed, weighted by plastically deformed volume in the different layers [#Arrazat_2010]_. References ############# .. [#Argatov_2016] `Argatov I.I. and Sabina F.J., "Small-scale indentation of an elastic coated half-space: The effect of compliant substrate" (2016). `_ .. [#Arrazat_2010] `Arrazat B. et al., "Nano indentation de couches dures ultra minces de ruthénium sur or" (2010). `_ .. [#Bhattacharya_1988] `Bhattacharya A.K. and Nix W.D., "Analysis of elastic and plastic deformation associated with indentation testing of thin films on substrates" (1988). `_ .. [#Bec_2006] `Bec S. et al., "Improvements in the indentation method with a surface force apparatus" (2006). `_ .. [#Bückle_1961] Bückle H., "VDI Berichte" (1961). .. [#Bull_2005] `Bull S.J., "Nanoindentation of coatings" (2005). `_ .. [#Bull_2011] `Bull S.J., "Mechanical response of atomic layer deposition alumina coatings on stiff and compliant substrates" (2011). `_ .. [#Bull_2014] `Bull S.J., "A simple method for the assessment of the contact modulus for coated systems." (2014). `_ .. [#Burnett1_1987] `Burnett P.J. and Rickerby D.S., "The mechanical properties of wear-resistant coatings: I: Modelling of hardness behaviour." (1987). `_ .. [#Burnett2_1987] `Burnett P.J. and Rickerby D.S., "The mechanical properties of wear-resistant coatings: II: Experimental studies and interpretation of hardness." (1987). `_ .. [#ChenVlassak2001] `Chen X. and Vlassak J.J., "Numerical study on the measurement of thin film mechanical properties by means of nanoindentation." (2001). `_ .. [#Chen_2004] `Chen S. et al., "Nanoindentation of thin-film-substrate system: Determination of film hardness and Young's modulus" (2004). `_ .. [#ChicotLesage_1995] `Chicot D. and Lesage J., "Absolute hardness of films and coatings" (1995). `_ .. [#DoernerNix_1986] `Doerner M.F. and Nix W.D., "A method for interpreting the data from depth-sensing indentation instruments" (1986). `_ .. [#Engel1992] `Engel P.A. et al., "Interpretation of superficial hardness for multilayer platings*" (1992). `_ .. [#Fernandes_2000] `Fernandes J.V. et al., "A model for coated surface hardness" (2000). `_ .. [#FischerCripps_2004] `Fischer-Cripps, A.C., "Nanoindentation 3rd Ed." (2011) `_ .. [#Gao_1992] `Gao H. et al., "Elastic contact versus indentation modeling of multi-layered materials" (1992). `_ .. [#Han_2005] `Han S.M. et al., "Combinatorial studies of mechanical properties of Ti–Al thin films using nanoindentation" (2005). `_ .. [#Han_2006] `Han S.M. et al., "Determining hardness of thin films in elastically mismatched film-on-substrate systems using nanoindentation" (2006). `_ .. [#Hay_2011] `Hay J. and Crawford B., "Measuring substrate-independent modulus of thin films" (2011). `_ .. [#He_1996] `He J.L. et al., "Hardness measurement of thin films: Separation from composite hardness" (1996). `_ .. [#Hu_1998] `Hu X.Z. and Lawn B. R. "A simple indentation stress–strain relation for contacts with spheres on bilayer structures" (1998). `_ .. [#Iost_1996] `Iost A. and Bigot R., "Hardness of coatings" (1996). `_ .. [#Iost_2005] `Iost A. et al., "Dureté des revêtements : quel modèle choisir ?" (2005). `_ .. [#Jennett_2001] `Jennett N. M. and Bushby A. J., "Adaptive Protocol for Robust Estimates of Coatings Properties by Nanoindentation" (2001). `_ .. [#Jönsson_1984] `Jönsson B. and Hogmark S., "Hardness measurements of thin films" (1984). `_ .. [#Joslin_1990] `Joslin D.L. and Oliver W.C., "A new method for analyzing data from continuous depth-sensing microindentation tests" (1990). `_ .. [#Jung_2004] `Jung Y.-G. et al. "Evaluation of elastic modulus and hardness of thin films by nanoindentation" (2004). `_ .. [#Kao_1981] `Kao P.-W. and Byrne J. G., "Ion Implantation Effects on Fatigue and Surface Hardness" (1981). `_ .. [#King_1987] `King R.B., "Elastic analysis of some punch problems for a layered medium" (1987). `_ .. [#Korsunsky_1998] `Korsunsky A.M. et al. "On the hardness of coated system" (1998). `_ .. [#Korsunsky_2009] `Korsunsky A.M. and Constantinescu A., "The influence of indenter bluntness on the apparent contact stiffness of thin coatings" (2009). `_ .. [#Li_2010] `Li H. et al., "New methods of analyzing indentation experiments on very thin films" (2010). `_ .. [#Li_2015] `Li Y. et al., "Models for nanoindentation of compliant films on stiff substrates" (2015). `_ .. [#Menčík_1997] `Menčík J. et al., "Determination of elastic modulus of thin layers using nanoindentation" (1997). `_ .. [#Mercier_2010] `Mercier D. et al., "Young's modulus measurement of a thin film from experimental nanoindentation performed on multilayer systems" (2010). `_ .. [#Mercier_2013] `Mercier D., "Behaviour laws of materials used in electrical contacts for « flip chip » technologies" (2013). `_ .. [#Nix_1997] `Nix W.D., "Elastic and plastic properties of thin films on substrates: nanoindentation techniques" (1997). `_ .. [#Page_1998] `Page T.F. et al., "Nanoindentation Characterisation of Coated Systems: P:S2 - A New Approach Using the Continuous Stiffness Technique" (1998). `_ .. [#PaillerMattei_2008] `Pailler-Mattei C. et al., "In vivo measurements of the elastic mechanical properties of human skin by indentation tests" (2008). `_ .. [#Perriot_2004] `Perriot A. and Barthel E., "Elastic contact to a coated half-space: Effective elastic modulus and real penetration" (2004). `_ .. [#Pharr_1992] `Pharr G.M. and Oliver W.C., "Measurement of Thin Film Mechanical Properties Using Nanoindentation" (1992). `_ .. [#PuchiCabrera_2015] `Puchi-Cabrera. E.S. et al., "A description of the composite elastic modulus of multilayer coated systems" (2015). `_ .. [#Rahmoun_2009] `Rahmoun K. et al., "A multilayer model for describing hardness variations of aged porous silicon low-dielectric-constant thin films" (2009). `_ .. [#Rar_2002] `Rar A. et al., "Assessment of new relation for the elastic compliance of a film–substrate system." (2002). `_ .. [#Rother_1996] `Rother B. and Jehn H.A., "Effects of the substrate on the determination of thin film mechanical properties by coating and interface characterization by depth-sensing indentation experiments" (1996). `_ .. [#Saha_2002] `Saha R. and Nix W.D., "Effects of the substrate on the determination of thin film mechanical properties by nanoindentation" (2002). `_ .. [#Sargent_1986] `Sargent P.M., "Use of the Indentation Size Effect on Microhardness for Materials Characterization" (2002). `_ .. [#Simunkova_2003] `Šimůnková Š., "Mechanical properties of thin film–substrate systems" (2003). `_ .. [#Tsui_2003] `Tsui T.Y. et al., "A method for making substrate-independent hardness measurements of soft metallic films on hard substrates by nanoindentation" (2003). `_ .. [#VanVliet_2001] `Van Vliet K.J. and Gouldstone A., "Mechanical Properties of Thin Films Quantified Via Instrumented Indentation" (2001). `_ .. [#Veprek_2016] `Veprek-Heijman M.G.J. and Veprek S., "The deformation of the substrate during indentation into superhard coatings: Bückle's rule revised" (2016). `_ .. [#Xu_2006] `Xu H. and Pharr G.M., "An improved relation for the effective elastic compliance of a film/substrate system during indentation by a flat cylindrical punch." (2006). `_ .. [#Zhang2007] `Zhang F. et al., "Indentation of a hard film on a soft substrate: Strain gradient hardening effects" (2007). `_ .. [#Zhou2008] `Zhou X. et al., "Investigation on methods for dealing with pile-up errors in evaluating the mechanical properties of thin metal films at sub-micron scale on hard substrates by nanoindentation technique." (2008). `_