# 6 Oliver Pharr with ODR

For a definition of the Oliver Pharr method see [11].
This method is shown unless the software was compiled without Fortran support.

## 6.1 Window

The window consists of several blocks:

• Info displays the maximum depth and force during the indentation

• Parameters shows the selected range in nm and in % of the maximum force, the correction $\beta$, and the radial correction.

• -

The fitting range can be selected either using the mouse or typing in the range entries. The range can be defined either in nm or in percent of the maximum force. It is often recommended to use the range 40–98 % F${}_{\mathrm{max}}$ for the fit, see section 6.2.

• -

The parameter $\beta$ accounts for any deviations from the axisymmetric case and is used in the calculation of the reduced modulus in equation (10). Currently, the default value is no correction $\beta=1.0$. The value supplied by the user is saved in the settings and can be reset to its default value.

• -

Optional radial correction of calculated hardness and modulus values according to ISO 14577-1:2015. The angle $\alpha$ denotes the cone-equivalent angle between the tip side and axis. For a Berkovich tip, $\alpha=70.3^{\circ}$.

• Set defaults button resets the fitting range, the beta value, and radial correction settings according to recommendations provided in ISO 14577-1.

• Fit button, see section 6.2 for details of the calculation.

• Results displays all results in the following order: the residual depth $h_{\mathrm{p}}$, the power $m$ of the power law function, the parameter $\varepsilon$, the contact depth $h_{\mathrm{c}}$, the slope $S$, the contact area $A_{\mathrm{p}}(h_{\mathrm{c}})$, the indentation hardness $H_{IT}$, the contact modulus $E_{r}$, the indentation modulus $E_{IT}$ and the ranges used for the fitting procedure. The variables are described in detail in section 6.2. If the fittings procedure failed a warning is shown.

• Uncertainties show the uncertainty analysis window, see section 15.2.1.

• Show log Show the report about the fitting procedure in a separate window. The reports are saved to files fit.log.op.err and fit.log.op.rpt.

• Save save parameters and results to given file.

• Graph display the unloading curve and the fitted curves. Stepwise zooming/unzooming can be performed by selecting a range with the mouse and pressing the Zoom/ Unzoom buttons. The graph is restored to its original size by the Restore button.

## 6.2 Procedure

This is a slight modification of the standard Oliver Pharr method described in section 5.2 using a better fitting procedure.

1. 1.

Fit the upper part of the unloading curve with a power law function

 $F=\alpha(h-h_{\mathrm{p}})^{m}.$

using orthogonal least squares as implemented in the package ODRPACK95 [13]. The range should be approx. 40–98 % F${}_{\mathrm{max}}$. All three parameters are fitted.

2. 2.

Same as steps 3–4 in 5.2.

3. 3.

If the radial correction is on, the hardness and reduced modulus are found by iterating the following relations

 $\displaystyle H^{n}$ $\displaystyle=$ $\displaystyle H^{0}\left(1+K\frac{H^{n-1}}{E_{IT}^{n-1}}\right)^{2}$ (12) $\displaystyle E_{r}^{n}$ $\displaystyle=$ $\displaystyle E_{r}^{0}\left(1+K\frac{H^{n-1}}{E_{IT}^{n-1}}\right)$ (13)

where $H^{0}$ and $E_{r}^{0}$ are the initial estimates obtained from step 4 in section 5.2. The iteration limit is set to a relative difference $0.5\cdot 10^{-4}$ Young’s modulus is updated in each step according to (11) and $K$ is the radial correction factor computed from Poisson’s ratio and the angle between tip sides and axis $\alpha$

 $K=\frac{1-2\nu}{2(1+\nu)}\sin\alpha.$ (14)

The area function is corrected with the final values of hardness and modulus as

 $A=A^{0}\left(1+K\frac{H}{E_{IT}}\right)^{2}.$ (15)