Link Search Menu Expand Document
Theory Manual Version 3.6
 Section 3.1: Weak formulation for Solid Materials Up Section 3.1: Weak formulation for Solid Materials Subsection 3.1.2: Discretization 

3.1.1 Linearization

Equation (3.1-1) is the starting point for the nonlinear finite element method. It is highly nonlinear and any method attempting to solve this equation, such as the Newton-Raphson method, necessarily has to be iterative.
To linearize the finite element equations, the directional derivative of the virtual work in eq.(3.1-1) must be calculated. In an iterative procedure, the quantity will be approximated by a trial solution . Linearization of the virtual work equation around this trial solution gives The directional derivative of the virtual work will eventually lead to the definition of the stiffness matrix. In order to proceed, it is convenient to split the virtual work into an internal and external virtual work component: where and The result is listed here without details of the derivation – see [28] for details. The linearization of the internal virtual work is given by Notice that this equation is symmetric in and . This symmetry will, upon discretization, yield a symmetric tangent matrix.
The external virtual work has contributions from both body forces and surface tractions. The precise form of the linearized external virtual work depends on the form of these forces. For surface tractions, normal pressure forces may be represented in FEBio. The linearized external work for this type of traction is given by Discretization of this equation will also lead to a symmetric component of the tangent matrix.
FEBio currently supports gravity as a body force, . Since this force is independent of the geometry, the contribution to the linearized external work is zero. Another type of body force implemented in FEBio is the centrifugal force. For a body rotating with a constant angular speed , about an axis passing through the point and directed along the unit vector , the body force is given by , where is the vector distance from a point to the axis of rotation,
figure ../Figures/FigCentrifugalBodyForce.png
The resulting linearized external work is given by which produces a symmetric expression that will yield a symmetric matrix.
 Section 3.1: Weak formulation for Solid Materials Up Section 3.1: Weak formulation for Solid Materials Subsection 3.1.2: Discretization