Solution Of Contact Problem Research Articles

Modeling Solutions of Axisymmetric Contact Problems Using Dimensionality Reduction Method with or without Adhesion

Axisymmetric contact problems of rigid bodies with elastomeric semi three-dimensional environment have been studied. It has been shown that three-dimensional contact problems with or without adhesion can be mapped with the dimensionality reduction method as one-dimensional problems with minimum error. The efficiency of the method has been illustrated by numerical calculations of practical problems. This also confirms the correctness of the analytical estimators.

Численное решение динамической задачи контакта упругопластических тел

Численное решение динамической задачи контакта упругопластических тел

An analytical solution of contact problem of indentation for multi-layered material

An analytical solution of contact problem of indentation for multi-layered material

Contact Problem for a Viscoelastic Half-Space and a Rough Rigid Cylinder under the Conditions of Hydrodynamic Lubrication

A model has been proposed for use in studying the combined effect of the roughness of a rigid punch and the viscous properties of a base separated by a thin lubricant layer exerted on the characteristics of contact interactions and the sliding friction force. The problem of the motion of a thin a lubricant layer between a fixed rigid cylinder with a regular relief, as well as the surface of a moving viscoelastic half-space, the rheological properties of which are described by an integral operator with an exponential creep kernel, has been considered. The pressure and thickness of the lubricating layer, as well as the deformation component of the frictional force depending on the sliding velocity, have been analyzed. A comparison of the results of solutions of contact problems for viscoelastic and elastic rough bodies in the presence of a lubricant has been presented.

Approximated Analytical Solution of a Problem on Indentation of an Electro-Elastic Half-Space with Inhomogeneous Coating by a Conductive Punch

The paper addresses to the development of the bilateral asymptotic method of solution of contact problems taking into account electro-mechanical fields coupling. Approximated analytical solutions are provided for contact problems on indentation of an electroelastic piezoelectric half-space with inhomogeneous in depth piezoelectric coating by a rigid conductive circular flat, spherical or conical punch.

The projected Barzilai–Borwein method with fall-back for strictly convex QCQP problems with separable constraints

A variant of the projected Barzilai–Borwein method for solving the strictly convex QCQP problems with separable constraints is presented. The convergence is enforced by a combination of the fall-back strategy and the fixed step-length gradient projection. Using the recent results on the decrease of the convex quadratic function along the projected-gradient path, we prove that the algorithm enjoys the R-linear convergence. The algorithm is plugged into our scalable TFETI based domain decomposition algorithm for the solution of contact problems and its performance is demonstrated on the solution of contact problems, including a frictionless problem and the problems with the isotropic and orthotropic Tresca friction.

A stable and optimally convergent LaTIn‐CutFEM algorithm for multiple unilateral contact problems

SummaryIn this paper, we propose a novel unfitted finite element method for the simulation of multiple body contact. The computational mesh is generated independently of the geometry of the interacting solids, which can be arbitrarily complex. The key novelty of the approach is the combination of elements of the CutFEM technology, namely, the enrichment of the solution field via the definition of overlapping fictitious domains with a dedicated penalty‐type regularisation of discrete operators and the LaTIn hybrid‐mixed formulation of complex interface conditions. Furthermore, the novel P1‐P1 discretisation scheme that we propose for the unfitted LaTIn solver is shown to be stable, robust, and optimally convergent with mesh refinement. Finally, the paper introduces a high‐performance 3D level set/CutFEM framework for the versatile and robust solution of contact problems involving multiple bodies of complex geometries, with more than 2 bodies interacting at a single point.

Solution of contact problems for Gao beam and elastic foundation

This paper presents mathematical formulations and a solution for contact problems that concern the nonlinear beam published by Gao (Nonlinear elastic beam theory with application in contact problems and variational approaches, Mech Res Commun 1996; 23: 11–17) and an elastic foundation. The beam is subjected to a vertical and also axial loading. The elastic deformable foundation is considered at a distance under the beam. The contact is modeled as static, frictionless and using the normal compliance contact condition. In comparison with the usual contact problem formulations, which are based on variational inequalities, we are able to derive for our problem a nonlinear variational equation. Solution of this problem is realized by means of the so-called control variational method. The main idea of this method is to transform the given contact problem to an optimal control problem, which can provide the requested solution. Finally, some results including numerical examples are offered to illustrate the usefulness of the presented solution method.

