点群10mm十次对称二维准晶中的两类接触问题

采用复变函数法探讨了在一个刚性压头作用下十次对称二维准晶材料的两类接触问题,即具有有限摩擦的接触问题以及粘结接触问题.特别地对于平底压头,获得了表征声子场和相位子场的全纯函数的显式表达式,以及在压头上的接触应力分布.结果显示,对于具有有限摩擦的接触问题,接触应力在接触区边缘具有实指数奇异性-1／2±β,其中β由准晶体的材料常数及静摩擦系数确定;而对于粘结接触问题,接触应力在接触区边缘具有振荡型奇异性-1／2±iε,其中ε由准晶体的材料常数确定.     

Two Kinds of Contact Problems in Dodecagonal Quasicrystals of Point Group 12 mm

In this paper,two kinds of contact problems in 2-D dodecagonal quasicrystals were discussed using the complex variable function method:one is the finite frictional contact problem,the other is the adhesive contact problem.The analytic expressions of contact stresses in the phonon and phason fields were obtained for a flat rigid punch,which showed that:(1) for the finite frictional contact problem,the contact stress exhibited power-type singularities at the edge of the contact zone;(2) for the adhesive contact problem,the contact stress exhibited oscillatory singularities at the edge of the contact zone.The distribution regulation of contact stress under punch was illustrated;and the low friction property of quasicrystals was verified graphically.     

八次对称二维准晶材料接触问题

本文通过引入位移函数和应用Fourier分析与对偶积分方程理论圆满解决了在一个刚性平头冲头作用下八次对称二维准晶材料的接触问题，得到了此材料接触问题应力与位移的解析表达式。结果表明，如果接触位移在接触区域内为一常数，则接触应力在接触边缘具有1／2阶奇异性，这为准晶材料的接触变形提供了重要的力学量。     

Axisymmetric contact problem of cubic quasicrystalline materials

The axisymmetric elasticity theory of cubic quasicrystal was developed in Ref. [1]. The axisymmetric elasticity problem of cubic quasicrystal is reduced to a single higher-order partial differential equation by introducing a displacement function, based on which, the exact analytic solutions for the elastic field of an axisymmetric contact problem of cubic quasicrystalline materials are obtained for universal contact stress or contact displacement. The result shows that if the contact stress has order −1/2 singularity on the edge of the contact domain, the contact displacement is a constant in the contact domain. Conversely, if the contact displacement is a constant, the contact stress must have order −1/2 singularity on the edge of the contact domain. Project supported by the National Natural Science Foundation of China (No. 19972011).     

Steady sliding frictional contact problem for a 2d elastic half-space with a discontinuous friction coefficient and related stress singularities

The steady sliding frictional contact problem between a moving rigid indentor of arbitrary shape and an isotropic homogeneous elastic half-space in plane strain is extensively analysed. The case where the friction coefficient is a step function (with respect to the space variable), that is, where there are jumps in the friction coefficient, is considered. The problem is put under the form of a variational inequality which is proved to always have a solution which, in addition, is unique in some cases. The solutions exhibit different kinds of universal singularities that are explicitly given. In particular, it is shown that the nature of the universal stress singularity at a jump of the friction coefficient is different depending on the sign of the jump.     

一般各向异性单侧接触界面上波的反射和折射

研究简谐弹性波在一般各向异性介质单侧接触界面上的反射和折射问题．利用Fouier分析方法将非线性Coulomb摩擦接触边界波动问题化为一组代数方程．给出了确定局部分离、滑移和粘着区的思路和方法及各区域的解；讨论了出现界面局部分离和滑移的条件．对特定材料组合情况进行了详细数值计算，给出了界面力、相对滑移速度、张开位移、高频谐波的反射折射系数等特征参量；考察了平面和反平面波动的耦合及整体滑移等．其中关于高频谐波的结果可对已有实验结果给出很好的定性解释．在大多数情况下，即使对摩擦系数无穷大的粘滞接触界面，分离区端部也总是存在一个很小的滑移区。     

铁路车轴过盈配合面微动损伤分析及有限元仿真

