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

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

Interfacial slip and creep in rolling contact incorporating a cylinder with an elastic layer

A systematic approach for investigating the interfacial behaviour of tyred systems is proposed. A two-dimensional contact model of an elastic strip, shrink-fitted onto a wheel, and subjected to different rolling contact conditions, has been adopted to illustrate the method. The model combines existing techniques to explore individual elastic contact problems and it enables us to characterise the behaviour at the strip/substrate interface caused by loads induced by a quasi-static application of stationary and moving loads on the surface of the layer. The solution is compared to the stationary load case and regimes of local slip, full stick, separation and frictional creep are identified and collated for a variety of loading conditions, materials and geometries. Further, this article presents an investigation of frictional shakedown for layered systems subjected to periodic contact loading. The term shakedown is here referred to as the possibility of developing interfacial residual stresses at the layer/substrate interface such that frictional slip, originally activated by the applied external contact load, ceases after a few loading cycles. The possible applicability of the Melan’s theorem for elastic frictional system is investigated and preliminary results presented.     

水工圆形隧洞围岩衬砌摩擦滑动接触的新解法

在实际工程中,围岩和衬砌接触时,它们之间并非完全光滑,也并非可以承受任意大的摩擦力.如果围岩与衬砌之间的剪应力大于所能承受的最大静摩擦力,接触面间将发生切向滑动,定义接触面上产生最小滑动量的状态为衬砌的真实工作状态,这种接触即为摩擦滑动接触.以库仑摩擦模型模拟围岩和衬砌之间的摩擦滑动接触,在考虑支护滞后效应的前提下,利用平面弹性复变函数方法列出了应力边界条件、应力连续条件以及位移连续条件的方程, 再结合最优化理论,建立了具有一般性的摩擦滑动接触解法.在利用混合罚函数法求解最优化问题的过程中,减少了设计变量的个数,极大地简化了优化模型,提升了优化过程的迭代速度以及优化结果的精度.以此为基础,获得了围岩和衬砌相互作用下圆形水工隧洞的应力解析解.该方法可以求解光滑接触和完全接触两种极限情况,具有一般性.同时,利用一种精确的计算方法得到了不同情况下满足完全接触条件摩擦系数的阈值,还分析了衬砌和围岩边界上切向应力的变化规律.      

双材料单伸臂叠合梁的界面滑移与自振特性分析

将叠合梁划分为接触区和分离区,接触区界面间的摩擦作用会对叠合梁的滑移、刚度和自振频率产生影响。本文给出了单伸臂叠合梁在均布荷载和集中力作用下,考虑叠合界面摩擦作用的滑移应变和滑移分布的表达式;推导了考虑叠合界面间摩擦力及摩擦力产生的抵抗弯矩共同作用下的截面刚度。假设梁按等波长和等刚度两种形式自由振动,运用传递矩阵法推导出了单伸臂叠合梁自振频率的计算公式。最后,基于一个钢-混凝土组合与叠合试验梁的算例分析了滑移沿梁长的分布及自振频率。结果表明:随界面间摩擦系数的增大,接触区界面的滑移值逐渐减小;摩擦效应对叠合梁自振频率的影响随阶次的升高越来越明显。     

A Peierls perspective on mechanisms of atomicfriction

Stick-slip behavior observed from nanoscale asperity friction experiments are often modeled by the one-degree-of-freedom Tomlinson model, which is unable to explain the effects of lattice structure and interface defects, or by molecular simulations which suffer temporal limitations in modeling the velocity- and temperature-dependent frictional behavior. A Peierls-type model developed in this work views the atomic frictional process as the initiation and gliding passage of dislocations with diffused cores on the interface. As a consequence of loss-of-ellipticity instability, the occurrence of stick-slip behavior relies on the interaction among interface slip field, contact stress fields, and existing defects. The friction stress for commensurate interface under large contact area can be approximated from the Rice model of screw-dislocation nucleation from a planar crack tip. The spatially inhomogeneous nature of rate-limiting processes and the coupling effects between contact size and interface incommensurability are successfully determined, which cannot otherwise be tackled in the Tomlinsonmodel.     

