Nonlinear finite elements/Newton method for finite elements

Newton's method for nonlinear finite elements

The finite element approximation of a nonlinear boundary value problem results in a system of nonlinear equations of the form

$K(u)u=f$

Here, K is a nonlinear function relating the nodal degrees of freedom and the internal nodal forces and f is the externally applied nodal forces. Now the problem is to find the values of the nodal degrees of freedom such that the above equation is satisfied. In general the above equation does not have a close form solution (or some algorithm to find the closed form solution), hence we resort to numerical root finding methods to find an approximate solution. Newtons method is commonly used solution procedure for the solution of system simultaneous nonlinear equations. We rewrite the above equation in the form $g(x)=0$

$g(u):=K(u)u-f$

Then, the Newton iteration formula can be written as

$\begin{array}{cc}{u}_{r+1}& =u_r-[\frac{dg(u_r)}{du}{]}^{-1}g(u_r)\\ & =u_r-T^{-1}g(u_r)\end{array}$

The tangent stiffness is given by

$T=\frac{dg(u_r)}{du}$

The iterative procedure is terminated when either the residual is very small or the difference between successive solutions is less than a specified tolerance.

However, both the residual $𝐫$ and the solution $𝐮$ are vectors. We usually compare the $L_2$ (Euclidean) norm of the vectors with a tolerance $ϵ$. In symbolic form, we check the norm of the residual using

$\Vert 𝐫{\Vert }_2\le ϵ\text{or}\sqrt{𝐫\cdot 𝐫}\le ϵ\text{or}\sqrt{\sum _{i=1}^n{r}_i^2}\le ϵ.$

For the difference between successive solutions we check

$\frac{\Vert {𝐮}_{r+1}-{𝐮}_r{\Vert }_2}{\Vert {𝐮}_{r+1}{\Vert }_2}\le ϵ\text{or}\frac{\sqrt{({𝐮}_{r+1}-{𝐮}_r)\cdot ({𝐮}_{r+1}-{𝐮}_r)}}{\sqrt{{𝐮}_{r+1}\cdot {𝐮}_{r+1}}}\le ϵ\text{or}\frac{\sqrt{\sum _{i=1}^n(u_i^{r+1}-u_i^r{)}^2}}{\sqrt{\sum _{i=1}^n(u_i^{r+1}{)}^2}}\le ϵ.$

Template:Subpage navbar