A modified Newton‐type Koiter‐Newton method for tracing the geometrically nonlinear response of structures

The Koiter‐Newton (KN) method is a combination of local multimode polynomial approximations inspired by Koiter's initial postbuckling theory and global corrections using the standard Newton method. In the original formulation, the local polynomial approximation, called a reduced‐order model, is used to make significantly more accurate predictions compared to the standard linear prediction used in conjunction with Newton method. The correction to the exact equilibrium path relied exclusively on Newton‐Raphson method using the full model. In this paper, we proposed a modified Newton‐type KN method to trace the geometrically nonlinear response of structures. The developed predictor‐corrector strategy is applied to each predicted solution of the reduced‐order model. The reduced‐order model can be used also in the correction phase, and the exact full nonlinear model is applied only to calculate force residuals. Remainder terms in both the displacement expansion and the reduced‐order model are well considered and constantly updated during correction. The same augmented finite element model system is used for both the construction of the reduced‐order model and the iterations for correction. Hence, the developed method can be seen as a particular modified Newton method with a constant iteration matrix over the single KN step. This significantly reduces the computational cost of the method. As a side product, the method has better error control, leading to more robust step size adaptation strategies. Numerical results demonstrate the effectiveness of the method in treating nonlinear buckling problems.     

A new predictor–corrector approach for the numerical integration of coupled electromechanical equations

In this paper, a new approach for the numerical solution of coupled electromechanical problems is presented. The structure of the considered problem consists of the low‐frequency integral formulation of the Maxwell equations coupled with Newton–Euler rigid‐body dynamic equations. Two different integration schemes based on the predictor–corrector approach are presented and discussed. In the first method, the electrical equation is integrated with an implicit single‐step time marching algorithm, while the mechanical dynamics is studied by a predictor–corrector scheme. The predictor uses the forward Euler method, while the corrector is based on the trapezoidal rule. The second method is based on the use of two interleaved predictor–corrector schemes: one for the electrical equations and the other for the mechanical ones. Both the presented methods have been validated by comparison with experimental data (when available) and with results obtained by other numerical formulations; in problems characterized by low speeds, both schemes produce accurate results, with similar computation times. When high speeds are involved, the first scheme needs shorter time steps (i.e., longer computation times) in order to achieve the same accuracy of the second one. A brief discussion on extending the algorithm for simulating deformable bodies is also presented. An example of application to a two‐degree‐of‐freedom levitating device based on permanent magnets is finally reported. Copyright © 2015 John Wiley & Sons, Ltd.     

Postbuckling analysis stabilized by penalty springs and intermediate corrections

Whenever a critical point in a non-linear finite element analysis is reached, an implicit Newton procedure is prevented from proceeding until the stiffness matrix is stabilized. Most stabilization procedures result in a damped Newton scheme. This can cause a reduced convergence rate. Based on documented stabilization strategies, an iterative procedure to reduce negative impact on the convergence rate resulting from the damping is to be proposed here. This will be done by a corrector iteration carried out between two successive equilibrium iterative steps. The equilibrium iteration is a Newton–Raphson scheme. Since the stiffness matrix has already been factorized within the preceding standard equilibrium iterative step, the stiffness matrix will be held constant during the corrector iteration, which in turn allows for a computational efficient treatment of the additional iterations. The proposed procedure has been strongly driven by Wright and Gaylord’s (ASCE J. Struct. Div., 94, 1143–1163, 1968) investigation.     

A heterogeneous space–time full approximation storage multilevel method for molecular dynamics simulations

A heterogeneous space–time full approximation storage (HFAS) multilevel formulation for molecular dynamics simulations is developed. The method consists of a waveform Newton smoothing that produces initial space–time iterates and a coarse model correction. The formulation is coined as heterogeneous since it permits different interatomic potentials to be applied at different physical scales. This results in a flexible framework for physics coupling. Time integration is performed in windows using the implicit Newmark predictor–corrector method that permits larger time integration steps than the explicit method. The size of the time steps is governed by accuracy rather than by stability considerations of the algorithm. We study three different variants of the method: the Picard iteration, constrained dynamics and force splitting. Numerical examples show that FAS based on force splitting provides significant time savings compared to standard explicit methods and alternative implicit space–time schemes. Parallel studies of the Picard iteration on harmonic problems illustrate the time parallelization effect that leads to a superior parallel performance compared to explicit methods. Copyright © 2007 John Wiley & Sons, Ltd.     

A Koiter‐Newton approach for nonlinear structural analysis