基于全尺寸铁路车轴疲劳试验，观察并分析了微动区损伤形貌及损伤机理. 基于测量的磨损轮廓建立有限元模型，计算分析了微动磨损对过盈配合面微动参量及轴向应力的影响. 结果表明:轮座近加载侧存在1个宽度约为20 mm的微动损伤区，根据形貌特征可以分为3个区域. 仿真得到的微动滑移区宽度与损伤区宽度基本一致，张开区宽度略小于磨损区. 未磨损时，接触压应力、摩擦剪应力及轴向应力峰值均出现在接触最边缘；在微动磨损作用下，接触压应力、摩擦剪应力、轴向应力峰值出现在磨损-未磨损边界，且轴向应力数值在磨损区由负变正. 磨屑的存在为接触面提供承载平台，在一定程度上抑制应力集中向内部转移.      

十二次对称二维准晶中的无摩擦接触问题

利用积分变换的方法讨论了在一个刚性压头作用下十二次对称二维准晶的无摩擦接触问题. 通过引入位移势函数,将数量巨大而复杂的偏微分方程转化为两个独立的双调和方程,应用Fourier分析与对偶积分方程理论解决了十二次对称二维准晶材料的无摩擦接触问题,得到了相应的接触应力解析表达式,结果表明：如果接触位移是一常数,则接触应力在接触区域边缘具有-1/2阶奇异性;反之,如果接触应力在接触区域边缘具有-1/2阶的奇异性,则接触位移一定为一常数,这为准晶材料的接触变形提供了重要的力学参数.     

Approximation of thermoelasticity contact problem with nonmonotone friction

The paper presents the formulation and approximation of a static thermoelasticity problem that describes bilateral frictional contact between a deformable body and a rigid foundation. The friction is in the form of a nonmonotone and multivalued law. The coupling effect of the problem is neglected. Therefore, the thermic part of the problem is considered independently on the elasticity problem. For the displacement vector, we formulate one substationary problem for a non-convex, locally Lipschitz continuous functional representing the total potential energy of the body. All problems formulated in the paper are approximated with the finite element method.     

New analytical and numerical results for two-dimensional contact profiles

Recently, a generalized Coulomb law for elastic bodies in contact has been developed by the author, which assumes that the tangential traction is the difference of the slip stress of the contact and the stick area, whereby each stick area corresponds to a smaller contact area. It holds for multiple contact regions also. Several applications for elastic half planes, half spaces, thin and thick layers and impact problems have been published. For plane contact of equal bodies with friction, it provides exact solutions, and the interior stress field can be expressed with analytical results in closed form. In this article, a singular superposition of flat punch solutions is outlined, in which the punches are aligned with an edge of the contact area. It is shown that this superposition satisfies Coulomb's inequalities directly, and new results for the Muskhelishvili potentials of several profiles are presented. It is illustrated how problems of singularity and multi-valuedness of complex functions can be solved in closed form, and the Chebyshev approximation used by earlier authors can be avoided. For comparison, some previous solutions for symmetric profiles are appended. Some results for the interior stress field, the pressure, the frictional traction and the surface displacements are compared with FEM solutions of an equivalent problem. The small differences between both methods show characteristic features of the FEM model and the theoretical assumptions, and are shortly explained. Further, this example can be used as benchmark test for FEM and BEM programs.     

Slip-line modeling of machining with a rounded-edge tool—Part II: analysis of the size effect and the shear strain-rate

Part II of the present study quantitatively analyzes orthogonal metal cutting processes based on the new slip-line model proposed in Part I. The applicable range of the model is illustrated, followed by an explanation of the non-unique nature of the model. It is suggested that the tool edge roundness be comprehensively defined by four variables. Namely: tool edge radius, position of the stagnation point on the tool edge, tool-chip frictional shear stress above the stagnation point on the tool edge, and tool-chip frictional shear stress below the stagnation point on the tool edge. The effects of these four variables on eight groups of machining parameters are investigated. These include (1) cutting force, thrust force, resultant force, and the ratio of cutting force to thrust force; (2) ploughing force; (3) chip up-curl radius; (4) chip thickness; (5) tool-chip contact length; (6) thickness of the primary shear zone; (7) average shear strain in the primary shear zone; and (8) average shear strain-rate in the primary shear zone. The importance of tool edge roundness is further reinforced by a series of new research findings made in this paper. It is revealed that the size effect highly depends on the material constitutive behavior in machining. The dependence of the thickness of the primary shear zone and the dependence of the magnitude of shear strain-rate in the primary shear zone on the tool edge radius are well demonstrated. A surprisingly good agreement between theory and experiments is reached.     

Double frictional receding contact problem for an orthotropic layer loaded by normal and tangential forces

