Reduced augmented Lagrangian bi-conjugate gradient method for impact-contact problems

To avoid the numerical oscillation of the penalty method and non-compatibility with explicit operators of conventional Lagrange multiplier methods used in transient contact problems to enforce surface contact conditions, a new approach to enforcing surface contact constraints for the transient nonlinear finite element problems, referred to as “the reduced augmented Lagrangian bi-conjugate gradient method (ALCG)”, is developed in this paper. Based on the nonlinear constrained optimization theory and is compatible with the explicit time integration scheme, this approach can also be used in implicit scheme naturally. The new surface contact constraint method presented has significant advantages over the widely adopted penalty function methods and the conventional Lagrangian multiplier methods. The surface contact constraints are satisfied more accurately for each step by the algorithm, so the oscillation of numerical solution for the explicit scheme is depressed. Through the development of new iteration strategy for solving nonlinear equations, ALCG method improves the computational efficiency greatly. Project supported by State Education Commission Doctoral Foundation and Natural Science Foundation of Liaoning Province.     

A boundary element formulation for wear modeling on 3D contact and rolling-contact problems

The present work shows a new numerical treatment for wear simulation on 3D contact and rolling-contact problems. This formulation is based on the boundary element method (BEM) for computing the elastic influence coefficients and on projection functions over the augmented Lagrangian for contact restrictions fulfillment. The constitutive equations of the potential contact zone are Signorini’s contact conditions, Coulomb’s law of friction and Holm–Archard’s law of wear. The proposed methodology is applied to predict wear on different contact and rolling-contact problems. Results are validated with numerical solutions and semi-analytical models presented in the literature. The BEM considers only the degrees of freedom involved on these kind of problems (those on the solids surfaces), reducing the number of unknowns and obtaining a very good approximation on contact tractions using a low number of elements. Together with the formulation, an acceleration strategy is presented allowing to reduce the times of resolution.     

3D frictional contact of anisotropic solids using BEM

A numerical method to study three-dimensional (3D) contact problems in solids with anisotropic elastic behavior is developed in this work. This formulation is based on the Boundary Element Method (BEM) for computing the elastic influence coefficients and on projection functions over the augmented Lagrangian for contact restrictions fulfillment. The constitutive equations of the potential contact zone are Signorini’s contact conditions and Coulomb’s law of friction. The formulation uses a recently introduced explicit approach for fundamental solutions evaluation, which are valid for general anisotropic behavior meanwhile mathematical degeneracies are allowed. The accuracy and robustness of the proposed method is illustrated by solving some examples previously presented in the literature. This approach is further applied to study the influence of solids anisotropy on the contact problem.     

Finite deformations of Ogden's materials under impact loading

The present paper is devoted to the modeling of finite deformations of a hyperelastic body described by Ogden's model under contact/impact conditions. Frictional contact problems are solved by means of the bi-potential method. The first order algorithm is applied to integrate the equation of motion. The total Lagrangian formulation is adopted to describe the geometrically non-linear behavior. For the finite element implementation, the explicit expression of the tangent operator is derived including the case of repeated eigenvalues. A numerical example is given to illustrate efficiency and accuracy of the method.     

Contact stresses: a short survey of models and methods of computations

Contact stresses are identified as normal and tangential forces between contacting solids. The normal stresses are modeled using unilateral and complementary conditions, elastic response and normal compliance. Friction laws describe the tangential traction. Friction of materials depends on pressure, sliding velocity, surface temperature, time of contact, surface roughness and presence of wear debris. Phenomenological, micro-mechanical and atomic-scale models as well as non-classical models of anisotropic and heterogeneous friction are important steps in the development of friction modeling. Sophisticated friction models are desirable in vibrating systems, materials processing, rolling contacts, rubber and polymers, geomechanics, bioengineering and living systems. Main numerical methods in contact mechanics are: finite element method, boundary element method and discrete element method. To include specific contact constraints, the following computing techniques are applied: Lagrange multipliers, penalty function, perturbated and augmented Lagrangian methods, mathematical programming methods. The advances of adhesion and impact modeling are outlined in this paper.     

A frictional mortar contact approach for the analysis of large inelastic deformation problems

