A Domain Decomposition Method for Frictional Contact Problem

The bilateral or unilateral contact problem with Coulomb friction between two elastic bodies is considered [1]. An algorithm is introduced to solve the resulting finite element system by a non‐overlapping domain decomposition method. The global problem is transformed to a smaller problem on the contact surface. The solution is obtained by using a successive approximation method, in each step of this algorithm we solve two intermediate problems the first with prescribed tangential pressure and the second with prescribed normal pressure.     

An H‐matrix Type Preconditioner For Frictional Contact Problems

The bilateral or unilateral contact problem with Coulomb friction between two elastic bodies is considered [1]. An algorithm is introduced to solve the resulting finite element system by a non‐overlapping domain decomposition method [2, 3]. The global problem is transformed to a independant local problems posed in each bodie and a problem posed on the contact surface (the interface problem). The solution is obtained by using a successive approximation method, in each step of this algorithm we solve two intermediate problems the first with prescribed tangential pressure and the second with prescribed normal pressure [8]. Our preconditioner construction is based on the application of the H‐matrix technique [6, 7] together with the representation of the H^1/2 seminorm by a sum of partial seminorms [4]. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

求解摩擦接触问题的一个非内点光滑化算法

给出了一个求解三维弹性有摩擦接触问题的新算法,即基于NCP函数的非内点光滑化算法．首先通过参变量变分原理和参数二次规划法,将三维弹性有摩擦接触问题的分析归结为线性互补问题的求解;然后利用NCP函数,将互补问题的求解转换为非光滑方程组的求解;再用凝聚函数对其进行光滑化,最后用NEWTON法解所得到的光滑非线性方程组．方法具有易于理解及实现方便等特点．通过线性互补问题的数值算例及接触问题实例证实了该算法的可靠性与有效性．     

Mixed L 2-Wasserstein Optimal Mapping Between Prescribed Density Functions

A time-dependent minimization problem for the computation of a mixed L ²-Wasserstein distance between two prescribed density functions is introduced in the spirit of Ref. 1 for the classical Wasserstein distance. The optimum of the cost function corresponds to an optimal mapping between prescribed initial and final densities. We enforce the final density conditions through a penalization term added to our cost function. A conjugate gradient method is used to solve this relaxed problem. We obtain an algorithm which computes an interpolated L ²-Wasserstein distance between two densities and the corresponding optimal mapping.     

A domain decomposition method for two-body contact problems with Tresca friction

The paper analyzes a continuous and discrete version of the Neumann-Neumann domain decomposition algorithm for two-body contact problems with Tresca friction. Each iterative step consists of a linear elasticity problem for one body with displacements prescribed on a contact part of the boundary and a contact problem with Tresca friction for the second body. To ensure continuity of contact stresses, two auxiliary Neumann problems in each domain are solved. Numerical experiments illustrate the performace of the proposed approach.     

Contact interaction of composite shells,subjected to follower loads,with a rigid convex foundation

Based on a 7-parameter shell model, a numerical algorithm is developed for solving the contact problem for a multilayered composite shell lying on a rigid convex foundation, which is subjected to a follower pressure and undergoes arbitrarily large rotations. A new geometrically exact solid shell element is formulated, which permits one to solve the nonlinear deformation problem for thin-walled composite structures under unilateral contact constraints by using a small number of load steps. The calculation of a homogeneous ring and an angle-ply toroidal shell interacting with plane and cylindrical foundations is considered.     

Pressure of a plane stamp of nearly circular cross section on an elastic half-space : PMM vol. 40, n 6, 1976, pp. 1143–1145

An approximate method of solving the contact problem of impressing a plane stamp of nearly circular cross section into an elastic half-space is suggested. The friction of the contact surface is neglected. A numerical algorithm for the method is produced. An elliptical and rectangular stamps are considered as examples.There is no general method of solving the problems for stamps of nearly circular cross section. Apart from the classical problem of a plane elliptical stamp, the literature gives solutions for the problems of polygonal stamps, with each problem however requiring a different approach. An approximate solution for the problem of impressing a stamp of nearly circular cross section into an elastic half-space is given in [1]. The method makes it possible to use the same approach to solve the contact problem for an arbitrary region of contact, and to construct an universal numerical algorithm. The program can be adapted to each particular case by making the corresponding changes in the procedure of computing the Fourier coefficients of the equation of the boundary of the area of contact. Below a numerical algorithm for the approximate method in question is given. A more effective formulation of the solution is given for the case of the elliptical stamp.     

