Automated solution procedures for negotiating abrupt non‐linearities and branch points

In the inelastic stability analysis of plated structures, incremental‐iterative finite element methods sometimes encounter prohibitive solution difficulties in the vicinity of sharp limit points, branch points and other regions of abrupt non‐linearity. Presents an analysis system that attempts to trace the non‐linear response associated with these types of problems at minor computational cost. Proposes a semi‐heuristic method for automatic load incrementation, termed the adaptive arc‐length procedure. This procedure is capable of detecting abrupt non‐linearities and reducing the increment size prior to encountering iterative convergence difficulties. The adaptive arc‐length method is also capable of increasing the increment size rapidly in regions of near linear response. This strategy, combined with consistent linearization to obtain the updated tangent stiffness matrix in all iterative steps, and with the use of a “minimum residual displacement” constraint on the iterations, is found to be effective in avoiding solution difficulties in many types of severe non‐linear problems. However, additional procedures are necessary to negotiate branch points within the solution path, as well as to ameliorate convergence difficulties in certain situations. Presents a special algorithm, termed the bifurcation processor, which is effective for solving many of these types of problems. Discusses several example solutions to illustrate the performance of the resulting analysis system.

[1]  Jose Couto Marques,et al.  Stress computation in elastoplasticity , 1984 .

[2]  P. Stehlin,et al.  Postbuckling analysis using a general purpose code , 1986 .

[3]  Isaac Fried,et al.  Orthogonal trajectory accession to the nonlinear equilibrium curve , 1984 .

[4]  Peter Wriggers,et al.  A general procedure for the direct computation of turning and bifurcation points , 1990 .

[5]  M. Crisfield A FAST INCREMENTAL/ITERATIVE SOLUTION PROCEDURE THAT HANDLES "SNAP-THROUGH" , 1981 .

[6]  Ekkehard Ramm,et al.  Strategies for Tracing the Nonlinear Response Near Limit Points , 1981 .

[7]  J. C. Simo,et al.  Consistent tangent operators for rate-independent elastoplasticity☆ , 1985 .

[8]  Peter Wriggers,et al.  Consistent linearization for path following methods in nonlinear FE analysis , 1986 .

[9]  J. C. Simo,et al.  A J 2 -flow theory exhibiting a corner-like effect and suitable for large-scale computation , 1987 .

[10]  Donald W. White,et al.  Accurate and efficient nonlinear formulation of a nine-node shell element with spurious mode control , 1990 .

[11]  John W. Hutchinson,et al.  Analytical and Numerical Study of the Effects of Initial Imperfections on the Inelastic Buckling of a Cruciform Column , 1976 .

[12]  K. Bathe,et al.  ON THE AUTOMATIC SOLUTION OF NONLINEAR FINITE ELEMENT EQUATIONS , 1983 .

[13]  A. Chulya,et al.  An improved automatic incremental algorithm for the efficient solution of nonlinear finite element equations , 1987 .

[14]  Siu-Lai Chan Geometric and material non‐linear analysis of beam‐columns and frames using the minimum residual displacement method , 1988 .

[15]  Peter Wriggers,et al.  A simple method for the calculation of postcritical branches , 1988 .

[16]  Donald W. White,et al.  Large displacement formulation of a three-dimensional beam element with cross-sectional warping , 1992 .