A multibody frictional mortar contact formulation (Gitterle et al., 2010) is extended for the simulation of solids undergoing finite strains with inelastic material behavior. The framework includes the modeling of finite strain inelastic deformation, the numerical treatment of frictional contact conditions and specific finite element technology. Several well-established and recent models are employed for each of these building blocks to capture the distinct physical aspects of the deformation behavior. The approach is based on a mortar formulation and the enforcement of contact constraints is realized with dual Lagrange multipliers. The introduction of nonlinear complementarity functions into the frictional contact conditions combined with the global equilibrium leads to a system of nonlinear equations, which is solved in terms of the semi-smooth Newton method. The resulting method can be interpreted as a primal–dual active set strategy (PDASS) which deals with contact nonlinearities, material and geometrical nonlinearities in one iterative scheme. The consistent linearization of all building blocks of the framework yields a robust and highly efficient approach for the analysis of metal forming problems. The effect of finite inelastic strains on the solution behavior of the PDASS method is examined in detail based on the complementarity parameters. A comprehensive set of numerical examples is presented to demonstrate the accuracy and efficiency of the approach against the traditional node-to-segment penalty contact formulation.     

A finite element finite-strain formulation for modeling colliding blocks of Gent materials

The present paper is devoted to the analysis of the contact/impact problems with Coulomb friction and large deformation between two hyperelastic bodies of Gent model. The total Lagrangian formulation is adopted to describe the geometrically non-linear behavior. For the finite element implementation, the explicit expression of the incremental law of Gent model is derived. A first order algorithm is applied for the numerical integration of the time-discretized equation of motion. Efficiency and accuracy of the resulting method is illustrated on a two-dimensional static contact problem and a three-dimensional dynamic contact problem as compared with ANSYS simulations.     

An implicit method for frictional contact, impact and rolling

In this paper an implicit method for frictional contact, impact and rolling is suggested. A nonclassical formulation of a two-dimensional hyperelastic body unilaterally constrained to rigid supports is proposed by following the ideas of Moreau and Jean. A total Lagrangian formulation of the system is given. The elastic properties are defined by coupling the second Piola–Kirchhoff stress to the Green–Lagrange strain via the Kirchhoff–St. Venant law. The equation of motion is written in the spirit of Moreau by using the mean value impulses introduced by Jean. The mean value impulses appear explicitly in the equation of motion. In such manner the treatment of nonconstant kinematic transformation matrices becomes straightforward. The rigid supports are described by smooth functions. By utilizing these functions and the mean value impulses, new contact/impact laws of Signorini and Coulomb type are formulated. The governing equations are solved by a nonsmooth Newton method. This is performed by following the augmented Lagrangian approach and deriving the consistent stiffness matrix as well as the contact stiffness matrices. Three two-dimensional examples are solved by the method: a contact problem, an impact problem and a rolling contact problem.     

Contact-impact formulation for multi-body systems using component mode synthesis

The efficiency and accuracy are two most concerned issues in the modeling and simulation of multi-body systems involving contact and impact. This paper proposed a formulation based on the component mode synthesis method for planar contact problems of flexible multi-body systems. A flexible body is divided into two parts： a contact zone and an un-contact zone. For the un-contact zone, by using the fixed-interface substructure method as reference, a few low-order modal coordinates are used to replace the nodal coordinates of the nodes, and meanwhile, the nodal coordinates of the local impact region are kept unchanged, therefore the total degrees of freedom （DOFs） are greatly cut down and the computational cost of the simulation is significantly reduced. By using additional constraint method, the impact constraint equations and kinematic constraint equations are derived, and the Lagrange equations of the first kind of flexible multi-body system are obtained. The impact of an elastic beam with a fixed half disk is simulated to verify the efficiency and accuracy of this method.     

A method for solving the factorized vorticity-stream function equations by finite elements

A new finite element method for solving the time-dependent incompressible Navier-Stokes equations with general boundary conditions is presented. The two second-order partial differential equations for the vorticity and the stream function are factorized, apart from the non-linear advection term, by eliminating the coupling due to the double specification on the stream function at (a part of) the boundary. This is achieved by reducing the no-slip boundary conditions to projection integral conditions for the vorticity field and by evaluating the relevant quantities involved according to an extension of the method of Glowinski and Pironneau for the biharmonic problem. Time integration schemes and iterative algorithms are introduced which require the solution only of banded linear systems of symmetric type. The proposed finite element formulation is compared with its finite difference equivalent by means of a few numerical examples. The results obtained using 4-noded bilinear elements provide an illustration of the superiority of the finite element based spatial discretization.     

