A starting-point strategy for interior-point algorithms for shakedown analysis of engineering structures

Lower-bound shakedown analysis leads to nonlinear convex optimization problems with large numbers of unknowns and constraints, the solution of which can be obtained efficiently by interior-point algorithms. The performance of these algorithms strongly depends on the choice of the starting point. In general, starting points should be located inside the feasible region. In addition, they should also be well centred. Although there exist several heuristics for the construction of suitable starting points, these are restricted, as long as only the mathematical procedure is considered without taking into account the nature of the underlying mechanical problem. Thus, in this article, a strategy is proposed for choosing appropriate starting points for interior-point algorithms applied to shakedown analysis. This strategy is based on both the mathematical characteristics and the physical meaning of the variables involved. The efficiency of the new method is illustrated by numerical examples from the field of power plant engineering.     

A unified mathematical programming formulation of strain driven and interior point algorithms for shakedown and limit analysis

A mathematical programming formulation of strain‐driven path‐following strategies to perform shakedown and limit analysis for perfectly elastoplastic materials in an FEM context is presented. From the optimization point of view, standard arc‐length strain‐driven elastoplastic analyses, recently extended to shakedown, are identified as particular decomposition strategies used to solve a proximal point algorithm applied to the static shakedown theorem that is then solved by means of a convergent sequence of safe states. The mathematical programming approach allows: a direct comparison with other non‐linear programming methods, simpler convergence proofs and duality to be exploited. Owing to the unified approach in terms of total stresses, the strain‐driven algorithms become more effective and less non‐linear with respect to a self‐equilibrated stress formulation and easier to implement in the existing codes performing elastoplastic analysis. The elastic domain is represented avoiding any linearization of the yield function so improving both the accuracy and the performance. Better results are obtained using two different finite elements, one with a good behavior in the elastic range and the other suitable for performing elastoplastic analysis. The proposed formulation is compared with a specialized implementation of the primal–dual interior point method suitable to solve the problems at hand. Copyright © 2011 John Wiley & Sons, Ltd.     

Ellipsoidal load‐domain shakedown analysis with von Mises yield criterion: A robust optimization approach

The static formulation of elastic shakedown analysis, based upon Melan's lower‐bound theorem, can essentially be viewed as a robust optimization problem. This paper discusses an advantage that is enjoyed by taking this perspective. Specifically, we assume the von Mises yield criterion and an ellipsoidal load domain. In this setting, the shakedown analysis problem, which is viewed as robust second‐order cone programming, can be recast as semidefinite programming. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite element shakedown analysis of two‐dimensional structures

In the present paper, the formulation proposed by Casciaro and Garcea (Comput. Meth. Appl. Mech. Eng., 2002; 191 :5761–5792) and applied to the shakedown analysis of plane frames, is extended to the analysis of two‐dimensional flat structures in both the cases of plane‐stress and plane‐strain. The discrete formulation is obtained using a mixed finite element in which both stress and displacement fields are interpolated. The material is assumed to be elasto‐plastic and a linearization of the elastic domain is performed. The result is a versatile iterative scheme well suited to implementation in general purpose FEM codes. An extensive series of numerical tests is presented showing the reliability of the proposed formulation. Copyright © 2005 John Wiley & Sons, Ltd.     

A cycle-oriented incremental analysis of shakedown problems

An effective method for determination of a quasi-static shakedown loading and of a steady-state response is proposed. The classical optimization problem based on Melan's theorem is suitably reformulated to meet the requirements of incremental analysis. The main attention is focused on the reduction of a computational time required for a completion of a transient elastic–plastic phase of deformation. The method differs significantly from classical incremental analyses. Here, a load cycle is approximated by the finite number of loading systems covering the cycle. Each system is then combined with a separate domain of the structure in which the load system can be treated as a dominant one. In this manner, the structure consists of parts, each of them undergoing suitably chosen one-parameter loading only. Such a modification allows us to build a set of non-linear equations for all loading systems covering the whole load cycle. As a consequence the structure can be treated as the one in which the transient plastic phase of deformation is analysed load cycle by load cycle without making load increments inside the considered cycle. Due to this innovation a significant reduction of the computational time required for the solution of the steady-state response of the structure is obtained what is illustrated on 3-D frames. © 1997 John Wiley & Sons, Ltd.     

