# Equation solving methods for nonlinear calculations

### Technical Article

This article explains the equation solver for a nonlinear calculation with a Newton-Raphson iteration.

#### Introduction

The finite element method is always used when mechanical problems cannot be solved analytically. Often, nonlinear effects such as failure under compression (geometric nonlinearity), plastification (material nonlinearity), and contact or kinematic degrees of freedom are also taken into account. These effects, especially in the case of nonlinear material models, can be taken into account by means of an iterative calculation method.

#### FEM formulation

Basic steps of an FEM formulation (further information can be found in [1] ):

- Weak form of equilibrium

${\int}_{\mathrm{B}}\underset{\mathrm{nichtlinearer}\mathrm{Anteil}}{\underset{\u23df}{{\mathrm{grad\delta}}_{\mathrm{u}}\cdots \mathrm{\sigma dV}}}={\int}_{\mathrm{B}}{\mathrm{\delta}}_{\mathrm{u}}\xb7\mathrm{\rho bdV}+{\int}_{\partial \mathrm{B}}{\mathrm{\delta}}_{\mathrm{u}}\xb7{\mathrm{t}}_{\mathrm{o}}\mathrm{dA}$

**δ**_{u}Virtual (test) displacement **t**_{0}Initial load factor **σdV**Internal forces **ρbdV**Solid forces **B**Integrated area - Converting to Voigt Notation with Tensor 4. Level

${\int}_{\mathrm{B}}\mathrm{\delta \epsilon}\cdots \mathrm{\u2102}\cdots \mathrm{\epsilon dV}={\int}_{\mathrm{B}}\underset{\mathrm{Volumenkr\xe4fte}}{\underset{\u23df}{{\mathrm{\delta u}}^{\mathrm{T}}\mathrm{bdV}}}+{\int}_{\partial \mathrm{B}}\underset{\mathrm{Oberfl\xe4chenkr\xe4fte}}{\underset{\u23df}{{\mathrm{\delta u}}^{\mathrm{T}}\xb7{\mathrm{t}}_{0}\mathrm{dA}}}$

**C**Stiffness matrix **δε**Variation of the strain state **B**Integrated area **ε**Strain **U**Deformation

This notation is used to solve the approximate solution of the FE approach for nonlinear materials. - For this, the displacement field is multiplied by the application functions.

$\mathrm{u}(\mathrm{x},\mathrm{t})=\mathrm{H}\left(\mathrm{x}\right)\hat{\mathrm{u}}\left(\mathrm{t}\right)$

**u(x,t)**Displacement over time (load increment) **H**Form function **û**Nodal displacement - Using the derivative of the displacement in the weak form. Numerical integration is used to calculate the nodal displacement, and in post-processing, the stresses and strains are calculated using the material law.

#### Sequence of a Newton-Raphson iteration

Due to the nonlinear material behavior, the material matrix C in Equation 2 above changes with each strain step. The standard calculation method to solve this problem is the so-called Newton-Raphson iteration. It is used to linearize the function in a starting point. In the iteration, the stiffness matrix C of the step is always used. In the linearized iteration step, a tangent is placed at the zero point of the function.

The equations belonging to the flow diagram in the figure above:

- Break load into load steps.

${}^{\mathrm{t}+\u2206\mathrm{t}}\mathrm{f}_{\mathrm{ext}}={}^{\mathrm{t}}\mathrm{f}_{\mathrm{ext}}+\u2206\mathrm{f}$

- Predictor step

${}_{0}{}^{\mathrm{t}}\mathrm{K}{\u2206}_{0}\mathrm{\phi}={}^{\mathrm{t}+\u2206\mathrm{t}}\mathrm{f}_{\mathrm{ext}}-{}_{0}{}^{\mathrm{t}}\mathrm{f}_{\mathrm{int}}$