A new approach termed the Koiter‐Newton is presented for the numerical solution of a class of elastic nonlinear structural response problems. It is a combination of a reduction method inspired by Koiter's post‐buckling analysis and Newton arc‐length method so that it is accurate over the entire equilibrium path and also computationally efficient in the presence of buckling. Finite element implementation based on element independent co‐rotational formulation is used. Various numerical examples of buckling sensitive structures are presented to evaluate the performance of the method. The examples demonstrate that the method is robust and completely automatic and that it outperforms traditional path‐following techniques. This improved efficiency will open the door for the direct use of detailed nonlinear finite element models in the design optimization of next generation flight and launch vehicles. Copyright © 2013 John Wiley & Sons, Ltd.     

An algorithm for the matrix‐free solution of quasistatic frictional contact problems

A contact enforcement algorithm has been developed for matrix‐free quasistatic finite element techniques. Matrix‐free (iterative) solution algorithms such as non‐linear conjugate gradients (CG) and dynamic relaxation (DR) are desirable for large solid mechanics applications where direct linear equation solving is prohibitively expensive, but in contrast to more traditional Newton–Raphson and quasi‐Newton iteration strategies, the number of iterations required for convergence is typically of the same order as the number of degrees of freedom of the model. It is therefore crucial that each of these iterations be inexpensive to per‐form, which is of course the essence of a matrix free method. In applying such methods to contact problems we emphasize here two requirements: first, that the treatment of the contact should not make an average equilibrium iteration considerably more expensive; and second, that the contact constraints should be imposed in such a way that they do not introduce spurious energy that acts against the iterative solver. These practical concerns are utilized to develop an iterative technique for accurate constraint enforcement that is suitable for non‐linear conjugate gradient and dynamic relaxation iterative schemes. Copyright © 1999 John Wiley & Sons, Ltd.     

Efficient non‐linear model reduction via a least‐squares Petrov–Galerkin projection and compressive tensor approximations

A Petrov–Galerkin projection method is proposed for reducing the dimension of a discrete non‐linear static or dynamic computational model in view of enabling its processing in real time. The right reduced‐order basis is chosen to be invariant and is constructed using the Proper Orthogonal Decomposition method. The left reduced‐order basis is selected to minimize the two‐norm of the residual arising at each Newton iteration. Thus, this basis is iteration‐dependent, enables capturing of non‐linearities, and leads to the globally convergent Gauss–Newton method. To avoid the significant computational cost of assembling the reduced‐order operators, the residual and action of the Jacobian on the right reduced‐order basis are each approximated by the product of an invariant, large‐scale matrix, and an iteration‐dependent, smaller one. The invariant matrix is computed using a data compression procedure that meets proposed consistency requirements. The iteration‐dependent matrix is computed to enable the least‐squares reconstruction of some entries of the approximated quantities. The results obtained for the solution of a turbulent flow problem and several non‐linear structural dynamics problems highlight the merit of the proposed consistency requirements. They also demonstrate the potential of this method to significantly reduce the computational cost associated with high‐dimensional non‐linear models while retaining their accuracy. Copyright © 2010 John Wiley & Sons, Ltd.     

An improved numerical scheme for the sine‐Gordon equation in 2+1 dimensions

A rational approximant of order 4, which is applied to a three‐time‐level recurrence relation, is used to transform the initial/boundary‐value problem associated with the two‐dimensional sine‐Gordon (SG) equation arising in the Josephson junctions problem. The resulting non‐linear system, which is analyzed for stability, is solved using an appropriate predictor–corrector (P–C) scheme, in which an explicit scheme of order 2 is used as predictor. For the implementation of the corrector, in order to avoid extended matrix evaluations, an auxiliary vector was successfully introduced. In this P–C scheme, a modification in the corrector has been proposed according to which the already evaluated corrected values are considered. The behavior of this P–C scheme is tested numerically on line and ring solitons known from the bibliography regarding the SG equation and conclusions for both undamped and damped problems are derived. Copyright © 2008 John Wiley & Sons, Ltd.     

Three‐dimensional finite element computations for frictional contact problems with non‐associated sliding rule

This paper presents an algorithm for solving anisotropic frictional contact problems where the sliding rule is non‐associated.The algorithm is based on a variational formulation of the complex interface model that combine the classical unilateral contact law and an anisotropic friction model with a non‐associated slip rule. Both the friction condition and the sliding potential are elliptical and have the same principal axes but with different semi‐axes ratio. The frictional contact law and its inverse are derived from a single non‐differentiable scalar‐valued function, called a bi‐potential. The convexity properties of the bi‐potential permit to associate stationary principles with initial/boundary value problems. With the present formulation, the time‐integration of the frictional contact law takes the form of a projection onto a convex set and only one predictor–corrector step addresses all cases (sticking, sliding, no‐contact). A solution algorithm is presented and tested on a simple example that shows the strong influence of the slip rule on the frictional behaviour. Copyright 2004 John Wiley & Sons, Ltd.     