With the increasing research in the field of contact mechanics, different types of contact models have been investigated by many researchers by employing various complex material models. To ascertain the orthotropy effect and modeling parameters on a receding contact model, the double frictional receding contact problem for an orthotropic bilayer loaded by a cylindrical punch is taken into account in this study. Assuming plane strain sliding conditions, the governing equations are found analytically using Fourier integral transformation technique. Then, the resulting singular integral equations are solved numerically using an iterative method. The weight function describing the asymptotic behavior of the stresses are investigated in detail and powers of the stress singularities are provided. To control the trustworthiness and correctness of the analytical formulation and to compare the resulting stress distributions and contact boundaries, a numerically efficient finite element method was employed using augmented Lagrange contact algorithm. The aim of this paper is to investigate the orthotropy effect, modeling parameters and coefficients of friction on the surface and interface stresses, surface and interface contact boundaries, powers of stress singularities, weight function and to provide highly parametric benchmark results for tribological community in designing wear resistant systems.

     

First mushroom-shaped adhesive microstructure: A review

This letter reviews the adhesive and frictional properties of the first mushroom-shaped adhesive microstructure (MSAMS), which has come a long way from inspiration by the attachment devices evolved in beetles to a large-scale industrial production. It was shown to have an that about twice higher pull-off force compared to a smooth control made from the same material measured on smooth substrates. Pull-off forces measured underwater are even higher than those in air. Moreover, it retained adhesive performance over thousands of attachment cycles and initial adhesive capability could be recovered by washing after being contaminated. In shearing, MSAMS exhibits reduced and stabilized friction in comparison with a smooth control, which demonstrated pronounced stick-slip motion, and shows zero pull-off force in a sheared state, allowing the adhesion to be switched on and off. The presence of a fluid in the contact zone showed adhesion enhancement on both smooth and rough substrates. All these features lead us to conclude that MSAMS may have practical potential in a variety of applications.     

Multiple-zone sliding contact with friction on an anisotropic thermoelastic half-space

The paper studies a class of multiple-zone sliding contact problems. This class is general enough to include frictional and thermal effects, and anisotropic response of the indented material. In particular, a rigid die (indenter) slides with Coulomb friction and at constant speed over the surface of a deformable and conducting body in the form of a 2D half-space. The body is assumed to behave as a thermoelastic transversely isotropic material. Thermoelasticity of the Green–Lindsay type is assumed to govern. The solution method is based on integral transforms and singular integral equations. First, an exact transform solution for the auxiliary problem of multiple-zone (integer n > 1) surface tractions is obtained. Then, an asymptotic form for this auxiliary problem is extracted. This form can be inverted analytically, and the result applied to sliding contacts with multiple zones. For illustration, detailed calculations are provided for the case of two (n = 2) contact zones. The solution yields the contact zone width and location in terms of sliding speed, friction, die profile, and also the force exerted. Calculations for the hexagonal material zinc illustrate effects of speed, friction and line of action of the die force on relative contact zone size, location of maximal values for the temperature and the compressive stress, and the maximum temperature for a given maximum stress. Finally, from our general results, a single contact zone solution follows as a simple limit.     

Two-dimensional sliding frictional contact of functionally graded materials

A multi-layered model for sliding frictional contact analysis of functionally graded materials (FGMs) with arbitrarily varying shear modulus under plane strain-state deformation has been developed. Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into several sub-layers and in each sub-layers the shear modulus is assumed to be a linear function while the Poisson's ratio is assumed to be a constant. In the contact area, it is assumed that the friction is one of Coulomb type. With this model the fundamental solutions for concentrated forces acting perpendicular and parallel to the FGMs layer surface are obtained. Then the sliding frictional contact problem of a functionally graded coated half-space is investigated. The transfer matrix method and Fourier integral transform technique are employed to cast the problem to a Cauchy singular integral equation. The contact stresses and contact area are calculated for various moving stamps by solving the equations numerically. The results show that appropriate gradual variation of the shear modulus can significantly alter the stresses in the contact zone.     

摩擦接触问题的比例边界等几何B可微方程组方法

摩擦接触问题是计算力学领域最具挑战性的问题之一,接触系统的泛函具有非线性、非光滑的特点,导致接触算法的收敛性与精确性难以保证.因此将比例边界等几何分析(scaled boundary isogeometric analysis,SBIGA)与B可微方程组(B dierential equation,BDE)相结合,提出了求解二维摩擦接触问题的比例边界等几何B可微方程组方法.在比例边界等几何坐标变换的基础上,通过虚功原理推导了关于边界控制点变量的接触平衡方程,表示成B可微方程组形式的接触条件可被严格满足,求解B可微方程组的算法的收敛性有理论保证.此比例边界等几何B可微方程组方法(SBIGA-BDE)只需在接触体边界进行等几何离散,使问题降低一维,能精确描述接触边界,并可通过节点插入算法进行真实接触区域的识别.此外,由于几何建模和数值分析使用相同的基函数,节约了划分网格的时间.以赫兹接触问题和悬臂梁摩擦接触问题为例,通过与解析解及数值计算软件ANSYS计算结果进行对比,验证了该方法求解二维摩擦接触问题的有效性及高精度等特点.      