A non‐linear programming approach to kinematic shakedown analysis of composite materials

Using a Representative volume element (RVE) to represent the microstructure of periodic composite materials, this paper develops a non‐linear numerical technique to calculate the macroscopic shakedown domains of composites subjected to cyclic loads. The shakedown analysis is performed using homogenization theory and the displacement‐based finite element method. With the aid of homogenization theory, the classical kinematic shakedown theorem is generalized to incorporate the microstructure of composites. Using an associated flow rule, the plastic dissipation power for an ellipsoid yield criterion is expressed in terms of the kinematically admissible velocity. By means of non‐linear mathematical programming techniques, a finite element formulation of kinematic shakedown analysis is then developed leading to a non‐linear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the plastic dissipation power which is to be minimized and an upper bound to the shakedown load of a composite is then obtained. An effective, direct iterative algorithm is proposed to solve the non‐linear programming problem. The effectiveness and efficiency of the proposed numerical method have been validated by several numerical examples. This can serve as a useful numerical tool for developing engineering design methods involving composite materials. Copyright © 2005 John Wiley & Sons, Ltd.     

Kinematical shakedown analysis with temperature‐dependent yield stress

This paper describes a direct shakedown analysis of structures subjected to variable thermal and mechanical loading. The classical kinematical shakedown theorem is modified to be implemented with any displacement‐based finite elements. The plastic incompressibility condition is imposed by the penalty function method. The shakedown limit is found via a non‐linear mathematical programming procedure. Two numerical shakedown methods are developed and implemented to provide alternative numerical means. The temperature‐dependent material model is included in theoretical and numerical calculation in a simple way. Its effect on shakedown limit is investigated. The numerical examination for some pressure vessel structures subjected to thermal and mechanical loading shows a satisfying precision and efficiency of the methods presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Application of h‐adaptive FEM and Zarka's approach to analysis of shakedown problems

The Zarka shakedown approach and the h‐adaptive finite element method are applied to evaluate residual stresses resulting from arbitrary cyclic loading. Two error indicators are used to refine the mesh: the explicit residual one which controls accuracy of the momentum balance and the interpolation error indicator which controls approximation of the modified back stresses. Validation tests performed for the Zarka method of simplified shakedown analysis suggest that such an approach may be used to obtain a quick estimate of residual states with the error acceptable for engineering purposes. Thus, it has been applied to compute residual stresses arising from service load in railroad rails. Copyright © 2004 John Wiley & Sons, Ltd.     

Computation of limit and shakedown loads using a node‐based smoothed finite element method

This paper presents a novel numerical procedure for computing limit and shakedown loads of structures using a node‐based smoothed FEM in combination with a primal–dual algorithm. An associated primal–dual form based on the von Mises yield criterion is adopted. The primal‐dual algorithm together with a Newton‐like iteration are then used to solve this associated primal–dual form to determine simultaneously both approximate upper and quasi‐lower bounds of the plastic collapse limit and the shakedown limit. The present formulation uses only linear approximations and its implementation into finite element programs is quite simple. Several numerical examples are given to show the reliability, accuracy, and generality of the present formulation compared with other available methods. Copyright © 2011 John Wiley & Sons, Ltd.     

A generalised method for ratchet analysis of structures undergoing arbitrary thermo‐mechanical load histories

A novel approach is presented based upon the Linear Matching Method framework in order to directly calculate the ratchet limit of structures subjected to arbitrary thermo‐mechanical load histories. Traditionally, ratchet analysis methods have been based upon the fundamental premise of decomposing the cyclic load history into cyclic and constant components, respectively, in order to assess the magnitude of additional constant loading a structure may accommodate before ratcheting occurs. The method proposed in this paper, for the first time, accurately and efficiently calculates the ratchet limit with respect to a proportional variation between the cyclic primary and secondary loads, as opposed to an additional primary load only. The method is a strain‐based approach and utilises a novel convergence scheme in order to calculate an approximate ratchet boundary based upon a predefined target magnitude of ratchet strain per cycle. The ratcheting failure mechanism evaluated by the method leads to less conservative ratchet boundaries compared with the traditional Bree solution. The method yields the total and plastic strain ranges as well as the ratchet strains for various levels of loading between the ratchet and limit load boundaries. Two example problems have been utilised in order to verify the proposed methodology. Copyright © 2015 John Wiley & Sons, Ltd.     

