A Numerical Three-Dimensional Model for the Contact of Rough Surfaces by Variational Principle

1996 ◽  
Vol 118 (1) ◽  
pp. 33-42 ◽  
Author(s):  
Xuefeng Tian ◽  
Bharat Bhushan

A new numerical method for the analysis of elastic and elastic-plastic contacts of two rough surfaces has been developed. The method is based on a variational principle in which the real area of contact and contact pressure distribution are those which minimize the total complementary potential energy. The present variational approach guarantees the uniqueness of the solution of the contact problem and significantly reduces the computation time as compared with the conventional matrix inversion method, and thus, makes it feasible to solve 3-D contact problem with large number of contact points. The model is extended to elastic-perfectly plastic contacts. The model is used to predict contact statistics for computer generated surfaces.

2000 ◽  
Vol 123 (2) ◽  
pp. 330-342 ◽  
Author(s):  
Wei Peng ◽  
Bharat Bhushan

A new numerical model for the three-dimensional contact analysis of a layered elastic–perfectly plastic half space with another rough surface is presented. The model is based on a variational principle in which the real area of contact and contact pressure distribution are those which minimize the total complementary potential energy. A quasi-Newton method is used to find the minimum. The influence coefficients matrix is determined using the Papkovich–Neuber potentials with fast Fourier transformation. The model is extended to elastic–perfectly plastic contacts in dry and wet conditions. Contact analyses have been performed to predict contact statistics of layered elastic/plastic solids with rough surfaces using this model. The effects of the stiffness of the layer and the substrate, layer thickness, as well as normal load are studied. Optimum layer parameters are identified to provide low friction/stiction and wear.


Author(s):  
Shaobiao Cai ◽  
Bharat Bhushan

A numerical three-dimensional contact model is presented to investigate the contact behavior of multilayered elastic-perfectly plastic solids with rough surfaces. The surface displacement and contact pressure distributions are obtained based on the variational principle with fast Fourier transform (FFT)-based scheme. Von Mises yield criterion is used to determine the onset of yield. The effective hardness is modeled and plays role when the local displacement meet the maximum displacement criterion. Simulations are performed to obtain the contact pressures, fractional total contact area, fractional plastic contact area, and surface/subsurface stresses. These contact statistics are analyzed to study the effects of the layer-to-substrate ratios of stiffness and hardness, surface roughness, and layers thickness of rough, two-layered elastic/plastic solids. The results yield insight into the effects of stiffness and hardness of layers and substrates, surface roughness, and applied load on the contact performance. The layer parameters leading to low friction, stiction, and wear are investigated and identified.


2018 ◽  
Vol 226 ◽  
pp. 04027 ◽  
Author(s):  
Alexander I. Sukhinov ◽  
Valentina V. Sidoryakina ◽  
Sofya V. Protsenko

This article is devoted to the confirmation the need for using a set of 3D dynamics models describing the various hydrophysical characteristics of the studied object to solve practical problems associated with the assessment of the ecological state of the water reservoirs. The present paper is devoted to the study of the three-dimensional model of transport and sedimentation of suspended matter in the coastal zone. The model takes into account such parameters as water movement, diffusionconvection, complicated bottom and shoreline geometry, lifting, transport and sedimentation of slurry. The existence and uniqueness of the solution of the corresponding indicated model of the initial-boundary value problem haas been envestigateded for two typical bottom boundary condirions. Also solution stability of the boundary-value problem in depend of functions: initial condition, boundary conditions and the righthand side in the norm L2 for any moment of time 0 < T < +∞, and also in the time-integral norm L2 has been proved. The model may be basis for the construction of hydrophysics models used to describe processes in the extraction of minerals from the seabed, in the dissemination of suspensions in shelf regions.


Author(s):  
Balázs Magyar ◽  
Bernd Sauer

In this paper, the influence of surface roughness on the local tribological load with a dry sliding contact is studied. First, three artificial rough surfaces with similar structure but different asperity heights are generated and projected on a smooth ball. After that, a contact pattern is determined between a rough ball and a smooth surface taking into account the elastic only as well as the linear elastic-perfectly plastic material description. On the basis of the calculated contact pressure distribution, the subsurface stresses and a three-dimensional temperature distribution in the sliding contact are calculated. The solutions show that a low surface roughness not necessarily results in low local tribological load of the surface.


Author(s):  
A. Sepehri ◽  
K. Farhang