Multibody elastic contact analysis by quadratic programming

A quadratic programming method for contact problems is extended to a general problem involving contact ofn elastic bodies. Sharp results of quadratic programming theory provide an equivalence between the originaln-body contact problem and the simplex algorithm used to solve the quadratic programming problem. Two multibody examples are solved to illustrate the technique.     

Contact interaction of cylindrical shells of different lengths

We solve a problem of determination of the contact rigidity of a structural joint of a steel thin-walled pipe under the action of internal pressure with a shroud made of a composite material using a special technology. A mathematical model of the contact interaction of the pipe and shroud modeled by cylindrical shells of different length was constructed using the classical Kirchhoff–Love theory of shells. We obtain an analytic solution of the contact problem under conditions of ideal contact of the elements of the structural joint by the conjugation method. A numerical analysis of the influence of geometric and physicomechanical characteristics of the shroud on the contact pressure and rigidity of the shrouded pipe is presented.     

On the contact interaction between a strip-shaped punch and a linearly-deformable base through a covering of variable thickness

The contact problem of the frictionless penetration of a punch with strip-shaped section into the surface of a linearly-deformable base protected by a thin elastic layer (covering) of variable thickness, the stiffness of which is comparable to or smaller than that of the supporting elastic body, is investigated. A Fredholm integral equation of the second kind is obtained for the unknown contact pressure with a coefficient in front of the leading term that is a fairly arbitrary function of the longitudinal coordinate. To solve it the Bubnov-Galerkin projection method is used in which the coordinate elements are chosen to be a system of orthogonal polynomials and delta-shaped functions [1, 2] (variational-difference method), together with an algorithm for the required asymptotic expansions [3] when the above-mentioned coefficient is small. In the special case of an elastic half-space protected by a covering of constant thickness, the results obtained are compared with the corresponding characteristics given in [4].     

Shape and topology optimization of contact problems by level set methods

This paper deals with the numerical solution of a topology and shape optimization problems of an elastic body in unilateral contact with a rigid foundation. The contact problem with the prescribed friction is considered. The structural optimization problem consists in finding such shape of the boundary of the domain occupied by the body that the normal contact stress along the contact boundary of the body is minimized. In the paper shape as well as topological derivatives formulae of the cost functional are provided using material derivative and asymptotic expansion methods, respectively. These derivatives are employed to formulate necessary optimality condition for simultaneous shape and topology optimization. Level set based numerical algorithm for the solution of the shape optimization problem is proposed. Level set method is used to describe the position of the boundary of the body and its evolution on a fixed mesh. This evolution is governed by Hamilton – Jacobi equation. The speed vector field driving the propagation of the boundary of the body is given by the shape derivative of a cost functional with respect to the free boundary. Numerical examples are provided. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The calibration of volatility for option pricing models with jump diffusion processes

This paper is devoted to calibrate smooth local volatility surface under jump-diffusion processes. This calibration problem is posed as an inverse problem: given a finite set of observed European option prices, find a local volatility function such that the theoretical option prices matches the observed ones optimally with respect to a prescribed performance criterion. Firstly, we obtain an Euler-Lagrange equation for the calibration problem using Tikhonov regularization method. Then we solve the Euler–Lagrange equation using an iterative algorithm and obtain the volatility. Finally, numerical experiments show the effectiveness of the proposed method.     

Uzawa block relaxation domain decomposition method for a two-body frictionless contact problem

We propose a Uzawa block relaxation domain decomposition method for a two-body frictionless contact problem. We introduce auxiliary variables to separate subdomains representing linear elastic bodies. Applying a Uzawa block relaxation algorithm to the corresponding augmented Lagrangian functional yields a domain decomposition algorithm in which we have to solve two uncoupled linear elasticity subproblems in each iteration while the auxiliary variables are computed explicitly using Kuhn–Tucker optimality conditions.     

A simple method for spherical indentation of an elastic thin layer with surface tension bonded to a rigid substrate

An alternative method is proposed to solve the spherical indentation problem of an elastic thin layer with surface tension bonded to a rigid substrate. Based on the Kerr model, we establish a simple modified governing equation incorporating the surface tension effects for describing the relationship between the pressure and downward deflection of the impressed surface of the layer. This modified governing equation holds both inside and outside the contact zone, making it possible to analyze the whole layer by a unified differential equation. Numerical results are presented for the contact pressure inside the contact zone, the surface deflection of the elastic layer and the load-contact zone width relation to illustrate the present method. The validity and accuracy of the present method are demonstrated by comparing our results with those available in the existing literature.     

