# FEM model¶

## Finite Element Modelling (FEM) of conical indentation¶

The present model is a simulation of the conical nanoindentation process, using the FEM software ABAQUS.

The Matlab function used to generate a Python script for ABAQUS is : python_abaqus

The model is axisymmetric with a geometry dependent mesh and restricted boundaries conditions.

### Geometry of the (multilayer) sample¶

Each layer of the sample is characterized by its thickness (\(t_\text{i}\)).

The thickness of the substrate (\(t_\text{sub}\)) is set as 2 times the highest thin film thickness.

The width (\(w\)) of the sample is calculated as a function of the substrate thickness or the indenter tip defect.

No delamination is allowed between thin films or between thin film and substrate.

Note

Dimensions are in nm.

### Geometry of the indenter¶

The indenter is defined as a rigid cono-spherical indenter. A spherical part is defined at the apex of the conical indenter (see Figure 32).

The radius \(R\) of the spherical part is calculated from the tip defect \(h_\text{tip}\) and the cone angle \(\alpha\), using the following equation. For Berkovich, Vickers and Cube-Corner indenters, the equivalent cone angle is used to set the cone angle.

\[R = \frac{h_\text{tip}}{\frac{1}{sin(\alpha)} - 1}\]

In case of a perfect conical indenter (\(h_\text{tip} = 0\text{nm}\)), a tip defect of 1nm giving a radius of \(R = 1.6\text{nm}\) is set into the Python file. Defining a cono-spherical tip avoids the geometrical singularity at the apex of the perfect conical indenter, which would imply an infinite stress at the contact interface. But, any non-zero positive value for the tip radius can be set into the GUI.

The transition depth \(h_\text{trans}\) between the spherical and the conical parts of a cono-spherical indenter, is calculated from the followig equation:

\[h_\text{trans} = R(1-sin(\alpha))\]

Find here the Matlab function to calculate the tip radius: tipRadius.m.

Find here the Matlab function to calculate the transition depth: transitionDepth.m.

Find here the Matlab function to calculate the tip defect: tipDefect.m.

### Mesh¶

The multilayer sample is divided by default into solid elements with eight nodes and axisymmetric deformation element CAX8R is adopted.

It is possible to divide the sample into solid elements with four nodes and with axisymmetric deformation element CAX4R, by changing the value of the variable “linear_elements” in the Matlab function python4abaqus from 0 (quadratic elements) to 1 (linear elements).

Note

- CAX4R: A 4-node bilinear axisymmetric quadrilateral, reduced integration, hourglass control.
- CAX3: A 3-node linear axisymmetric triangle.
- CAX8R: An 8-node biquadratic axisymmetric quadrilateral, reduced integration.
- CAX6M: A 6-node modified quadratic axisymmetric triangle.

### Material properties¶

For each layers of the multilayer sample, the material properties (Young’s modulus and Poisson’s ratio) are defined using the inputs given by the user from the GUI. Material properties are considered by default to be isotropic. The density is set by default to 1.0.

Note

Young’s moduli are in GPa.

### Contact definition¶

The contact is defined by default frictionless for the tangential behavior and hard for the normal behavior.

The external surface of the indenter is defined as the “master” region and the top surface of the (multilayer) sample is defined as the “slave” region.

### Boundaries conditions¶

Nodes are constrained along the rotation axis from moving in the radial direction (\(x\)). The nodes on the bottom surface of the sample are constrained along the radial axis from moving in the radial (\(x\)) and vertical (\(z\)) directions (see Figure 34 and Figure 35).

Indentation process is simulated by imposing a vertical displacement to the rigid indenter along the (\(z\)) axis (see Figure 34 and Figure 35). A value of 200nm for the indentation depth is set by default.

Warning

Indentation displacement is given in nanometers and is negative.

## Generation of the Python script for ABAQUS¶