LATIN: A new view and an extension to wave propagation in nonlinear media

The LATIN (acronym of LArge Time INcrement) method was originally devised as a non‐incremental procedure for the solution of quasi‐static problems in continuum mechanics with material nonlinearity. In contrast to standard incremental methods like Newton and modified Newton, LATIN is an iterative procedure applied to the entire loading path. In each LATIN iteration, two problems are solved: a local problem, which is nonlinear but algebraic and miniature, and a global problem, which involves the entire loading process but is linear. The convergence of these iterations, which has been shown to occur for a large class of nonlinear problems, provides an approximate solution to the original problem. In this paper, the LATIN method is presented from a different viewpoint, taking advantage of the causality principle. In this new view, LATIN is an incremental method, and the LATIN iterations are performed within each load step, similarly to the way that Newton iterations are performed. The advantages of the new approach are discussed. In addition, LATIN is extended for the solution of time‐dependent wave problems. As a relatively simple model for illustrating the new formulation, lateral wave propagation in a flat membrane made of a nonlinear material is considered. Numerical examples demonstrate the performance of the scheme, in conjunction with finite element discretization in space and the Newmark trapezoidal algorithm in time. Copyright © 2017 John Wiley & Sons, Ltd.     

Incremental elastoplastic analysis for active macro‐zones

In this paper a strategy to perform incremental elastoplastic analysis using the symmetric Galerkin boundary element method for multidomain type problems is shown. The discretization of the body is performed through substructures, distinguishing the bem‐elements characterizing the so‐called active macro‐zones, where the plastic consistency condition may be violated, and the macro‐elements having elastic behaviour only. Incremental analysis uses the well‐known concept of self‐equilibrium stress field here shown in a discrete form through the introduction of the influence matrix (self‐stress matrix). The nonlinear analysis does not use updating of the elastic response inside each plastic loop, but at the end of the load increment only. This is possible by using the self‐stress matrix, both, in the predictor phase, for computing the stress caused by the stored plastic strains, and, in the corrector phase, for solving a nonlinear global system, which provides the elastoplastic solution of the active macro‐zones. The use of active macro‐zones gives rise to a nonlocal and path‐independent approach, which is characterized by a notable reduction of the number of plastic iterations. The proposed strategy shows several computational advantages as shown by the results of some numerical tests, reported at the end of this paper. These tests were performed using the Karnak.sGbem code, in which the present procedure was introduced as an additional module.Copyright © 2012 John Wiley & Sons, Ltd.     

Adaptive time stepping strategy for solidification processes based on modified local time truncation error

Adaptive time step methods provide a computationally efficient means of simulating transient problems with a high degree of numerical accuracy. However, choosing appropriate time steps to model the transient characteristics of solidification processes is difficult. The Gresho–Lee–Sani predictor–corrector strategy, one of the most commonly applied adaptive time step methods, fails to accurately model the latent heat release associated with phase change due to its exaggerated time steps while the apparent heat capacity method is applied. Accordingly, the current study develops a modified local time truncation error‐based strategy designed to adaptively adjust the size of the time step during the simulated solidification procedure in such a way that the effects of latent heat release are more accurately modeled and the precision of the computational solutions correspondingly improved. The computational accuracy and efficiency of the proposed method are demonstrated via the simulation of several one‐dimensional and two‐dimensional thermal problems characterized by different phase change phenomena and boundary conditions. The feasibility of the proposed method for the modeling of solidification processes is further verified via its applications to the enthalpy method. Copyright © 2009 John Wiley & Sons, Ltd.     

Pseudo‐transient continuation for nonlinear transient elasticity

This paper demonstrates how pseudo‐transient continuation improves the efficiency and robustness of a Newton iteration within a non‐linear transient elasticity simulation. Pseudo‐transient continuation improves efficiency by enabling larger time steps than possible with a Newton iteration. Robustness improves because pseudo‐transient continuation recovers the convergence of Newton's method when the initial iterate is not within the region of local convergence. We illustrate the benefits of pseudo‐transient continuation on a non‐linear transient simulation of a buckling cylinder, including a comparison with a line search‐based Newton iteration. Copyright © 2008 John Wiley & Sons, Ltd.     

Stabilization of bi‐linear mixed finite elements for tetrahedra with enhanced interpolation using volume and area bubble functions

