Elastodynamic Contact Problems for Interface Cracks under Harmonic Loading

2-D fracture dynamics' problems for elastic bimaterials with cracks located at the bonding interface under the oblique time harmonic wave are considered in the study. The system of boundary integral equations for displacements and tractions is derived from Somigliana identity taking the contact interaction of the opposite crack faces into account. For the numerical solution the collocation method with piecewise constant approximation is used. The numerical results are obtained for various values of the angle of the wave incidence and the wave frequency taking the friction effects into account. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A periodic system of electrically permeable cracks at the interface between two piezoelectric materials

We constructed a closed solution for a bimaterial plane consisting of two dissimilar piezoelectric half-planes with a periodic system of electrically permeable cracks at the interface between these materials. The presence of zones of smooth contact of the crack lips near their tips was taken into account. By representing the characteristics of electromechanical fields via piecewise analytic functions, we reduced the problem to a Dirichlet–Riemann periodic problem, which was solved exactly. As a result of numerical analysis of the derived solution, we studied the dependence of the relative length of the contact zones and stress intensity factors on the ratio between the crack length and period for different combinations of piezoelectric materials.     

The solution of a problem in contact fracture mechanics on the nucleation and development of a bridged crack in the hub of a friction pair

The contact deformation of the hub of a plunger pair is considered. It is assumed that, during the repeated reciprocating motion of the plunger, a crack is initiated and fracture of the materials of the elements of the contact pair occurs. The problem of the equilibrium of the hub of a friction pair with a crack nucleus reduces to solving a system of non-linear singular integrodifferential equations with a Cauchy-type kernel. The normal and shear forces in the zone where the crack originates are found from the solution of this system of equations. The condition for the appearance of a crack is formulated, taking account of the criterion of the limit traction of the bonds in the material. A problem for the plunger of a friction pair as applied to a borehole sucker rod pump is considered as an example. In conclusion, the case when there are several arbitrarily distributed rectilinear bridged cracks, with bonds between the crack faces in the end zone, close to the contact surface of the hub is investigated.     

Development of generalized Iwan model to simulate frictional contacts with variable normal loads

Most friction models are originally proposed to predict restoring forces in mechanical contacts with constant normal load. In practice the contact interface kinematics may involve normal motion in addition to the tangential displacements, leading to variation of the contact normal load. This phenomenon is observed most strongly in contacts with high lateral vibration amplitudes and is known as slap. The current study establishes a general friction model to account for variation in the normal load and enables one to predict the behavior of a contact more precisely. Iwan model (1966) [5] is a suitable candidate for contact interface modeling and is able to represent the stick-micro/macro slip behavior involved in a friction contact. This physical based model is employed in the current work and its physical parameters are generalized to include the normal load variation effects. The model is characterized by a slippage distribution density function and a linear stiffness at stick state. Both these parameters, defined in presence of constant normal load in the original model, are derived considering normal load variation leading to generalization of the contact model. Conventional models with constant normal loads produce symmetric contact interface hysteresis loops, but the developed generalized Iwan model is capable of generating asymmetric hysteresis loops similar to those frequently seen in experiments. The generalized contact model is employed to simulate the measured behavior of a beam with frictional support observed in an experimental test set-up. The contact slippage distribution function is first identified in a constant normal load condition. Next in low levels of contact preloads where variation of the normal load is significant, the identified distribution function in generalized form is employed to predict the experimental observations.     

Problem of thermoelectromagnetoelasticity for a piezoelectric/piezomagnetic bimaterial with interface crack

We consider a plane strain problem for a piezoelectric/piezomagnetic bimaterial space with a crack in the region of the interface of the materials. At infinity, tensile and shear stresses and heat, electric, and magnetic flows are set. Using representations for all mechanical, thermal, and electromagnetic factors in terms of piecewise analytic functions, we formulate problems of linear conjugation that correspond to a model of an open crack and models taking into account the contact zone in the vicinity of a crack tip. Exact analytic solutions of the indicated problems are constructed. Expressions for stresses, the electric and magnetic inductions, jumps of derivatives of displacements, and electric and magnetic potentials on the interface are written. The coefficients of intensities of the indicated factors are presented. We derive a transcendental equation for the determination of the real length of the contact zone. The dependences of this length and the coefficients of intensity on the set external influences are investigated.     

Contact Problems for Interface Cracks under Harmonic Loading

The linear crack located at the bimaterial interface between dissimilar homogeneous linearly elastic materials under normally incident harmonic tension-compression wave is considered in the study. The problem is solved by the method of boundary integral equations using an iterative algorithm. The dynamic stress intensity factors are computed as functions of the loading frequency taking the contact interaction of the opposite crack's faces into account. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The calculation of the energy losses in a sliding elastic contact with friction and wear (the plane problem)