**K**Stiffness matrix of the previous time step ^{t+}Δ^{t}f_{ext}External forces increased by a further load step _{0}^{t}f_{int}Internal forces of the previous time step **ϕ**Strain - In point 3 Iteration of the flow diagram, the total distortion minus the plastic distortion is calculated (corrector step).

The goal of the iterative calculation is always that the sum of the loads results in zero. However, this is not possible numerically. Therefore, an abort rate ε is defined at which the calculation is aborted with sufficient accuracy.

$\mathrm{R}=\left|{}^{\mathrm{t}+\u2206\mathrm{t}}\mathrm{f}_{\mathrm{ext}}-{}_{\mathrm{n}}{}^{\mathrm{t}+\u2206\mathrm{t}}\mathrm{f}_{\mathrm{int}}\right|\mathrm{\epsilon}$

R |
Abortion rate |

f_{ext} |
External forces |

f_{int} |
internal forces |

ε |
Epsilon abortion rate |

t |
Time step |

In the program, you can control the termination rate under the calculation parameters.

The following figure shows the process of a Newton-Raphson iteration. In the first iteration

${}^{\mathrm{t}+\u2206\mathrm{t}}\mathrm{R}-{}^{\mathrm{t}+\u2206\mathrm{t}}\mathrm{F}^{0}$

R |
Abortion rate |

t |
Time step |

F |
Force |

the abort rate R or ε is not reached. Also in the second iteration (red), the tolerance limit is not reached. Only in the third iteration is the distance of the tangent stiffness so small that convergence is achieved.

As already mentioned, the deformation is continuously added up during the iteration.

#### Conclusion

The Newton-Raphson iteration has the consistency or convergence order 2. The number of "correct" locations of the iteration doubles for each step. As a result, a Newton-Raphson iteration converges quadratically and the accuracy increases with each iteration as the method converges. However, if the method does not converge, the error goes to infinity and the calculation is aborted.

The causes of the error are, for example, a too steep slope of the load-deformation curve and a too flat slope of the curve in the plastic area. If the load-deformation curve in the figure above shows a too strong buckling in the second iteration step, the material tangent and thus the stiffness matrix would be the slope of the elastic range. In this case, the slope of the zero point would be incorrect for the plastic area. This is one of the reasons why increasing the load steps is accompanied by improved convergence.

#### Keywords

Convergence Load step Time step

#### Reference

#### Links

### Write Comment...

### Write Comment...

#### Contact us

Do you have questions or need advice?

Contact our free e-mail, chat, or forum support or find various suggested solutions and useful tips on our FAQ page.

###### Recommended Events

###### Videos

###### Models to Download

###### Knowledge Base Articles

New

#### Modeling of a glass railing | Tips for entering data in RFEM

The architectural requirements for railings are still very high and usually require a high degree of transparency.

###### Screenshots

###### Product Features Articles

#### Automatic Determination of Number of Load Increments

If the check box "Number of load increments" is deactivated, the number of load increments will be determined automatically in RFEM to solve nonlinear tasks efficiently.

The used method is based on a heuristic algorithm.

###### Frequently Asked Questions (FAQ)

- Which Dlubal programs do I need to calculate membrane and textile structures?
- Does a user-defined layer composition in RF-LAMINATE need to be shared separately with other users along with the RFEM model file or is this information stored directly in the RFEM file?
- I have modeled two nodes on a surface in RFEM. The results do not show the displacement at these two nodes. How can I obtain it?
- Are the Dlubal Software programs suitable for the analysis and design of hydraulic steel structures?
- Which add-on modules and programs do I need for laminate and sandwich structures as well as cross-laminated timber (CLT)?

- I would like to analyze cable and tensile structures. Is it possible to use RFEM or RSTAB for this?
- How do I display the nodal displacements graphically?
- Are silos/tanks possible to export from RFEM into Revit?
- In connection with the calculation according to Theory III. Order, I get significantly smaller deformations than when calculating according to Theory I. or II. Order. How is this possible?
- Is it possible to also have RFEM output the area of each finite element along with its location?