微液滴在PDMS软基体表面的动态摩擦行为研究

通过固液界面摩擦力测试装置研究了微液滴在PDMS软基体表面运动时的动态摩擦学行为,并对微液滴体积、滑动速度及软基体力学性能对固液界面动态摩擦行为的影响进行了分析. 结果表明:微液滴在软基体表面运动时表现出最大静摩擦力和动态摩擦力. 最大静摩擦力与微液滴黏度和速度梯度呈正比,动态摩擦力与微液滴体积、滑动速度和基体力学性能有关. 随着微液滴体积的增加,三相接触线长度增加,动态摩擦力增加;随着相对滑动速度增加,三相接触线长度及接触角滞后增加,动态摩擦力增加;随着软基体弹性模量降低,固液界面黏附力增加,固液界面运动能量耗散增加,动态摩擦力增加. 研究结果可为PDMS软基体表面微液滴的精确驱动和运动参数优化提供理论指导,也可进一步丰富固液界面摩擦理论.      

垂直入射弹性波在压电体单侧接触界面上的传播特性

弹性波与压电材料接触界面的相互作用问题是工程应用中常见而复杂的问题,入射波足够强会引起界面出现滑移和分离,但滑移和分离的边界未知,边界条件具有非线性特性。通过Fourier分析,将混合边值问题的求解转化为非线性代数方程,利用软件通过迭代修正的方法进行了求解;给出3种状态边界的求解,分析入射波强度、外加应力及电场对界面状态的影响,并对高频谐波的特性进行分析,通过实例对理论推导进行验证,结果显示:入射波强度、外加荷载和电场的大小及摩擦因数均会影响到界面,通过改变这些条件可以控制界面状态,另外检测高频谐波的信号也可以反映界面状态。     

Re-polarization of elastic waves at a frictional contact interface—II. Incidence of a P or SV wave

This is Part II of a two-part paper which analyses the re-polarization of elastic waves at a frictional contact interface between two solids. The re-polarization of SH waves was solved in Part I by the use of the Fourier analysis. Here, in Part II, we consider the re-polarization of P or SV waves. It is assumed that the two solids are pressed together and, at the same time, loaded by anti-plane and in-plane shearing traction. If the incident wave is sufficiently strong, localized separation and slip may take place at the interface. As a result, the incident in-plane wave is re-polarized at the interface so that the anti-plane waves (SH waves) are induced. Using the method similar to that of Part I and considering the boundary conditions involving separation and slip, we manage to reduce the problem to a set of algebraic equations coupled with simple integral equations. An iterative method is developed based on the solution to the perfectly bonded interface. The locations and sizes of the separation and slip zones, the interface traction, the slip velocities, the global sliding velocities and the energy dissipation and partition are displayed for the case of two identical materials. It is found that the separation zones and the gaps are independent of the induced waves.     

Shakedown limits for an oscillating,elastic rolling contact with Coulomb friction

We examine experimentally and theoretically the effect of frictional shakedown of a three-dimensional elastic rolling contact. Small oscillations of the local normal forces lead to incremental sliding processes within the area of contact. Consequently, this causes a macroscopic slip motion of the two contacting bodies. If the oscillation amplitude is sufficiently small, the frictional slip ceases after the first few loading periods and a safe shakedown occurs. Otherwise the slip motion is continued and the contact fails.     

基于边界元法求解三维弹性摩擦接触问题

边界元法求解三维摩擦接触问题,其中一个关键点在于如何确定滑移方向。即当出现相对滑移时,滑移方向如何确定。当前常采用的方法是,粘结点利用切向面力得到滑移方向,滑移点利用切向相对位移得到滑移方向。不过该方法难以保证收敛性。针对这一问题,本文采用滑移方向预测技术得到滑移方向。即以后出现相对滑移时,滑移方向采用预测技术中得到的滑移方向。由于摩擦接触问题和历史加载相关,本文采用增量法求解。不同摩擦系数下的数值结果都证明了本文算法的有效性和收敛性及滑移方向预测技术的有效性。     