The plane problem of the sliding contact of a punch with an elastic foundation when there is friction and wear is considered. Assuming the existence of a steady solution in a moving system of coordinates, relations are derived between the sliding velocity, the wear, the contact stresses and the displacements for an arbitrary dependence of the wear rate on the contact pressure. Taking into account the presence of a deformation component of the friction force, an equation is written for the balance of the mechanical energy for the punch - elastic base system considered. It is shown that the equality of the work of the external force in displacing the punch to the losses due to friction and the change in the shape of the foundation due to wear is satisfied when the work done by the contact stresses on the increments of the boundary displacements is equal to zero, and the frictional losses must be determined taking into account the non-uniformity of the distributions of the shear contact stresses and the sliding velocity in the contact area. Two special cases of the foundation in the form of a wide and narrow strip are considered, for which the total coefficient of friction is calculated, taking into account the deformation component of the friction force.     

带裂纹十次对称二维准晶平面弹性的无摩擦接触问题

借助经典平面弹性复变函数方法，研究了单个刚性凸基底压头作用下，带任意形状裂纹十次对称二维准晶半平面弹性的无摩擦接触问题.利用十次对称二维准晶位移、应力的复变函数表达式, 带任意形状裂纹的准晶半平面弹性无摩擦接触问题被转换为可解的解析函数复合边值问题，进而简化成一类可解的Riemann边值问题.通过求解Riemann边值问题，得到了应力函数的封闭解, 并给出了裂纹端点处应力强度因子和压头下方准晶体表面任意点处接触应力的显式表达式.从压头下方接触应力的表达式可以看出, 接触应力在压头边缘和裂纹端点处具有奇异性.当忽略相位子场影响时, 该文所得结论与弹性材料对应结果一致.数值算例分别给出了单个平底刚性压头无摩擦压入带单个垂直裂纹和水平裂纹的十次对称二维准晶下半平面的结果.该文所得结论为准晶材料的应用提供了理论参考.     

Investigation of the prefracture zone at the tip of a mode I crack reaching a nonsmooth interface of an elastic media

We perform calculation of the initial prefracture zone at the tip of a mode I crack that reaches a nonsmooth interface of two dissimilar elastic media at its corner point by the Wiener–Hopf method. The zone is modeled by a line of normal displacement fracture on a crack continuation. Expressions for the length of the prefracture zone and the potential energy accumulated in it are obtained. Their numerical values are compared with the corresponding values for the prefracture zone in a bonding material on the interface of the media, on the basis of which we draw a conclusion on the possible direction of development of the zone.     

Contact interaction with friction of two elastic wedge-shaped bodies of different materials

The pleine static contact problem of elasticity theory concerning the impression of one wedge-shaped body into another of different material along sections of the side surfaces is examined. The abutting sections of both wedges start from the vertices. The problem is solved taking friction into account. In the case of greatest interest for applications, when the aperture angle of one of the wedges is π, an exact closed solution is constructed in the form of Cauchy-type integrals. However, the method of solution can be used for any wedge aperture angle.     

The asymptotic expansion of elastic fields in the vicinity of the contour of a plane crack at the interface of two materials

The stress-strain state in the neighbourhood of the front of a plane crack at the interface of two dissimilar half-spaces of ideally elastic isotropic materials is investigated. The form of the asymptotic expansions of the projections of the displacement vector onto the axis, directed along the tangent, the principal normal and the binormal to the crack contour is obtained. It is shown that asymptotic expansions of the projections of the displacement vector onto directions corresponding to the tangent and principal normal, beginning with the second-order term of the expansion, include both terms with half-integer and complex powers of the distances to the crack contour. This indicates that these projections of the solutions of the three-dimensional problem have singularities, defined by the solutions of both the antiplane and plane strain problems of cracks at the interface of materials. The singularities of the projection of the displacement vector on to the binormal correspond to the singularities of the solution of the plane strain problem.     

Two types of interface defects

The solution of a plane problem in the theory of elasticity for a two-component body with an interface, a finite part of which is either weakly distorted or is a weakly curved crack is constructed using the perturbation method. In the first case, it is assumed that the discontinuities in the forces and displacements at the interface are known, and, in the second case, the non-equilibrium nature of the load in the crack is taken into account. General quadrature formulae are derived for the complex potentials, which enable any approximation to be obtained in terms of elementary functions in many important practical cases. An algorithm is indicated for calculating each approximation. Families of defects are studied, the form of which is determined by power functions. The effect of the amplitude of the distortion and the shape of the interface crack on the Cherepanov–Rice integral as well as the shape of the distorted part of the interface on the stress concentration is investigated in the first approximation. An analysis of the applicability of the oscillating solution for a distorted interface crack is carried out. The results of the calculations are shown in the form of graphical relations.     

The indentation of a spherical punch into an ideally plastic half-space when there is contact friction

Calculations are presented of the indentation of a spherical punch into an ideally plastic half-space under condition of complete plasticity and taking account of contact friction, which is modelled according to Prandtl and Coulomb. Friction leads to the formation of a rigid zone at the centre of the punch when there is slipping of the material on the remaining part of the contact boundary. Limit values of the friction coefficients are obtained for which the rigid zone extends over the whole of the contact boundary. The dependence of the indentation force on the radius of the plastic area is in good agreement with experimental data.     

