The paper deals with approximations and the numerical realization of a class of hemivariational inequalities used for modeling of delamination and nonmonotone friction problems. Assumptions guaranteeing convergence of discrete models are verified and numerical results of several model examples computed by a nonsmooth variant of Newton method are presented.

In this paper we prove that every weak and strong local minimizer $u\in W^{1,2}(\Omega ,{ℝ}^3)$ of the functional $I\left(u\right)={\int }_{\Omega }{\left|Du\right|}^2+f\left(\mathrm{Adj}Du\right)+g\left(\det Du\right),$ where $u:\Omega \subset {ℝ}^3\to {ℝ}^3$, f grows like ${\left|\mathrm{Adj}Du\right|}^p$, g grows like ${\left|\det Du\right|}^q$ and 1<q<p<2, is $C^{1,\alpha }$ on an open subset ${\Omega }_0$ of Ω such that $\mathrm{𝑚𝑒𝑎𝑠}(\Omega \setminus {\Omega }_0)=0$. Such functionals naturally arise from nonlinear elasticity problems. The key point in order to obtain the partial regularity result is to establish an energy estimate of Caccioppoli type, which is based on an appropriate choice of the test functions. The limit case $p=q\le 2$ is also treated for weak local minimizers. ...