An implicit method for frictional contact, impact and rolling

In this paper an implicit method for frictional contact, impact and rolling is suggested. A nonclassical formulation of a two-dimensional hyperelastic body unilaterally constrained to rigid supports is proposed by following the ideas of Moreau and Jean. A total Lagrangian formulation of the system is given. The elastic properties are defined by coupling the second Piola–Kirchhoff stress to the Green–Lagrange strain via the Kirchhoff–St. Venant law. The equation of motion is written in the spirit of Moreau by using the mean value impulses introduced by Jean. The mean value impulses appear explicitly in the equation of motion. In such manner the treatment of nonconstant kinematic transformation matrices becomes straightforward. The rigid supports are described by smooth functions. By utilizing these functions and the mean value impulses, new contact/impact laws of Signorini and Coulomb type are formulated. The governing equations are solved by a nonsmooth Newton method. This is performed by following the augmented Lagrangian approach and deriving the consistent stiffness matrix as well as the contact stiffness matrices. Three two-dimensional examples are solved by the method: a contact problem, an impact problem and a rolling contact problem.     

Subsonic semi-infinite crack with a finite friction zone in a bimaterial

Propagation of a semi-infinite crack along the interface between an elastic half-plane and a rigid half-plane is analyzed. The crack advances at constant subsonic speed. It is assumed that, ahead of the crack, there is a finite segment where the conditions of Coulomb friction law are satisfied. The contact zone of unknown a priori length propagates with the same speed as the crack. The problem reduces to a vector Riemann–Hilbert problem with a piece-wise constant matrix coefficient discontinuous at three points, 0, 1, and ∞. The problem is solved exactly in terms of Kummer's solutions of the associated hypergeometric differential equation. Numerical results are reported for the length of the contact friction zone, the stress singularity factor, the normal displacement u₂, and the dynamic energy release rate G. It is found that in the case of frictionless contact for both the sub-Rayleigh and super-Rayleigh regimes, G is positive and the stress intensity factor K_II does not vanish. In the sub-Rayleigh case, the normal displacement is positive everywhere in the opening zone. In the super-Rayleigh regime, there is a small neighborhood of the ending point of the open zone where the normal displacement is negative.     

解三维摩擦接触问题的一个二阶锥线性互补法

针对三维摩擦接触问题的求解,给出了一种基于参变量变分原理的二阶锥线性互补法. 首先,基于三维Coulomb摩擦锥在数学表述上属于二阶锥的事实,利用二阶锥规划对偶理论,建立了三维Coulomb摩擦接触条件的参变量二阶锥线性互补模型,它是二维Coulomb摩擦接触条件参变量线性互补模型在三维情形下的自然推广;随后,利用参变量变分原理与有限元方法,建立了求解三维摩擦接触问题的二阶锥线性互补法. 较之于将三维Coulomb摩擦锥进行显式线性化的线性互补法,该方法无需对三维Coulomb摩擦锥进行线性化,因而在保证精度的前提下所解问题的规模要小很多. 最后通过算例展示了该方法的特点.      

New models for a class of adhesive grippers. The hemivariational inequality approach

Summary A hemivariational inequality model for adhesive grasping problems is proposed and studied in this paper. The unilateral frictionless and frictional contact effects between the fingertips and the grasping object that lead to linear complementarity problems with singular matrices for the study of static equilibrium of the gripper-object system are generalized here to cover adhesive multifingered grippers. Adhesive effects are modelled by appropriately defined, generally nonconvex, yield sets in the space of contact stresses, friction stresses, gaps or frictional slips and their combinations. The hemivariational inequality problem that arises may involve copositive plus, symmetric matrices and nonempty closed sets for the frictionless gripper problem and copositive plus, nonsymmetric matrices with starshaped sets for the frictional case. Solvability conditions that guarantee the existence of a solution to the gripper problem are given. They specify the conditions which are required to hold between external forces, fingertip mechanical behavior and finger placement in order to solve the gripper problem.