https://hal.archives-ouvertes.fr/hal-01513346/document

Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement

Large linear systems of saddle point type arise in a wide variety of applications throughout computational science and engineering. Due to their indefiniteness and often poor spectral properties, such linear systems represent a significant challenge for solver developers. In recent years there has been a surge of interest in saddle point problems, and numerous solution techniques have been proposed for this type of system. The aim of this paper is to present and discuss a large selection of solution methods for linear systems in saddle point form, with an emphasis on iterative methods for large and sparse problems.

/pdf/numerical-solution-of-saddle-point-problems-l583iob5c2.pdf

Numerical solution of saddle point problems

Introduction to linear and nonlinear programming

IN two previous notes1, Prof. Max Born and I have shown that one can obtain a theory of superconductivity by taking account of the fact that the interaction of the electrons with the ionic lattice is appreciable only near the boundaries of Brillouin zones, and particularly strong near the corners of these. This leads to the criterion that the metal should be superconductive if a set of corners of a Brillouin zone is lying very near the Fermi surface, considered as a sphere, which limits the region in the momentum space completely filled with electrons.

On the theory of superconductivity

This paper describes modern techniques used to solve contact problems within Computational Mechanics. On the basis of a continuum description of contact, the mathematical structure of the contact formulation for finite element methods is derived. Emphasis is also placed on the constitutive behavior at the contact interface for normal and tangential (frictional) contact. Furthermore, different discretization schemes currently applied to solve engineering problems are formulated for small and finite strain problems. These include isoparametric interpolations, node-to-segment discretizations and also mortar and Nitsche techniques. Furthermore, several algorithms related to spatial contact search and fulfillment of the inequality constraints at the contact interface are discussed. Here, especially the penalty and Lagrange multiplier schemes are considered and also SQP- and linear-programming methods are reviewed.


Keywords:

contact mechanics;
friction;
penalty method;
Lagrange multiplier method;
contact algorithms;
finite element method;
finite deformations;
discretization methods

Computational Contact Mechanics

This paper is focussed on path following methods which are derived from consistent linearizations. The linearization procedure leads to some well-known constraint equations—like the constant arc length in the load-displacement space—and to different formulations than those given in the literature. A full Newton scheme for the unknown quantities (displacements and load parameter) can be formulated. A comparison of the derived algorithms with other path following methods is included to show advantages and limits of the methods. Using the linearization technique together with scaling a family of path following methods is introduced. Here, the scaling bypasses physical inconsistencies associated with mixed quantities like displacements and rotations in the global vector of the unknowns. Several possible scaling procedures are derived from a unified formulation. A discussion of these methods by means of numerical examples shows that up to now the choice of the scaling procedure is problem-dependent. If the arc-length methods are combined with a modified Newton method, an enhancement of the algorithms is achieved by line search techniques. Here, a simple but efficient line search was implemented and compared with a numerical relaxation technique. Both methods improve the convergence rate considerably.

Consistent linearization for path following methods in nonlinear FE analysis

A new locking-free brick element technique for large deformation problems in elasticity ☆

The practical behaviour of problems exhibiting bifurcation with secondary branches cannot be studied in general by using standard path‐following methods such as arc‐length schemes. Special algorithms have to be employed for the detection of bifurcation and limit points and furthermore for branch‐switching. Simple methods for this purpose are given by inspection of the determinant of the tangent stiffness matrix or the calculation of the current stiffness parameter. Near stability points, the associated eigenvalue problem has to be solved in order to calculate the number of existing branches. The associated eigenvectors are used for a perturbation of the solution at bifurcation points. This perturbation is performed by adding the scaled eigenvector to the deformed configuration in an appropriate way. Several examples of beam and shell problems show the performance of the method.

A simple method for the calculation of postcritical branches

http://yoksis.bilkent.edu.tr/pdf/files/10.1016-j.cma.2011.10.014.pdf

Three-dimensional mortar-based frictional contact treatment in isogeometric analysis with NURBS

New finite elements are needed as well in research as in industry environments for the development of virtual prediction techniques. The design and implementation of novel finite elements for specific purposes is a tedious and time consuming task, especially for nonlinear formulations. The automation of this process can help to speed up this process considerably since the generation of the final computer code can be accelerated by order of several magnitudes. This book provides the reader with the required knowledge needed to employ modern automatic tools like AceGen within solid mechanics in a successful way. It covers the range from the theoretical background, algorithmic treatments to many different applications. The book is written for advanced students in the engineering field and for researchers in educational and industrial environments

Peter Wriggers

Papers

Consistent linearization for path following methods in nonlinear FE analysis

A new locking-free brick element technique for large deformation problems in elasticity ☆

A simple method for the calculation of postcritical branches

Three-dimensional mortar-based frictional contact treatment in isogeometric analysis with NURBS

Automation of Finite Element Methods