In order to overcome the oscillatory effects of the mixed bi‐linear Galerkin formulation for tetrahedral elements, a stabilization approach is presented. To this end the mixed method of incompatible modes and the mixed method of enhanced strains are reformulated, thus giving both the interpretation of a mixed finite element method with stabilization terms. For non‐linear problems, these are non‐linearly dependent on the current deformation state and therefore are replaced by linearly dependent stabilization terms. The approach becomes most attractive for the numerical implementation, since the use of quantities related to the previous Newton iteration step, typically arising for mixed‐enhanced elements, is completely avoided. The stabilization matrices for the mixed method of incompatible modes and the mixed method of enhanced strains are obtained with volume and area bubble functions. Various numerical examples are presented, which illustrate successfully the stabilization effect for bi‐linear tetrahedral elements. Copyright © 2007 John Wiley & Sons, Ltd.     

Computational techniques for optimization of problems involving non‐linear transient simulations

Many optimization problems in engineering require coupling a mathematical programming process to a numerical simulation. When the latter is non‐linear, the resulting computer time may become unaffordably large because three sequential procedures are nested: the outer loop is associated to the optimization process, the middle one corresponds to the time marching scheme and the innermost loop is required for solving iteratively the non‐linear system of equations at each time step. We propose four techniques for reducing CPU time. First, derive the initial values of state variables at each time (innermost loop) from those computed at the previous optimization iteration (outermost loop). Second, select time increment on the basis of those used for the previous optimization iteration. Third, define convergence criteria for the simulation problem on the basis of the optimization process, so that they are only as stringent as really needed. Finally, computations associated to the optimization are shown to be greatly reduced by adopting Newton–Raphson, or a variant, for solving the simulation problem. The effectiveness of these techniques is illustrated through application to three examples involving automatic calibration of non‐linear groundwater flow problems. The total number of iterations is reduced by a factor ranging between 1·7 and 4·6. Copyright © 1999 John Wiley & Sons, Ltd.     

Optimization‐based stability analysis of structures under unilateral constraints

This paper discusses an optimization‐based technique for determining the stability of a given equilibrium point of the unilaterally constrained structural system, which is subjected to the static load. We deal with the three problems in mechanics sharing the common mathematical properties: (i) structures containing no‐compression cables; (ii) frictionless contacts; and (iii) elastic–plastic trusses with non‐negative hardening. It is shown that the stability of a given equilibrium point of these structures can be determined by solving a maximization problem of a convex function over a convex set. On the basis of the difference of convex functions optimization, we propose an algorithm to solve the stability determination problem, at each iteration of which a second‐order cone programming problem is to be solved. The problems presented are solved for various structures to determine the stability of given equilibrium points. Copyright © 2008 John Wiley & Sons, Ltd.     

A rate‐consistent and intrinsically path‐dependent incremental method for plastic buckling problems

This paper presents an incremental predictor–corrector method which is able to handle the continuous spreading of elastic unloading and, therefore, is particularly well suited to solve plastic bucking problems. The method, which deals explicitly with rate variables and equations, is (i) rate consistent, because it leads to the ‘true’ tangent matrix, and (ii) intrinsically path‐dependent, because it enables an adequate identification and characterization of the onset of elastic unloading within each incremental step. Copyright © 2000 John Wiley & Sons, Ltd.     

Fatigue control of thin‐walled structures

The control of the ultimate fatigue behaviour of slender thin‐walled structures is treated in the paper. A macro‐ and micromechanical simulation model adopting the neural network approach is used for the numerical analysis of the problem. The numerical treatment of the non‐linear problems is made using the updated Lagrangian formulation of motion combined with the pseudo‐force technique in the FETM‐approach. Each step of the iteration approaches the solution of the linear problem and the feasibility of the parallel processing and neural network numerical techniques is established. The application on the actual slender bridge is made in order to demonstrate the efficiency of the procedures suggested. Copyright © 2003 John Wiley & Sons, Ltd.     

A novel numerical method of high‐order accuracy for flow in unsaturated porous media

We construct finite volume schemes of very high order of accuracy in space and time for solving the nonlinear Richards equation (RE). The general scheme is based on a three‐stage predictor–corrector procedure. First, a high‐order weighted essentially non‐oscillatory (WENO) reconstruction procedure is applied to the cell averages at the current time level to guarantee monotonicity in the presence of steep gradients. Second, the temporal evolution of the WENO reconstruction polynomials is computed in a predictor stage by using a global weak form of the governing equations. A global space–time DG FEM is used to obtain a scheme without the parabolic time‐step restriction caused by the presence of the diffusion term in the RE. The resulting nonlinear algebraic system is solved by a Newton–Krylov method, where the generalized minimal residual method algorithm of Saad and Schulz is used to solve the linear subsystems. Finally, as a third step, the cell averages of the finite volume method are updated using a one‐step scheme, on the basis of the solution calculated previously in the space–time predictor stage. Our scheme is validated against analytical, experimental, and other numerical reference solutions in four test cases. A numerical convergence study performed allows us to show that the proposed novel scheme is high order accurate in space and time. Copyright © 2011 John Wiley & Sons, Ltd.