Deflation for semismooth equations

Variational inequalities can in general support distinct solutions. In this paper we study an algorithm for computing distinct solutions of a variational inequality, without varying the initial guess supplied to the solver. The central idea is the combination of a semismooth Newton method with a deflation operator that eliminates known solutions from consideration. Given one root of a semismooth residual, deflation constructs a new problem for which a semismooth Newton method will not converge to the known root, even from the same initial guess. This enables the discovery of other roots. We prove the effectiveness of the deflation technique under the same assumptions that guarantee locally superlinear convergence of a semismooth Newton method. We demonstrate its utility on various finite- and infinite-dimensional examples drawn from constrained optimization, game theory, economics and solid mechanics.Comment: 24 pages, 3 figure

Smooth finite strain plasticity with non-local pressure support

The aim of this work is to introduce an alternative framework to solve problems of finite strain elastoplasticity including anisotropy and kinematic hardening coupled with any isotropic hyperelastic law. After deriving the constitutive equations and inequalities without any of the customary simplifications, we arrive at a new general elasto-plastic system. We integrate the elasto-plastic algebraico-differential system and replace the loading–unloading condition by a Chen–Mangasarian smooth function to obtain a non-linear system solved by a trust region method. Despite being non-standard, this approach is advantageous, since quadratic convergence is always obtained by the non-linear solver and very large steps can be used with negligible effect in the results. Discretized equilibrium is, in contrast with traditional approaches, smooth and well behaved. In addition, since no return mapping algorithm is used, there is no need to use a predictor. The work follows our previous studies of element technology and highly non-linear visco-elasticity. From a general framework, with exact linearization, systematic particularization is made to prototype constitutive models shown as examples. Our element with non-local pressure support is used. Examples illustrating the generality of the method are presented with excellent results

A new semi-implicit formulation for multiple-surface ow rules in multiplicative plasticity

We propose new integration scheme

Continuation method for nonlinear complementarity problems via normal maps

Cataloged from PDF version of article.In a recent paper by Chen and Mangasarian (C. Chen, O.L. Mangasarian, A class of smoothing functions for nonlinear and mixed complementarity problems, Computational Optimization and Applications 2 (1996), 97±138) a class of parametric smoothing functions has been proposed to approximate the plus function present in many optimization and complementarity related problems. This paper uses these smoothing functions to approximate the normal map formulation of nonlinear complementarity problems (NCP). Properties of the smoothing function are investigated based on the density functions that de®nes the smooth approximations. A continuation method is then proposed to solve the NCPs arising from the approximations. Su cient conditions are provided to guarantee the boundedness of the solution trajectory. Furthermore, the structure of the subproblems arising in the proposed continuation method is analyzed for di erent choices of smoothing functions. Computational results of the continuation method are reported. Ó 1999 Elsevier Science B.V. All rights reserved

Differential-Algebraic Equations and Beyond: From Smooth to Nonsmooth Constrained Dynamical Systems

The present article presents a summarizing view at differential-algebraic equations (DAEs) and analyzes how new application fields and corresponding mathematical models lead to innovations both in theory and in numerical analysis for this problem class. Recent numerical methods for nonsmooth dynamical systems subject to unilateral contact and friction illustrate the topicality of this development.Comment: Preprint of Book Chapte

An MPCC Formulation and Its Smooth Solution Algorithm for Continuous Network Design Problem

Continuous network design problem (CNDP) is searching for a transportation network configuration to minimize the sum of the total system travel time and the investment cost of link capacity expansions by considering that the travellers follow a traditional Wardrop user equilibrium (UE) to choose their routes. In this paper, the CNDP model can be formulated as mathematical programs with complementarity constraints (MPCC) by describing UE as a non-linear complementarity problem (NCP). To address the difficulty resulting from complementarity constraints in MPCC, they are substituted by the Fischer-Burmeister (FB) function, which can be smoothed by the introduction of the smoothing parameter. Therefore, the MPCC can be transformed into a well-behaved non-linear program (NLP) by replacing the complementarity constraints with a smooth equation. Consequently, the solver such as LINDOGLOBAL in GAMS can be used to solve the smooth approximate NLP to obtain the solution to MPCC for modelling CNDP. The numerical experiments on the example from the literature demonstrate that the proposed algorithm is feasible.</p