After material properties are configured (Young’s moduli and Poisson’s ratios) and the model geometry is given (thickness for each thin films), a Python script for ABAQUS can be generated by pressing the ‘FEM’ button.

The python script is saved in the folder where your nanoindentation results are stored.

To generate the FEM model in ABAQUS, apply the following procedure:

- start ABAQUS
- select the folder containing input files : ‘File’ ==> ‘Set Work Directory…’
- select and run the Python file containing the FEM model (*.py) : File’ ==> ‘Run Script’

Note

Dimensions are in nm and Young’s moduli are in GPa, implying that load is in nN.

## References¶

[1] | Atkins A.G. and Tabor D., “Plastic indentation in metals with cones” (1965). |

[2] | DiCarlo A. et al., “Prediction of stress–strain relation using cone indentation: effect of friction” (2004). |

[3] | Harsono E. et al., “The effect of friction on indentation test results” (2008). |

[4] | Johnson K.L., “Contact Mechanics” (1987), ISBN - 9780521347969. |

[5] | Wang. Y., “Effects of indenter angle and friction on the mechanical properties of film materials” (2016). |

- ABAQUS documentation
- Charleux L., “Abapy Documentation”
- Bhattacharya A.K. and Nix W.D., “Analysis of elastic and plastic deformation associated with indentation testing of thin films on substrates” (1988).
- Bhattacharya A.K. and Nix W.D., “Finite element analysis of cone indentation.” (1991).
- Madsen D.T. et al. “Finite Element Simulation of Indentation Behavior of thin Films.” (1991).
- Sun Y. et al., “Finite element analysis of the critical ratio of coating thickness to indentation depth for coating property measurements by nanoindentation.” (1995).
- Bolshakov A. and Pharr G.M., “Influences of pile-up on the measurement of mechanical properties by load and depth sensing indentation techniques.” (1998)
- Lichinchi M. et al., “Simulation of Berkovich nanoindentation experiments on thin films using finite element method.” (1998).
- Chen X. and Vlassak J.J., “A Finite Element Study on the Nanoindentation of Thin Films.” (2000).
- Chen X. and Vlassak J.J., “Numerical study on the measurement of thin film mechanical properties by means of nanoindentation.” (2001).
- Bouzakis K.-D. et al., “Thin hard coatings stress–strain curve determination through a FEM supported evaluation of nanoindentation test results.”, (2001).
- Hubert N. et al. “Identification of elastic-plastic material parameters from pyramidal indentation of thin films.” (2002).
- Shan Z. and Sitaraman S.K., “Elastic–plastic characterization of thin films using nanoindentation technique.” (2003).
- He J.L. and Veprek S., “Finite element modeling of indentation into superhard coatings.” (2003).
- Xu Z.-H. and Rowcliffe D., “Finite element analysis of substrate effects on indentation behaviour of thin films.” (2004).
- Panich N. and Sun Y. “Effect of penetration depth on indentation response of soft coatings on hard substrates: a finite element analysis” (2004)
- Cai X. and Bangert H., “Hardness measurements of thin films-determining the critical ratio of depth to thickness using FEM.” (2005).
- Bressan J.D. et al., “Modeling of nanoindentation of bulk and thin film by finite element method.” (2005).
- Pelletier H. et al., “Characterization of mechanical properties of thin films using nanoindentation test.” (2006)
- Charleux L., “Micromécanique de l’essai d’indentation : expériences et simulations.”, PhD thesis, Institut National Polytechnique de Grenoble (2006).
- Xu H., “A Finite Element Study of the Contact Stiffness of Homogenous Materials and Thin Films.” PhD thesis - University of Tennessee - Knoxville (2007)
- Antunes J.M. et al., “On the determination of the Young’s modulus of thin films using indentation tests” (2007).
- Chen S.H. et al., “Small scale, grain size and substrate effects in nano-indentation experiment of film–substrate systems.” (2007).
- Huang X. and Pelegri A.A., “Finite element analysis on nanoindentation with friction contact at the film/substrate interface.” (2007).
- Kopernik M. and Pietrzyk M., “2d numerical simulation of elasto-plastic deformation of thin hard coating systems in deep nanoindentation test with sharp indenter.” (2007).
- Pelegri A.A. and Huang X., “Nanoindentation on soft film/hard substrate and hard film/soft substrate material systems with finite element analysis.” (2008).
- Wittler O. et al., “Mechanical characterisation of thin metal layers by modelling of the nanoindentation experiment.” (2008).
- Sakharova N.A. et al., “Comparison between Berkovich, Vickers and conical indentation tests: A three-dimensional numerical simulation study.” (2009).
- Rao C.S. and Reddy C.E., “Finite element modeling of nanoindentation to extract load-displacement characteristics of bulk materials and thin films.” (2009).
- Gonzalez M. et al., “Modeling the substrate effects on nanoindentation mechanical property measurement.” (2009).
- Chen C., “2-D Finite Element Modeling for Nanoindentation and Fracture Stress Analysis.”, PhD thesis - University of South Florida (2009).
- Moore S.W. et al., “Nanoindentation in elastoplastic materials: insights from numerical simulations.” (2010).
- Guo W., “On the influence of indenter tip geometry on the identification of material parameters in indentation testing.”, PhD thesis, University of Liège (2010).
- Wang J.S. et al., “Identification of elastic parameters of transversely isotropic thin films by combining nanoindentation and FEM analysis.” (2010).
- Jiang L.M. et al., “Elastic-plastic properties of thin film on elastic-plastic substrates characterized by nanoindentation test.” (2010).
- Miller D.C. et al., “Thermo-mechanical properties of alumina films created using the atomic layer deposition technique.” (2010).
- Dowhan L. et al., “Investigation of thin films by nanoindentation with doe and numerical methods.” (2011).
- Isselé H. et al., “Determination of the Young’s modulus of a TiN Thin Film by nanoindentation: analytical models and numerical FEM simulation.” (2012).
- Phiciato’s blog (2013)
- Karimpour M. et al., “An inverse method for extracting the mechanical properties of the constituent materials of a multilayer from nanoindentation data.” (2013).
- Fourcade T., “Études des propriétés mécaniques de matériaux métalliques en couches minces”, PhD thesis - Université de Toulouse (2013).
- Moćko W. et al., “Simulation of nanoindentation experiments of single-layer and double-layer thin films using finite element method.” (2014).
- Kibech S. et al., “Nanoindentation” (2014).
- Kopernik M. and Milenin A., “Numerical modeling of substrate effect on determination of elastic and plastic properties of TiN nanocoating in nanoindentation test.” (2014).
- Li Y. et al., “Models for nanoindentation of compliant films on stiff substrates” (2015).
- Moreau S. et al., “Application of Design of Computer Experiments (DoCE) method for the extraction of the elasto-plastic behavior law of ECD copper through nano-indentation tests” (2015).
- Gupta A. K. et al., “Evaluation of elasto-plastic properties of ITO film using combined nanoindentation and finite element approach” (2015).
- Veprek-Heijman M.G.J. and Veprek S., “The deformation of the substrate during indentation into superhard coatings: Bückle’s rule revised” (2016).
- Du Y. et al., “A novel tri-layer nanoindentation method to measure the mechanical properties of a porous brittle ultra-low-k dielectric thin film” (2017).
- Lofaj F. and Németh D., “The effects of tip sharpness and coating thickness on nanoindentation measurements in hard coatings on softer substrates by FEM” (2017).
- Chakroun N. et al., “A new inverse analysis method for identifying the elastic properties of thin films considering thickness and substrate effects simultaneously” (2017).
- Chakroun N. et al., “Measuring elastic properties of the constituent multilayer coatings for different modulation periods” (2018).
- Pöhl F., “Determination of unique plastic properties from sharp indentation” (2019).