Use of a modified torsional Kolsky bar to study frictional slip resistance in rock-analog materials at coseismic slip rates

The knowledge of shear resistance in earth faults is of fundamental importance to our understanding of the magnitude of stress drop and the associated energy release during typical seismic rupture events. In the present study a modified torsional Kolsky bar is employed to investigate frictional slip resistance in rock-analog materials (i.e., quartz and soda lime glass) at normal stresses of relevance to earthquake physics (30–80 MPa) and co-seismic slip rates. The results indicate the coefficient of kinetic friction to be in the range of 0.2 to 0.3. These values of the coefficient of friction are much lower when compared to those obtained in rocks at quasi-static slip rates. In all experiments slip weakening is observed and is preceded by slip strengthening. The slip weakening is understood to be due to thermal weakening induced by flash-heating at asperity contacts and requires a few mm of slip to be effective; the slip strengthening is understood to be due to an increase in the real area of contact at the asperity junctions due to localized plastic flow and subsequent coalescence and solidification of local softened/melt patches at the slip interface.     

微动磨损对过盈配合结构微动参量的影响

基于Archard磨损模型与Abaqus自适应网格技术建立了仿真分析过盈配合结构微动磨损的计算模型,对配合面轮廓随循环周次的变化进行了预测,并详细研究了微动磨损对接触压应力、摩擦剪切应力及微动滑移幅值等微动参量的影响.结果表明:该计算模型能够较为准确地对配合面轮廓随循环周次的变化进行预测;由于微动磨损的作用,配合边缘处接触压应力的峰值逐渐增大,且其位置逐渐向配合内部移动;摩擦剪切应力的最大值逐渐由黏着-滑移过渡位置向磨损-未磨损过渡位置移动;张开区域的宽度以及滑移区内各位置处的微动滑移幅值均随着循环周次的增加而增大.     

A New Method for Solution of 3D Elastic-Plastic Frictional Contact Problems

ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｔａｃｔｐｒｏｂｌｅｍｉｓｏｆｐａｒｔｉｃｕｌａｒｉｎｔｅｒｅｓｔｉｎｎｕｍｅｒｏｕｓｅｎｇｉｎｅｅｒｉｎｇａｐｐｌｉｃａｔｉｏｎｓｒａｎｇｉｎｇｆｒｏｍｓｔｒｕｃｔｕｒｅ＿ｓｔｒｕｃｔｕｒｅｉｎｔｅｒａｃｔｉｏｎｔｏｍａｃｈｉｎｅｄｅｓｉｇｎａｎｄｍｅｔａｌｆｏｒｍｉｎｇｕｎｄｅｒｖａｒｉｏｕｓｌｏａｄｃａｓｅｓ.Ｃｏｎｔａｃｔｐｒｏｂｌｅｍｓａｒｅｎｏｎｌｉｎｅａｒｆｏｒｔｈｅｍｏｖｉｎｇｂｏｕｎｄａｒｙａｎｄｔｈｅｅｘｉｓｔｅｎｃｅｏｆｆｒｉｃｔｉｏｎａｌｏｎ…     

柔性梁斜碰撞问题的非线性动力学建模和实验研究

本文研究柔性梁点面斜碰撞问题。用Hertz接触模型处理法向撞击力,分别用Hertz切向接触模型和Coulomb摩擦力模型处理粘滞状态和滑动状态的摩擦力。从精确的应变与位移的关系出发,用绝对节点坐标法建立了柔性梁的动力学方程。为了准确地处理斜碰撞切向运动的复杂状态,提出滑动-粘滞切换的准则,在此基础上,设计了斜碰撞实验,数值对比了法向撞击力和法向速度的时间历程的仿真计算结果与实验结果,验证了Hertz理论在斜碰撞情况下的正确性。另一方面切向速度的实验与理论的结果对照表明滑动-粘滞切换准则的有效性。     

A modified torsional kolsky bar for investigating dynamic friction