Numerical solution of smooth and rough contact problems

A general procedure for the numerical solution of three-dimensional normal and tangential contact problems for non-conforming bodies of arbitrary shape is presented. The proposed procedure is based on the interpolation of known analytical solutions relating surface displacements of an elastic half-space subjected to pyramidal distributions of normal and tangential surface tractions. These formulas are used to interpolate unknown pressure distributions on the contact zone of two elastic bodies in contact. The proposed interpolation significantly reduces the computational burden of the numerical procedure required to interpolate the actual pressure distribution and to determine the initially unknown contact region. It amounts to expressing the contact pressure as a function of the elastic parameters of the two bodies, the distance between their surfaces and the relative displacement between two far-points pertaining to the bodies in contact. The procedure has been validated by comparison with classical contact problems and the results show excellent agreement with existing analytical and numerical solutions.

Численное решение контактной задачи теории упругости с односторонними связями с помощью смешанной схемы метода конечных элементов

Численное решение контактной задачи теории упругости с односторонними связями с помощью смешанной схемы метода конечных элементов

New Solution of Axisymmetric Contact Problem of Elasticity

We consider two problems of elasticity, namely, the Boussinesq problem about the action of a lumped force on a half-space and the related problem about the interaction of the half-space with a cylindrical rigid punch with plane base. In the classical statement, these problems have singular solutions. In the Boussinesq problem, the displacement under the action of the force is infinitely large, and in the punch problem, the infinitely large variable is the pressure on the punch boundary. In the present paper, these problems are solved with the use of relations of generalized elasticity derived regarding a medium element of small but finite dimensions rather than a traditional infinitesimal element. The structure parameter of the medium contained in the solutions can be determined experimentally. The obtained generalized solutions of the problems under study are regular.

A Note on Massively Parallel Implementation of FETI for the Solution of Contact Problems

The paper deals with the solution of large multibody contact problems using massively parallel computers and domain decomposition methods. These methods can solve the problems discretized by billions of nodal variables at the cost nearly proportional to the number of variables using up to thousands cores before the communication costs start to dominate the computational costs. The paper describes the ingredients essential for efficient massively parallel implementation that increases the parallel scalability beyond the limit mentioned above. The improvements were enhanced into a new software package PERMON which is based on PETSc. The performance of the algorithm is demonstrated on the solution of an academic benchmark discretized by nearlybillion of nodal variables.

On the Efficient Reconstruction of Displacements in FETI Methods for Contact Problems

The final step in the solution of contact problems of elasticity by FETI-based domain decomposition methods is the reconstruction of displacements corresponding to the Lagrange multipliers for ''gluing'' of subdomains and non-penetration conditions. The rigid body component of the displacements is usually obtained by means of a well known but quite complex formula, the application of which requires reassembling and factorization of some large matrices. Here we propose a simple formula which is applicable to many variants of the FETI based algorithms for contact problems. The method takes a negligible time and avoids reassembling or factorization of any matrices.

Singular solutions of contact problems and block elements

In this work, we consider mixed problems of elasticity theory, in particular, contact problems for cases that are nontraditional. They include mixed problems with discontinuous boundary conditions in which the singularities in the behavior of contact stresses are not studied or the energy of the singularities is unbounded. An example of such mixed problems is contact problems for two rigid stamps approaching each other by rectilinear boundaries up to contact but not merging into one stamp. It has been shown that such problems, which appear in seismology, failure theory, and civil engineering, have singular components with unbounded energy and can be solved by topological methods with pointwise convergence, in particular, by the block element method. Numerical methods that are based on using the energy integral are not applicable to such problems in view of its divergence.

Semi-analytic contact technique in a non-linear parametric model order reduction method for gear simulations