A Newton method for three-dimensional fretting problems

The present paper concerns the numerical treatment of fretting problems using a finite element analysis. The governing equations resulting from a formal finite element discretization of an elastic body with a potential contact surface are considered in a quasi-static setting. The constitutive equations of the potential contact surface are Signorinis contact conditions, Coulombs law of friction and Archards law of wear. Using a backward Euler time discretization and an approach based on projections, the governing equations are written as an augmented Lagrangian formulation which is implemented and solved using a Newton algorithm for three-dimensional fretting problems of didactic nature. Details concerning the implementation are provided.     

On the distributed Lagrange multiplier/fictitious domain method for rigid‐particle‐laden flows: a proposition for an alternative formulation of the Lagrange multipliers

The distributed Lagrange multiplier/fictitious domain method proposed for the direct numerical simulation of particle‐laden flows is considered in this work. First, it is demonstrated that improved accuracy is obtained with a coupled numerical scheme, whereby the pressure and the Lagrange multiplier fields enforcing incompressibility and rigid body motion, respectively, are calculated and applied together. However, the convergence characteristics of the iterative solution of the coupled scheme are poor because symmetric but indefinite and poorly conditioned matrices are produced. An analysis is then presented, which suggests that the cause for the matrix pathologies lies in the interaction of the respective matrix operators enforcing incompressibility and rigid body motion. On the basis of this analysis, an alternative formulation is developed for the Lagrange multipliers, being now composed of a set of forces distributed only on the particle boundary together with a set of couples distributed within the particle core. The new formulation is tested with several types of flows with stationary or moving particles under creeping or finite Reynolds number conditions and it is demonstrated that it produces correct results and better conditioned matrices, thus enabling faster and more reliable convergence of the conjugate gradient method. The analysis and tests, therefore, support the expectation that the proposed formulation is promising and worthy of further study and improvement. Copyright © 2011 John Wiley & Sons, Ltd.     

Local and Global Injective Solution of the Rotationally Symmetric Sphere Problem

There are problems in the classical linear theory of elasticity whose closed form solutions, while satisfying the governing equations of equilibrium together with well-posed boundary conditions, predict the existence of regions, often quite small, inside the body where material overlaps. Of course, material overlapping is not physically realistic, and one possible way to prevent it combines linear theory with the requirement that the deformation field be injective. A formulation of minimization problems in classical linear elasticity proposed by Fosdick and Royer [3] imposes this requirement through a Lagrange multiplier technique. An existence theorem for minimizers of plane problems is also presented. In general, however, it is not certain that such minimizers exist. Here, the Euler–Lagrange equations corresponding to a family of three-dimensional problems is investigated. In classical linear elasticity, these problems do not have bounded solutions inside a body of anisotropic material for a range of material parameters. For another range of parameters, bounded solutions do exist but yield stresses that are infinite at a point inside the body. In addition, these solutions are not injective in a region surrounding this point, yielding unrealistic behavior such as overlapping of material. Applying the formulation of Fosdick and Royer on this family of problems, it is shown that both the displacements and the constitutive part of the stresses are bounded for all values of the material parameters and that the injectivity constraint is preserved. In addition, a penalty functional formulation of the constrained elastic problems is proposed, which allows to devise a numerical approach to compute the solutions of these problems. The approach consists of finding the displacement field that minimizes an augmented potential energy functional. This augmented functional is composed of the potential energy of linear elasticity theory and of a penalty functional divided by a penalty parameter. A sequence of solutions is then constructed, parameterized by the penalty parameter, that converges to a function that satisfies the first variation conditions for a minimizer of the constrained minimization problem when this parameter tends to infinity. This approach has the advantages of being mathematically appealling and computationally simple to implement.     

平面柔性多体系统正碰撞动力学建模理论研究