This paper introduces an experiment to investigate dry sliding resistance of frictional interfaces at normal pressures up to 100 MPa, slip speeds up to 10 m/s and slip distances of approximately 10 mm. This new apparatus involves a novel modification of the conventional torsional Kolsky bar apparatus, employed extensively in the past for investigating high strain rate behavior of engineering materials. The new experimental configuration represents a significant improvement over conventional tribology experiments because it uses elastic torsional waves with a superimposed static compressive force to control the interfacial traction. Moreover, the apparatus allows critical frictional parameters such as the interfacial sliding resistance, slip speeds and slip without the use of transducers at the frictional interface. The usefulness of the device is demonstrated by presenting results of high-speed friction on 6061-T6 Al/1018 steel and Carpenter Hampden tool steel/7075-T6 Al tribo pairs.     

单侧摩擦接触约束夹杂物与平面波的动力相互作用--时域边界元分析

应用时域边界元法研究了瞬态平面波对单侧摩擦约束夹杂物的散射问题，假设界面摩擦遵守库仑定理，当入射波足够强时界面会出现局部分离和滑移，由于边界上的区域(分离区、滑移区和粘着区)是未知的，所以该问题实际上是个复杂的边界非线性问题，为了确定未知区域，该文发展了一种有效的迭代技术，作为算例，计算了一个无限域中圆柱埋置夹杂物对瞬态平面波的散射问题。     

Shakedown in elastic contact problems with Coulomb friction

Elastic systems with frictional interfaces subjected to periodic loading are sometimes predicted to ‘shake down’ in the sense that frictional slip ceases after the first few loading cycles. The similarities in behaviour between such systems and monolithic bodies with elastic–plastic constitutive behaviour have prompted various authors to speculate that Melan’s theorem might apply to them – i.e., that the existence of a state of residual stress sufficient to prevent further slip is a sufficient condition for the system to shake down.In this paper, we prove this result for ‘complete’ contact problems in the discrete formulation (i) for systems with no coupling between relative tangential displacements at the interface and the corresponding normal contact tractions and (ii) for certain two-dimensional problems in which the friction coefficient at each node is less than a certain critical value. We also present counter-examples for all systems that do not fall into these categories, thus giving a definitive statement of the conditions under which Melan’s theorem can be used to predict whether such a system will shake down.     

Slip weakening in rocks and analog materials at co-seismic slip rates

Determination of co-seismic slip resistance in earth faults is critical for understanding the magnitude of shear-stress reduction and hence the near-fault acceleration that can occur during earthquakes. Also, knowledge of shear resistance dependency on slip velocity, slip distance, normal stress, and surface roughness is fundamental information for understanding earthquake physics and the energy released during such events. In the present study, plate-impact pressure-shear friction experiments are employed to investigate the frictional resistance in fine-grained Arkansas novaculite rock and an analog material (soda-lime glass), under relevant interfacial conditions. The results of the experiments indicate that a wide range of frictional slip conditions can exist during a single slip event. These conditions range from initial no-slip followed by slip weakening, strengthening, and seizure at the frictional interface. For the case of glass vs. glass experiments, the first slip-weakening, with μ in the range of 0.4-0.2, is understood to be most likely due to thermal-induced flash heating and incipient melting at asperity junctions, while the slip-strengthening, with μ in the range of 0.1-0.4, is understood to be a result of coalescence and solidification of local melt patches. For the case of the fine grain Arkansas Novaculite rocks a similar range of slip conditions is observed.     

A new computational algorithm for contact friction modeling of large plastic deformation in powder compaction processes

In this paper, the large deformation frictional contact of powder forming process is modeled based on a new computational algorithm by imposing the contact constraints and modifying the contact properties of frictional slip. A simple and efficient numerical algorithm is presented for imposing the contact constraints and frictional contact properties based on the node-to-surface contact technique to simulate the large deformation contact problem in the compaction process of powder. The Coulomb friction law is used to simulate the friction between the rigid punch and the workpiece by the use of penalty approach. A double-surface cap plasticity model is employed together with the nonlinear contact friction algorithm within the framework of large FE deformation in order to predict the non-uniform relative density distribution during large deformation of powder die-pressing. Finally, the numerical schemes are examined for accuracy and efficiency in modeling of a set of powder components.