In this work we present a novel method for the solution of gear contact problems in flexible multi-body. These problems are characterized by significant variation in the location and size of the contact area, typically requiring a high number of degrees of freedom to correctly capture deformation and stress fields. Therefore fully dynamic simulation is computationally prohibitive. To overcome these limitations, we exploit a combined analytic-numerical contact model within a parametric model order reduction (PMOR) scheme. The reduction space consists of a truncated set of eigenvectors augmented with a parameter dependent set of residual static shape vectors. Each static shape is computed by interpolating among a set of displacement modes of the interacting bodies, obtained from a series of precomputed static contact analyses. During the contact analyses, an analytic model based on the Hertz theory describes the teeth local deformation. We implement the proposed method in an in-house code and we apply it to spur and helical gears dynamic contact analyses. We compare the results with classical PMOR schemes highlighting how the combined use of the semi-analytic contact model allows to decrease further the model complexity as well as the computational burden, for both static and dynamic cases. Finally, we validate the methodology by means of a comparison with experimental data found in literature, showing that the numerical method is able to capture quantitatively the static transmission error measurements in case of both helical and spur geared transmission for different torque levels.

Lubricated point heavily loaded contacts of functionally graded materials. Part 1. Dry contacts

Over the last couple decades coatings attract more and more attention in practical applications. The present study addresses a question which is not well studied: how coated surfaces behave in lubricated contacts? In other words, this is a study of the effectiveness of functionally graded materials in heavily loaded point elastohydrodynamically lubricated contacts with straight lubricant entrainment. As a part of the study, some criteria of coating effectiveness are introduced and discussed. More specifically, the behavior of main parameters such as the lubrication film thickness and the frictional force in point heavily loaded lubricated contacts of functionally graded elastic materials are considered. The problem is studied based on the method of matched asymptotic expansions which allows us to split the problem into two separate problems: a dry contact problems for functionally graded elastic materials and an elastohydrodynamically lubricated problem for functionally graded materials. The elastohydrodynamically lubricated problem uses as input data not only the operational and physical parameters of the materials and lubricant but also the asymptotic behavior of the dry contact problem solution near the contact boundaries. Therefore, a sequence of two problems must be solved: the dry contact problems for functionally graded elastic materials and the elastohydrodynamically lubricated problem for functionally graded materials. Similar methods have been used for the analysis of an elastohydrodynamically lubricated problem for heavily loaded line contacts of functionally graded materials. The dry contact problem will be analyzed in Part 1 of the paper based on a semi-analytical bilateral method which produces correct asymptotic solutions for thin and thick coatings. The analytical expressions for contact pressure are obtained and analyzed for various combinations of coating thicknesses and elastic properties. The elastohydrodynamically lubricated problem will be considered in Part 2 of the paper based on the method of matched asymptotic expansions. In the analysis of the elastohydrodynamically lubricated problem, as in the case of homogeneous contact materials, it is shown that the whole contact region can be subdivided into three subregions: the central one which is adjacent to the other two regions occupied by the ends of the zones. The central region can be subdivided into the Hertzian region and then adjacent to it inlet and exit zones which, in turn, are adjacent to the inlet and exit boundaries of the contact, respectively. In the Hertzian region the elastohydrodynamically lubricated problem solution is very close to the solution of the corresponding dry (i.e. non-lubricated) contact problem for functionally graded elastic materials which have been analyzed. In the central region in the inlet and exit zones of a heavily loaded point elastohydrodynamically lubricated contact, the elastohydrodynamically lubricated problem for functionally graded elastic materials using certain scaling transforms can be reduced to asymptotically valid equations identical to the ones obtained in the inlet and exit zones of heavily loaded line elastohydrodynamically lubricated contacts for homogeneous elastic materials. Therefore, many of the well known properties of heavily loaded line elastohydrodynamically lubricated contacts for homogeneous elastic materials are also valid for heavily loaded point elastohydrodynamically lubricated contacts for functionally graded elastic materials. These asymptotically valid equations can be analyzed and numerically solved based on stable methods using a specific regularization approach, which were developed for lubricated line contacts. Also, this asymptotic analysis leads to an easy analytical derivation of formulas for the lubrication film thickness which take into account the inhomogeneity of the elastic materials. As a result of this analysis, some criteria for lubrication film thickness increase and friction force reduction are proposed. These criteria depend on lubricant properties as well as the properties of functionally graded elastic materials involved in lubricated contacts. Such a sequential solution of the elastohydrodynamically lubricated problem for functionally graded materials makes it perfectly clear what the dependence is of elastohydrodynamically lubricated contact parameters on the solid material (including the coating) and lubricant properties.

Lubricated point heavily loaded contacts of functionally graded materials. Part 2. Lubricated contacts

