A Class of Mathematical Programs with Equilibrium Constraints: A Smooth Algorithm and Applications to Contact Problems

We discuss a special mathematical programming problem with equilibrium constraints (MPEC), that arises in material and shape optimization problems involving the contact of a rod or a plate with a rigid obstacle. This MPEC can be reduced to a nonlinear programming problem with independent variables and some dependent variables implicity defined by the solution of a mixed linear complementarity problem (MLCP). A projected-gradient algorithm including a complementarity method is proposed to solve this optimization problem. Several numerical examples are reported to illustrate the efficiency of this methodology in practice.     

ELASTO-PLASTIC INTERFACE MODEL FOR 3D-FRICTIONAL ORTHOTROPIC CONTACT PROBLEMS

The present study deals with the solution of the fully three-dimensional contact/friction problem taking into account microstructural characteristics of the surfaces. An incremental non-associated hardening friction law model analogous to the classical theory of plasticity is used. Two different non-linear friction functions in the orthogonal directions are used to account for the orthotropic properties of the contacting bodies. A frontal solver processing unsymmetric matrices is adopted. Two numerical examples have been selected to show applicability of the method proposed. © 1997 by John Wiley & Sons, Ltd.     

Adaptive radial basis function–generated finite differences method for contact problems

This paper proposes an original adaptive refinement framework using radial basis function–generated finite differences method. Node distributions are generated with a Poisson disc sampling–based algorithm from a given continuous density function, which is altered during the refinement process based on the error indicator. All elements of the proposed adaptive strategy rely only on meshless concepts, which leads to great flexibility and generality of the solution procedure. The proposed framework is tested on four gradually more complex contact problems. First, a disc under pressure is considered and the computed stress field is compared to the closed-form solution of the problem to assess the behaviour of the algorithm and the influence of free parameters. Second, a Hertzian contact problem is studied to analyse the proposed algorithm with an ad hoc error indicator and to test both refinement and derefinement. A contact problem, typical for fretting fatigue, with no known closed-form solution is considered and solved next. It is demonstrated that the proposed methodology produces results comparable with finite element method without the need for manual refinement or any human intervention. In the last case, generality of the proposed approach is demonstrated by solving a three-dimensional Boussinesq's problem.     

SEMIANALYTICAL SINGULAR ELEMENT AND ITS APPLICATION TO STRESS CALCULATION AND OPTIMIZATION

A semianalytical singular element is proposed to model the stress fields in plane and axisymmetric elastic bodies in the vicinity of singular points. The element generation is performed taking into account the structure of analytical solution in the vicinity of such points. This renders the element well suited for any type of singular points. Introduction of specially generated shape functions in the approximating expressions for displacements enables one to ensure the displacement compatibility with the standard elements. A detailed consideration is given to the convergence of the finite-element procedure. The extensive possibilities of the proposed singular element are exemplified by solving particular elastic problems. The applicability of the elaborated singular element is discussed with respect to optimization of geometry and material properties in the vicinity of singular points. Two optimization problems are presented to illustrate a search for the body geometry providing the optimal stress distribution.     

全动操纵面系统中摩擦参数测量方法研究

通过解析方法和数值方法研究了库仑摩擦对结构动响应的影响规律。利用库仑摩擦在时域上能够傅里叶展开为方波函数的特点,得到含库仑摩擦的单自由度系统在简谐激励下的级数响应,结果表明,库仑摩擦的引入导致结构的频域响应中,除了激振频率成分外,还会有奇数倍频的成分,在激振频率为0.605倍固有频率时,奇数倍频成分与激振频率成分的幅值比例是摩擦力的线性函数。通过数值方法证明,上述规律能够推广到二自由度的全动操纵面模型中。利用这一规律,通过动响应试验的方法测量出全动操纵面系统中的库仑摩擦参数,并验证了这种方法的正确性。     

Energy consistent algorithms for frictional contact problems

In this paper, the energy and momentum conserving algorithmic paradigm is extended to encompass a phenomenon featuring physical dissipation: dynamic frictional contact. Whereas in other works dealing with conservative systems the chief aim is often the maintenance of numerical stability in the non-linear regime, in this investigation we seek to achieve not only this benefit but also the accurate algorithmic production of physical dissipation associated with frictional processes. The approach here features a product formula algorithm for the evolution of local frictional conditions, with the associated operator split guided by an a priori energy estimate. The resulting algorithm is characterized by exact conservation of energy during stick friction, and positive dissipation consistent with the frictional model used during slip. Effectiveness of the algorithm is demonstrated by a series of finite element simulations involving large deformations and frictional slip, complete with appropriate comparisons to more traditional schemes. © 1998 John Wiley & Sons, Ltd.     

A novel method for integrating first‐ and second‐order differential equations in elastohydrodynamic lubrication for the solution of smooth isothermal,line contact problems

This paper contains details of recent developments in the analysis of elastohydrodynamic lubrication problems using the finite element method. A steady state isothermal finite element formulation of the smooth line contact problem with Newtonian lubricant behaviour is presented containing both first‐ and second‐order formulations of the hydrodynamic equation. Previous problems with the limited range of applicability of both first‐ and second‐order finite difference solutions have been overcome by summing both the first‐ and second‐order equations' weighted contributions. Application of the method to a range of problems unattainable by either single first‐ or second‐order formulations is presented. Copyright © 1999 John Wiley & Sons, Ltd.