针对目前柔性多体系统碰撞动力学建模方法存在的不足,对影响碰撞动力学仿真的主要因素如柔性体建模和碰撞初始条件进行分析,建立起基于变约束的柔性体碰撞动力学方程。首先,为了解决子结构法在处理碰撞界面搜索时面临的难题,引入多体系统柔性体有限元描述方法,推导出凸形柔性体接触点间法向位移约束的二阶导数形式。其次,从碰撞引起的接触界面速度不连续机理出发,结合连续介质力学间断面理论,给出碰撞瞬时由物体本身物理性质决定的接触位置处速度跳跃公式。最后对两弹性圆盘低速碰撞问题进行数值仿真。结果表明本文提出的改进方法符合力学基本原理,仿真结果满足收敛性要求。     

A 3D FEM–BEM rolling contact formulation for unstructured meshes

This work presents a new formulation for solving 3D steady-state rolling contact problems. The convective terms for computing the tangential slip velocities involved in the rolling problem, are evaluated using a new approximation inspired in numerical fluid dynamics techniques for unstructured meshes. Moreover, the elastic influence coefficients of the surface points in contact are approached by means of the finite element method (FEM) and/or the boundary element method (BEM). The contact problem is based on an Augmented Lagrangian Formulation and the use of projection functions to establish the contact restrictions. Finally, the resulting nonlinear equations set is solved using the generalized Newton method with line search (GNMls), presenting some acceleration strategies as: a new and more simplified projection operator, which makes it possible to obtain a quasi-complementarity of the contact variables, reducing the number of contact problem unknowns, and using iterative solvers. The presented methodology is validated solving some rolling contact problems and analyzed for some unstructured mesh examples.     

共轭梯度法求解非线性多宗量稳态传热反问题

应用共轭梯度法求解非线性多宗量稳态热传导反问题。采用八节点的等参单元在空间上进行离散,建立了便于敏度分析的非线性正演和反演的有限元模型,可直接求导进行敏度分析。给出了相关的数值验证,对测量误差及测点数目的影响作了初步探讨,结果表明,采用的算法能够对非线性稳态热传导中导热系数和边界条件联合反问题进行有效的求解,并具有较高精度。     

Nonlinear vibration analysis of micro-plates based on strain gradient elasticity theory

In this study, non-linear free vibration of micro-plates based on strain gradient elasticity theory is investigated. A general form of Mindlin’s first-strain gradient elasticity theory is employed to obtain a general Kirchhoff micro-plate formulation. The von Karman strain tensor is used to capture the geometric non-linearity. The governing equations of motion and boundary conditions are obtained in a variational framework. The Homotopy analysis method is employed to obtain an accurate analytical expression for the non-linear natural frequency of vibration. For some specific values of the gradient-based material parameters, the general plate formulation can be reduced to those based on some special forms of strain gradient elasticity theory. Accordingly, three different micro-plate formulations are introduced, which are based on three special strain gradient elasticity theories. It is found that both geometric non-linearity and size effect increase the natural frequency of vibration. In a micro-plate having a thickness comparable with the material length scale parameter, the strain gradient effect on increasing the non-linear natural frequency is higher than that of the geometric non-linearity. By increasing the plate thickness, the strain gradient effect decreases or even diminishes. In this case, geometric non-linearity plays the main role on increasing the natural frequency of vibration. In addition, it is shown that for micro-plates with some specific thickness to length scale parameter ratios, both geometric non-linearity and size effect have significant role on increasing the frequency of non-linear vibration.     

Finite element implementation of a generalised non-local rate-dependent crystallographic formulation for finite strains

This work describes the finite element implementation of a generalised strain gradient and rate-dependent crystallographic formulation for finite strains and general anisothermal conditions based on a multiplicative decomposition of the deformation gradient. The implementation involved the development of both a novel finite element formulation to determine the spatial slip rate gradients at each material point, and an implicit numerical integration scheme at the constitutive level to update the stresses and solution dependent variables. The time-integration procedure uses a Newton–Raphson scheme with a single level of iteration to solve the incremental non-linear equations associated with the non-local constitutive formulation. Closed-form solutions for the relevant fourth-order Jacobian tensors are given. The proposed numerical scheme is formulated in a general form and hence should be applicable to most existing crystallographic models. The crystallographic formulation is then used to investigate the effect of the morphology and volume fraction of the reinforcing phase of a two-phase single crystal on its macroscopic behaviour.