In this paper we consider the contact between two rectangular rough surfaces that provide normal and tangential contact force as well as contact moment to counteract the net moment imposed by the applied forces. The surfaces are permitted to develop slight angular misalignment and thereby contact moment is derived. Through this scheme it is possible to also define elastic contribution to friction since the half-plane tangential contact force on one side of an asperity is no longer balanced by the half-plane tangential force component on the opposite side. The elastic friction force however is shown to be of a much smaller order than the contact normal force. Approximate closed form equations are found for contact force and moment as functions of separation, asperity radius of curvature sum, mean plane slope and nominal contact dimension. The approximate equations are shown to give error within seven percent.


Author(s):  
Tore Butlin ◽  
Alain Batailly

There is a drive towards minimising operating clearances within turbomachines in order to limit reverse leakage flows and hence improve their efficiency. This increases the likelihood of contact occurring between the blade and the casing, which can give rise to high amplitude vibration. Modelling this interaction represents a significant computational challenge. The non-linear contact precludes the use of well-established linear methods, and is also subject to uncertainties: the contact law is imprecisely known and the exact geometry of imperfections that trigger contact may be unknown. In this paper a novel approach is presented that aims to account for the uncertainties associated with the non-linearity in a non-probabilistic way. The worst case is sought, by framing the system as a constrained anti-optimisation problem. The target to be maximised represents a metric of the output of interest. The degrees of freedom of the anti-optimisation are the non-linear input forces (considered as external loads), and the constraints are designed to capture what is thought to be known about the non-linear contact law and geometry. A realistic three-dimensional model of a turbine blade is used to explore the approach, with contact considered at the leading and trailing edge. The blade dynamics are described in terms of a linear transfer function matrix and the target metric of interest is chosen to be the peak displacement of the contact points. The non-linearity is taken to result from an offset shaft, giving a sinusoidal clearance variation. The blade is driven at constant frequency and the scope of the study is limited to finding bounds on periodic solutions. A variety of constraint conditions are explored that describe aspects of the non-linearity. For example, only compressive forces are permitted (no tension from the contact), and the displacement must not exceed the clearance. The method yields encouraging initial results: constraints can be identified that give efficient estimates of the upper bound response of the system as a function of drive frequency. The results are compared with a benchmark time-domain simulation and are found to correctly over-predict the response without being overly conservative. Broad trends are also in agreement with the benchmark solution. The proposed method appears to be a promising approach for efficiently accounting for uncertainties associated with the non-linearity and thus improving blade design.


2021 ◽  
Vol 34 (04) ◽  
pp. 1148-1163
Author(s):  
Omid Amani Jafarlou ◽  
Afshin Ashja Ardalan ◽  
Ahmad Adib ◽  
Alireza Ganji ◽  
Soheila Bouzari

The airborne magnetic method is one of the most significantly applied approaches in mineral exploration. This magnetic method has been widely used in recent years due to its extensive coverage, lower cost than other methods and coverage of impassable areas and forests. For many years in Iran, the magnetic method has been utilized by authorities to survey ore deposits and to prepare structural and geological maps. This paper presents part of the Dehshir area exploration data of airborne surveys by helicopter with flight line spacing of 250 meters and flight height of 40 meters conducted by Geological Survey and Mineral Exploration of Iran. Furthermore, among the potential areas in the Dehshir area is the studied hereby area of Block C2 of Sourk mine. This region in Yazd province with a 4.1 km2 area is located at 95 km northwest of Yazd city. The purpose of this paper is to present a three-dimensional model of magnetic data and its validation using drilling data. The three-dimensional modeling of magnetic data has been processed by Inversion method. Next, the locations of boreholes were proposed for drilling. Then, the boreholes data is used to evaluate the geophysical method model by comparing the three-dimensional modeling of magnetic data and three-dimensional modeling of the proposed drilled boreholes data.


Author(s):  
D. L. Callahan

Modern polishing, precision machining and microindentation techniques allow the processing and mechanical characterization of ceramics at nanometric scales and within entirely plastic deformation regimes. The mechanical response of most ceramics to such highly constrained contact is not predictable from macroscopic properties and the microstructural deformation patterns have proven difficult to characterize by the application of any individual technique. In this study, TEM techniques of contrast analysis and CBED are combined with stereographic analysis to construct a three-dimensional microstructure deformation map of the surface of a perfectly plastic microindentation on macroscopically brittle aluminum nitride.The bright field image in Figure 1 shows a lg Vickers microindentation contained within a single AlN grain far from any boundaries. High densities of dislocations are evident, particularly near facet edges but are not individually resolvable. The prominent bend contours also indicate the severity of plastic deformation. Figure 2 is a selected area diffraction pattern covering the entire indentation area.


Sign in / Sign up

Export Citation Format

Share Document