The paper presents a study of the effectiveness of functionally graded materials in heavily loaded point elastohydrodynamically lubricated contacts with straight lubricant entrainment. As part of the study, some criteria of coating effectiveness are introduced and discussed. The approach used for achieving this goal is based on the recently developed asymptotic solution method for steady isothermal elastohydrodynamically lubricated problems of heavily loaded point contacts, as well as on a semi-analytical method for the solutions of contact problems for elastic functionally graded materials without the presence of lubrication in the contacts (dry contacts). For line elastohydrodynamically lubricated contacts such a problem has already been studied. The solution process can be subdivided into two stages: the solution of the contact problems for dry contacts and a solution of the problem for lubricated contacts. The solution of dry contact problems for functionally graded materials is described in detail in Part 1 of the paper. As in the case of homogeneous contact materials in the analysis of the elastohydrodynamically lubricated problem for functionally graded elastic materials, it is shown that the whole contact region can be subdivided into three subregions: the central one which is adjacent to the other two regions occupied by the ends of the zones. The central region can be subdivided into the Hertzian region and adjacent to it the inlet and exit zones which, in turn, are adjacent to the inlet and exit boundaries of the contact, respectively. In the Hertzian region the elastohydrodynamically lubricated problem solution is very close to the solution of the corresponding dry (i.e. non-lubricated) contact problem for functionally graded elastic materials which have been analyzed in Part 1 of the paper. In the central region in the inlet and exit zones of the heavily loaded point elastohydrodynamically lubricated contact, the elastohydrodynamically lubricated problem for functionally graded elastic materials using certain similarity transforms can be reduced to asymptotically valid equations identical to the ones obtained in the inlet and exit zones of the heavily loaded line elastohydrodynamically lubricated contacts for homogeneous elastic materials. Therefore, many of the well known properties of heavily loaded line elastohydrodynamically lubricated contacts for homogeneous elastic materials are also valid for heavily loaded point elastohydrodynamically lubricated contacts for functionally graded elastic materials. These asymptotically valid equations can be analyzed and numerically solved based on stable methods using a specific regularization approach which was developed for lubricated line contacts. Also, this asymptotic analysis leads to the easy analytical derivation of equations for the lubrication film thickness which take into account the inhomogeneity of the elastic materials. As a result of this analysis, criteria for an increase in lubrication film and friction force reduction are proposed. These criteria depend on lubricant properties, as well as the properties of functionally graded elastic materials involved in lubricated contacts.

Contact problem for an orthotropic half-space

Numerical and analytical solutions of the 3D contact problem of elasticity on the penetration of a rigid punch into an orthotropic half-space are obtained disregarding the friction forces.A numericalmethod ofHammerstein-type nonlinear boundary integral equations was used in the case of unknown contact region, which permits determining the contact region and the pressure in this region. The exact solution of the contact problem for a punch shaped as an elliptic paraboloid was used to debug the program of the numerical method. The structure of the exact solution of the problem of indentation of an elliptic punch with polynomial base was determined. The computations were performed for various materials in the case of the penetration of an elliptic or conical punch.

Implementation of asymptotics of the wear contact problem solution for identifying the wear law based on the results of tribological tests

The parameter of the wear law of bronze coupled with steel under dry friction is identified by comparing the results of mathematical modeling with experimental data. The wear tests were carried out according to the steel-ball–flat-bronze sample scheme under relative reciprocal sliding of the tested samples. The modeling of the wear of the sample was based on asymptotic behavior of the solution of the corresponding wear-contact problem. The effectiveness of using this asymptotic behavior to identify the wear law based on the tribological test results has been demonstrated.

Closed-form solution of axisymmetric contact problem for a viscoelastic base within cycle of increasing and decreasing of load on the indenter

The axisymmetric contact problem for a spherical indenter and viscoelastic base has been considered within a cycle of increasing and decreasing of load. The loading-unloading cycle has been performed at a constant speed of indenter movment. The quasi-static approximation has been used. The study has been carried out according to linear viscoelasticity theory, while the material model allows one to take into account an arbitrary number of relaxation times. The dependence of load on the indenter on time has been defined. Expressions for determining the unknown boundary of the contact area in the loading–unloading process has been given. The hysteretic loop area at different values of the material relaxation time has been analyzed. The applicability of the offered analytical model of elastomer indentation experiments data processing has been discussed.