Combined interior‐point method and semismooth Newton method for frictionless contact problems

In the present paper, a solution scheme is proposed for frictionless contact problems of linear elastic bodies, which are discretized using the finite element method with lower order elements. An approach combining the interior‐point method and the semismooth Newton method is proposed. In this method, an initial active set for the semismooth Newton method is obtained from the approximate optimal solution by the interior‐point method. The simplest node‐to‐node contact model is considered in the present paper, that is, pairs of matching nodes exist on the contact surfaces. However, the discussions can be easily extended to a node‐to‐segment or segment‐to‐segment contact model. In order to evaluate the proposed method, a number of illustrative examples of the frictionless contact problem are shown. The proposed combined method is compared with the interior‐point method and the semismooth Newton method. Two numerical examples that are difficult to solve using the semismooth Newton method are solved effectively using the proposed combined method. It is shown that the proposed method converges within far fewer iterations than the semismooth Newton methods or the interior‐point method. Copyright © 2009 John Wiley & Sons, Ltd.     

A genetic algorithm for sequencing type problems in engineering design

A genetic algorithm for engineering applications that involve sequencing of operations is proposed and demonstrated. Such applications are known as travelling salesman problems in operations research literature. The proposed algorithm uses some new operators that are different from those typically used in genetic algorithms. Some enhancements for improving performance of the algorithm are also described. Treatment of two salesmen in the problem is also discussed. Results for test problems, including a vehicle A-pillar subassembly welding sequence application, show performance of the proposed algorithm to be quite robust. © 1997 John Wiley & Sons, Ltd.     

A Koiter reduction technique for the nonlinear thermoelastic analysis of shell structures prone to buckling

The multi-modal Koiter method is a reduction technique for estimating quickly the nonlinear buckling response of structures under mechanical loads requiring a fine discretization. The reduced model is based on a quadratic approximation of the full model using a few linear buckling modes and their second order corrections, followed by the projection of the equilibrium equations onto the modal subspace. In this article, the method is reformulated for geometrically nonlinear thermoelastic analyses of shell structures. The starting point is an isogeometric solid-shell discretization with an accurate modeling of thermal strains, temperature-dependent materials and general temperature distributions. The equilibrium path is defined as a displacement versus temperature amplifier curve, while mechanical loads are kept constant. The strain energy nonlinearity with respect to the temperature amplifier is the main obstacle in the definition of an accurate reduced model. This task is achieved by coherent asymptotic expansions in mixed form using independent stress variables at the integration points and accurate linear buckling modes obtained by a two-point mixed eigenvalue problem. Structures made of isotropic and multi-layered composite materials, including variable angle tow laminates, are considered to show the accuracy of the proposed method.     

Large‐deformation and friction analysis of non‐linear elastic cable networks by second‐order cone programming

A new formulation is presented for equilibrium shape analysis of cable networks considering geometrical and material non‐linearities. Friction between cables and joint devices is also considered. The second‐order cone programming (SOCP) problem which has the same solution as that of minimization of total potential energy is solved to obtain the equilibrium configuration. The optimality conditions are derived to verify that the solution satisfies equilibrium conditions and friction laws. Since no assumption on stress state is needed in the proposed method, no process of trial and error is involved. No effort is required to develop any analysis software because SOCP can be solved by using the well‐developed software based on the interior‐point method. Copyright © 2002 John Wiley & Sons, Ltd.     

Optimization of large‐scale truss structures using sparse SAND formulations

Sparsity features of simultaneous analysis and design (SAND) formulations are studied and exploited for optimization of large‐scale truss structures. Three formulations are described and implemented with an existing analysis code. SAND formulations have large number of variables; however, gradients of the functions and Hessian of the Lagrangian are quite sparsely populated. Therefore, this structure of the problem is exploited and an optimization algorithm with sparse matrix capability is used to solve large‐scale problems. An existing analysis software is integrated with an optimizer based on a sparse sequential quadratic programming (SQP) algorithm to solve sample problems. The formulations and algorithms work quite well for the sample problems, and their performances are compared and discussed. For all the cases considered, the SAND formulations out perform the conventional formulation except one case. Further research is suggested to fully study and utilize sparse features of the alternative SAND formulations and to develop more efficient sparse solvers for large‐scale and more complex applications. Copyright © 2006 John Wiley & Sons, Ltd.