Block‐centered upwind multistep difference method and convergence analysis for numerical simulation of oil reservoir

The three‐dimensional displacement of two‐phase flow in porous media is a preliminary problem of numerical simulation of energy science and mathematics. The mathematical model is formulated by a nonlinear system of partial differential equations to describe incompressible miscible case. The pressure is defined by an elliptic equation, and the concentration is defined by a convection‐dominated diffusion equation. The pressure generates Darcy velocity and controls the dynamic change of concentration. We adopt a conservative block‐centered scheme to approximate the pressure and Darcy velocity, and the accuracy of Darcy velocity is improved one order. We use a block‐centered upwind multistep method to solve the concentration, where the time derivative is approximated by multistep method, and the diffusion term and convection term are treated by a block‐centered scheme and an upwind scheme, respectively. The composite algorithm is effective to solve such a convection‐dominated problem, since numerical oscillation and dispersion are avoided and computational accuracy is improved. Block‐centered method is conservative, and the concentration and the adjoint function are computed simultaneously. This physical nature is important in numerical simulation of seepage fluid. Using the convergence theory and techniques of priori estimates, we derive optimal estimate error. Numerical experiments and data show the support and consistency of theoretical result. The argument in the present paper shows a powerful tool to solve the well‐known model problem.     

The bending of a circular membrane on a linearly deformed foundation

The axisymmetric problem of the bending of a circular transversely-loaded membrane (i.e., a thin plate having no flexural stiffness), which lies without friction on a linearly deformed foundation, where there is contact over the whole area of the membrane, is considered. The problem is reduced to the combined investigation of a differential equation for the bending of the membrane and an integral equation of the first kind with an irregular kernel in the unknown contact pressure. The method of special orthonormalized polynomials and the regular asymptotic “large λ” method are used to solve the problem.     

Analysis and Numerical Approximation of a Contact Problem Involving Nonlinear Hencky-Type Materials with Nonlocal Coulomb’s Friction Law

A static frictional contact problem between an elasto-plastic body and a rigid foundation is considered. The material’s behavior is described by the nonlinear elastic constitutive Hencky’s law. The contact is modeled with the Signorini condition and a version of Coulomb’s law in which the coefficient of friction depends on the slip. The existence of a weak solution is proved by using Schauder’s fixed-point theorem combined with arguments of abstract variational inequalities. Afterward, a successive iteration technique, based on the Ka?anov method, to solve the problem numerically is proposed, and its convergence is established. Then, to improve the conditioning of the iterative problem, an appropriate Augmented Lagrangian formulation is used and that will lead us to Uzawa block relaxation method in every iteration. Finally, numerical experiments of two-dimensional test problems are carried out to illustrate the performance of the proposed algorithm.     

Steady-State Frictional Heat Generation on a Periodic Sliding Contact

We investigate the motion of a periodic system of rigid, isolated, parabolic dies along the surface of a half-space. The action of friction forces results in heat generation in the contact region. We assume that the surfaces of the dies are thermally insulated and the half-space is a heat conductor. We reduce the problem under consideration to a set of two integral equations for the contact temperature and pressure. We solve these equations numerically and investigate the influence of thermal deformation on the distributions of temperature and pressure and on the dimensions of the contact region.     

A matrix-free implementation of Riemannian Newton’s method on the Stiefel manifold

Newton’s method for unconstrained optimization problems on the Euclidean space can be generalized to that on Riemannian manifolds. The truncated singular value problem is one particular problem defined on the product of two Stiefel manifolds, and an algorithm of the Riemannian Newton’s method for this problem has been designed. However, this algorithm is not easy to implement in its original form because the Newton equation is expressed by a system of matrix equations which is difficult to solve directly. In the present paper, we propose an effective implementation of the Newton algorithm. A matrix-free Krylov subspace method is used to solve a symmetric linear system into which the Newton equation is rewritten. The presented approach can be used on other problems as well. Numerical experiments demonstrate that the proposed method is effective for the above optimization problem.     

求解“余数问题”的算法研究

求解"余数问题"可归结为一次同余式组x≡ri(mod pi)或一次不定方程组x=pixi+ri的求解.当方程的个数n与模pi(i=1,2,…,n)较大时,用同余式理论和孙子定理求解的过程非常繁琐.为此,运用试算分析法和辗转相除法,给出了求解上述问题的两种通用的计算机算法和程序.通过实践证明,该算法具有计算步骤简便,求解灵活快速,通用性强等优点.