Product: ABAQUS/Standard
A major contribution to the computational effort involved in nonlinear analysis is the solution of the nonlinear equations (Equation 2.2.1–1). In most cases ABAQUS/Standard uses Newton's method to solve these equations, as described in Nonlinear solution methods in ABAQUS/Standard, Section 2.2.1. The principal advantage of Newton's method is its quadratic convergence rate when the approximation at iteration i is within the “radius of convergence”—that is, when the gradients defined by provide an improvement to the solution. The method has two major disadvantages: the Jacobian matrix has to be calculated, and this same matrix has to be solved. The calculation of the Jacobian matrix is a problem because, in many important cases, it is difficult to derive the form of the matrix algebraically. The solution of the Jacobian is a problem because of the computational effort involved: as the problem size increases, the direct solution of the linear equations can dominate the entire computational effort.
There are a number of important nonlinear applications where the Jacobian is symmetric, is fairly well conditioned, and does not change greatly from one iteration to the next. Examples are implicit dynamic time integration with small time increments relative to the periods of the natural vibrations that participate in the response or small-displacement elastic-plastic analysis where the yielding is confined (such as occurs in many practical fracture mechanics applications). In such cases, especially when the problem is large, it can be less expensive to use an alternative to the Newton approach to the solution of the nonlinear equations. The “quasi-Newton” methods are such an approach; and Matthies and Strang (1979) have shown that, for systems of equations with a symmetric Jacobian matrix, the BFGS (Broyden, Fletcher, Goldfarb, Shanno) method can be written in a simple form that is especially effective on the computer and is successful in such applications. This method is implemented in ABAQUS/Standard and is described in this section. The user must select this method explicitly: by default, ABAQUS/Standard uses the standard Newton method.
The basis of quasi-Newton methods is to obtain a series of improved approximations to the Jacobian matrix, , that satisfy the secant condition:
so that approaches as the iterations proceed. Equation 2.2.2–1 is the basic quasi-Newton equation.For convenience we define the change in the residual from one iteration to the next as
Matthies and Strang's implementation of the BFGS method is a computationally inexpensive method of creating a series of approximations to that satisfy Equation 2.2.2–1 and retain the symmetry and positive definiteness of . They accomplish this by updating to using a “product plus increment” form:
whereIn the actual implementation of this version of the BFGS method, each is not stored: rather, a “kernel” matrix, , is used (as the decomposition of ), and the update is accomplished by premultiplication of the kernel matrix by the terms
In general, the rate of convergence of the quasi-Newton method is slower than the quadratic rate of convergence of Newton's method, though faster than the linear rate of convergence of the modified Newton method.