A. Yu. Ishlinskii's researches in the field of rolling friction and their development

The results of the research by A. Yu. Ishlinskii and others on the theory of rolling contact of deformable bodies, taking into account their imperfect elasticity and partial slip in the area contact region, are presented.     

The effect of microdamage on fatigue-crack growth

The theory of fatigue-crack growth, based on a synthesis of fracture mechanics and continuum mechanics of microdamage accumulation, is applied to the problem of crack growth under cyclic loading, taking into account the plastic deformations in the tip zones. The model of a thin plastic zone, which is a region of considerable plastic deformations at the crack tip, is supplemented by taking into account the effect of microdamage on the value of the specific fracture work and the limit stresses in the tip zone. Governing equations which describe fatigue-crack growth taking these factors into account are derived. The effect of the material characteristics and the load parameters on the growth rate and the distribution of microdamage in the tip zone and on its extensions is investigated by a computational experiment. Particular attention is given to the initial stage when crack growth may occur abruptly and the growth rate depends substantially on the initial conditions     

The pressure of a punch with a rounded edge on an elastic half-space

The problem of the unilateral contact without friction for a punch, the face of which is characterized by a rapid change in the neighbourhood of the a priori unknown boundary of the contact area, is investigated. Asymptotic formulae are obtained for the function which describes the variation of the contact area and the contact-pressure density in the boundary-layer region. The problem of the behaviour of the contact pressures in the neighbourhood of a smoothed stress concentrator is considered.     

Variational approach to contact problems in nonlinear elasticity

We use variational methods to study problems in nonlinear 3-dimensional elasticity where the deformation of the elastic body is restricted by a rigid obstacle. For an assigned variational problem we first verify the existence of constrained minimizers whereby we extend previous results. Then we rigorously derive the Euler-Lagrange equation as necessary condition for minimizers, which was possible before only under strong smoothness assumptions on the solution. The Lagrange multiplier corresponding to the obstacle constraint provides structural information about the nature of frictionless contact. In the case of contact with, e.g., a corner of the obstacle, we derive a qualitatively new contact condition taking into account the deformed shape of the elastic body. By our analysis it is shown here for the first time rigorously that energy minimizers really solve the mechanical contact problem. Received: 20 October 2000 / Accepted: 7 June 2001 / Published online: 5 September 2002     

Three-dimensional modelling of discrete particles by superellipsoids

Discrete Elements are used for the simulation of granular materials (sand, ballast) as well as for molecular assemblies. Circles (2D) and spheres (3D) are often used in literature on the Discrete Element Method (DEM) however they represent a strong idealisation of the real geometry. Superellipsoids provide the opportunity to generate a wide variety of three-dimensional geometrical shapes (e.g. sphere, cube, cylinder). The motion of each particle is described by means of rigid body dynamics. Suitable numerical integration methods are necessary which are able to conserve the essential physical quantities like momentum energy etc.. Possible choices are e.g. the explicit Verlet-Leapfrog method for the translation and the explicit fourth order Runge-Kutta method for the rotation. The implemented contact formulation takes damping as well as friction into account. Efficient implementation of the contact search is the main aim of this part of the work. It is subdivided into the neighbourhood search and the local search. A bisection algorithm is used to calculate the gap between two superellipsoids within the search. For the neighbourhood search two binning algorithms were implemented and compared for several packages of particles. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Physico-mechanical properties of polymer materials and the friction problem

Aspects of the friction behavior of polymer materials associated with their molecular structure and responsible for the relatively small value of the modulus of elasticity as compared with the ultimate strength are examined. It is shown that, in view of the important influence of hydrostatic pressure on the mechanical properties of polymer materials, the mechanical characteristics obtained from uniaxial testing cannot be used in calculating the contact area and the forces of friction. Formulas are presented for calculating the mechanical characteristics under these conditions. The results of indentation experiments designed to simulate contact processes are discussed. It is shown that the resistance to deformation of the asperities on the surface of polymer materials is of the order of the yield stress, and not two or three times greater, as with metals. The results of contact creep studies are described and evaluated. The results of investigations of the mechanical and antifriction properties of filled polymer materials show that the forces of friction are inversely proportional to the modulus of elasticity, while the thermophysical characteristics are a function not only of the thermophysical characteristics of the filler and the base, but also depend to a great extent on the shape of the filler particles; thus, when a fibrous filler with a low coefficient of linear expansion is used, the thermal stability of the friction material can be considerably improved. The results of a study of the adhesion interaction of polymer materials under conditions of omnidirectional nonuniform compression and simultaneous deformation are presented. It is shown that the adhesion interaction is strong even at room temperature. Aspects of the mechanical properties of lubricants that determine their effectiveness in polymer friction are considered.Institute of Polymer Mechanics, Academy of Sciences of the Latvian SSR, Riga. Translated from Mekhanika Polimerov, Vol. 5, No. 2, pp. 303–314, March–April, 1969.