Effects of interfacial friction on flaw tolerant adhesion between two dissimilar elastic solids

Flaw tolerance refers to a state in which a pre-existing crack-like flaw does not propagate even as the material is stretched to failure near its theoretical strength. Such an optimal scenario can be achieved when the characteristic length scale is reduced to below a critical value. So far, the critical conditions to achieve flaw tolerance have been discussed mostly for homogeneous materials or for two dissimilar materials in frictionless or perfectly bonded adhesion. In this paper, we consider the role of friction in flaw tolerant adhesion between two dissimilar elastic solids. We adopt a frictional contact model in which slip is allowed wherever the shear stress along the interface reaches a threshold value defined as the friction strength. The critical length scale for flaw tolerance is derived analytically for a penny-shaped crack and for an external circular crack. Compared to the cases of frictionless contact, we find that interfacial friction can reduce the critical length scales for flaw tolerance by up to 12.5%.     

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

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

Hertzian fracture at unloading

Hertzian fracture through indentation of flat float glass specimens by steel balls has been examined experimentally. Initiation of cone cracks has been observed and failure loads together with contact and fracture radii determined at monotonically increasing load but also during unloading phases. Contact of dissimilar elastic solids under decreasing load may cause crack inception triggered by finite interface friction and accordingly the coefficient of friction was determined by two different methods. In order to make relevant predictions of experimental findings, a robust computational procedure has been developed to determine global and local field values in particular at unloading at finite friction. It was found that at continued loading it is possible to specify in advance how the contact domain divides into invariant regions of stick and slip. The maximum tensile stress was found to occur at the free surface just outside the contact contour, the relative distance depending on the different elastic compliance properties and the coefficient of friction. In contrast, at unloading invariance properties are lost and stick/slip regions proved to be severely history dependant and in particular with an opposed frictional shear stress at the contact boundary region. This causes an increase of the maximum tensile stress at the contour under progressive unloading. Predictions of loads to cause crack initiation during full cycles were made based on a critical stress fracture criterion and proved to be favourable as compared to the experimental results.     

Delamination of an elastic film from an elastic–plastic substrate during adhesive contact loading and unloading

Adhesive contact between a rigid sphere and an elastic film on an elastic–perfectly plastic substrate was examined in the context of finite element simulation results. Surface adhesion was modeled by nonlinear springs obeying a force-displacement relationship governed by the Lennard–Jones potential. A bilinear cohesive zone law with prescribed cohesive strength and work of adhesion was used to simulate crack initiation and growth at the film/substrate interface. It is shown that the unloading response consists of five sequential stages: elastic recovery, interface damage (crack) initiation, damage evolution (delamination), film elastic bending, and abrupt surface separation (jump-out), with plastic deformation in the substrate occurring only during damage initiation. Substrate plasticity produces partial closure of the cohesive zone upon full unloading (jump-out), residual tensile stresses at the front of the crack tip, and irreversible downward bending of the elastic film. Finite element simulations illustrate the effects of minimum surface separation (i.e., maximum compressive surface force), work of adhesion and cohesive strength of the film/substrate interface, substrate yield strength, and initial crack size on the evolution of the surface force, residual deflection of the elastic film, film-substrate separation (debonding), crack-tip opening displacement, and contact instabilities (jump-in and jump-out) during a full load–unload cycle. The results of this study provide insight into the interdependence of contact instabilities and interfacial damage (cracking) encountered in layered media during adhesive contact loading and unloading.     

Hertz contact at finite friction and arbitrary profiles

Axisymmetric contact at finite Coulomb friction and arbitrary profiles is examined analytically and numerically for dissimilar linear elastic solids. Invariance and generality are aimed at and an incremental procedure is developed resulting in a reduced benchmark problem corresponding to a rigid flat indentation of an elastic half-space. The reduced problem, being independent of loading and contact region, was solved by a finite element method based on a stationary contact contour and characterized by high accuracy. Subsequently, a tailored cumulative superposition procedure was developed to resolve the original problem to determine global and local field values. Save for the influence of the coefficients of friction and contraction ratio, it is shown that at partial slip the evolving relative stick-slip contour is independent of any convex and smooth contact profile at monotonic loading. For flat and conical profiles with rounded edges and apices, results are illustrated for relations between force, depth and contact contours together with surface stress distributions. The solution for dissimilar solids in a full space is transformed to a half-space problem and solved for a combination of material parameters in order to first determine interface traction distributions. Subsequently, full field values for the two solids were computed individually. In order to predict initiation of fracture and plastic flow, results are reported for the location and magnitude of maximum tensile stress and effective stress, respectively, for a range of geometrical and material parameters. In two illustrations, predicted results are compared with experimental findings related to initiation of brittle fracture and load-depth relations at nanoindentation.     

弹簧-金属丝网橡胶组合减振器迟滞力学模型及实验研究

金属丝网橡胶材料是一种完全由金属丝编织成的多孔复合材料,与传统螺旋卷制金属橡胶材料相比,其改进了成型工艺,剔除了制备过程中大量的手工工艺干扰,提高机械化程度,重合度更高,拥有更稳定的力学性能.由于金属丝网橡胶材料具有承载能力高、阻尼大、耐高温、耐低温、耐老化、抗油抗腐蚀等优良特性,在很多方面强于传统橡胶,多用于航空航天、船舶、军事武器等军工工业.弹簧$\!$-$\!$-$\!$金属丝网橡胶组合减振器具有可设计刚度和较高承载能力,但因其具有复杂的非线性迟滞特性,目前相关材料的本构模型还难以准确描述其力学特性.本文在弹簧$\!$-$\!$-$\!$金属丝网橡胶组合减振器静态迟滞力学性能实验的基础上,结合其干摩擦阻尼迟滞特性,提出了一种迟滞力学性能理论模型.根据减振器迟滞实验恢复力$\!$-$\!$-$\!$位移曲线特点,利用参数分离的方法将迟滞曲线分解为弹性恢复力和干摩擦阻尼力,分别建模求解等效刚度和干摩擦阻尼系数,以此建立了组合减振器理论模型,并与实验结果进行对比及进行误差分析,验证了理论模型的准确性.      

弹塑性微凸体侧向接触相互作用能耗

传统的结合面研究多基于光滑刚性平面与等效粗糙表面接触假设,忽略了结合面上微凸体侧向接触及相邻微凸体之间的相互作用,这导致理论模型与实际结合面存在较大出入.针对承受法向静、动态力的机械结合面,从微观上研究了微凸体侧向接触及相互作用的接触能耗.将法向静、动态力分解为法向分力和切向分力,获取弹性/弹塑性/塑性阶段考虑微凸体侧接触及相互作用的加、卸载法向分力-变形和切向分力-位移的关系.通过力的合成定理,从而获取加、卸载法向合力与总变形之间的关系,由于法向分力产生的塑性变形及切向分力产生的摩擦,导致加载、卸载法向合力-总变形曲线存在迟滞回线.通过对一个加、卸载周期内的法向合力-总变形曲线积分,获得一个周期的微凸体接触能耗,包括应变能耗及摩擦能耗.仿真分析表明:微凸体在3个阶段的能耗均随变形的增大而非线性增大.微凸体侧向接触角度越大,能耗越大,且在弹性阶段最为明显.在弹性阶段,仅存在侧向的摩擦能耗,故结合面在低载荷作用下必须采用双粗糙表面假设.在塑性阶段,由于微凸体接触能耗为应变能耗,且接触角对其能耗影响甚微,故结合面在大载荷作用下可采用单平面假设对其进行研究.相对于KE和Etsion模型,本文提出的模型与Bartier的实验结果更吻合.     

Contact stresses: a short survey of models and methods of computations

Contact stresses are identified as normal and tangential forces between contacting solids. The normal stresses are modeled using unilateral and complementary conditions, elastic response and normal compliance. Friction laws describe the tangential traction. Friction of materials depends on pressure, sliding velocity, surface temperature, time of contact, surface roughness and presence of wear debris. Phenomenological, micro-mechanical and atomic-scale models as well as non-classical models of anisotropic and heterogeneous friction are important steps in the development of friction modeling. Sophisticated friction models are desirable in vibrating systems, materials processing, rolling contacts, rubber and polymers, geomechanics, bioengineering and living systems. Main numerical methods in contact mechanics are: finite element method, boundary element method and discrete element method. To include specific contact constraints, the following computing techniques are applied: Lagrange multipliers, penalty function, perturbated and augmented Lagrangian methods, mathematical programming methods. The advances of adhesion and impact modeling are outlined in this paper.     

移动接触下求解二维闭合裂纹的双重迭代边界元法

使用子域边界元法对受移动接触弹性体作用下的二维闭合裂纹问题进行了数值计算。由于两弹性体的接触界面和裂纹表面的接触范围的大小和接触状态事先是未知的 ,对此 ,在两个接触表面同时采用迭代的方法进行了求解。在裂纹的每个裂尖上都采用了四分之一的奇异单元以保证裂尖位移场和应力场奇异性的满足。用我们编制的二维裂纹问题程序对一些中心裂纹问题进行了计算 ,计算结果与经典断裂力学的理论值比较吻合。在无摩擦的条件下 ,对一些具有不同角度且受移动接触弹性体作用下的闭合裂纹问题进行了数值计算 ,得到了一些耦合作用下的应力强度因子的计算结果     

Fretting contact of two elastic solids with graded coatings under torsion

In this paper, the fretting contact problem for two elastic solids with graded coatings is investigated. We assume a conventional axisymmetric Hertzian contact takes place between two elastic solids under the action of the normal pressure. The application of the torque produces an annulus of slip. It is assumed that the surface shear traction within the contact area is limited by Coulomb’s friction law and the torsion angel was produced within the central adhesion zone as a rigid body. The linear multi-layer model is used to model the functionally graded coating with arbitrarily varying shear modulus. This model divides the coating into a series of sub-layers with the elastic modulus varying linearly in each sub-layer and continuous on the sub-interfaces. By using the transfer matrix method and Hankel integral transform technique, this problem is formulated as the solution of the Cauchy singular integral equations. The contact tractions are calculated by solving the equations numerically. The results show that the appropriate gradual variation of the shear modulus can significantly alter the contact tractions. Therefore, graded coatings may have potential applications in improving the resistance to fretting contact damage at the contact surfaces.     

Detachment of a rigid solid from an elastic wavy surface: Theory

The mechanics of detachment between a wavy elastic half space and a rigid solid is considered. Solutions for the axisymmetric problem of a rigid sphere and the plane strain problem of a rigid cylinder detaching from a wavy surface are developed. The interacting solids are taken to be in complete contact over a finite area initially. It is shown that the surface waviness makes the detachment process unstable, with the interface separating in alternating stable and unstable segments. Each unstable segment dissipates mechanical energy, leading to an increase in the total work of separation compared to that of a flat surface. Further, waviness causes the maximum separation force or the pull-off force to increase during detachment, resulting in an apparent toughening of the interface. This mechanism provides an alternative explanation to the experimental observations in the literature that roughness can sometimes lead to increase in pull-off force. It also illustrates the role of roughness in the attachment capability of several insect feet possessing soft pads. The basic solution presented here can be used to analyze the detachment of surfaces with multiple scale roughness as well. The solution also suggests strategies to improve reversible adhesion of a soft material by designing optimal surface topographies.     

大变形软材料接触摩擦的在线测试装备与技术分析

我国高档密封件与液压件产品几乎全部依赖进口,国家重大技术装备及国防装备等配套零部件行业的许多问题与摩擦学有十分重要的关系,基础件摩擦学又是国际上的研究热点与发达国家的竞争高地. 在本文中阐述了国内外关键机械零部件的接触界面原位在线测试的研究进展,及其在大变形软材料密封系统中的应用,综述了摩擦润滑中迁移状态实时在线观测技术、聚合物密封界面在线观测技术和国内在摩擦过程的微区域原位研究等进展,还分析了特殊与极端工况条件下大变形聚合物密封的性能测试、可靠性寿命分析以及多工况联合测试的技术进展,分析和讨论了高性能密封件行业的基础共性难题,最后进行了总结与展望.      

Orientation-dependent adhesion strength of a rigid cylinder in non-slipping contact with a transversely isotropic half-space

Recently, Chen and Gao [Chen, S., Gao, H., 2007. Bio-inspired mechanics of reversible adhesion: orientation-dependent adhesion strength for non-slipping adhesive contact with transversely isotropic elastic materials. J. Mech. Phys. solids 55, 1001–1015] studied the problem of a rigid cylinder in non-slipping adhesive contact with a transversely isotropic solid subjected to an inclined pulling force. An implicit assumption made in their study was that the contact region remains symmetric with respect to the center of the cylinder. This assumption is, however, not self-consistent because the resulting energy release rates at two contact edges, which are supposed to be identical, actually differ from each other. Here we revisit the original problem of Chen and Gao and derive the correct solution by removing this problematic assumption. The corrected solution provides a proper insight into the concept of orientation-dependent adhesion strength in anisotropic elastic solids.     

Simulation of fatigue crack growth at plastically mismatched bi-material interfaces

The present paper reports the results of a study on fatigue crack growth in a system with an interface between two elastic–plastic solids of different yield strength. Fatigue crack growth analysis is conducted by the use of a cohesive zone model. Irreversibility of material separation processes is introduced through the load history dependent degradation of the cohesive strength. Fatigue crack growth is considered to occur along the direction perpendicular to the interface and along the interface. Crack growth rate acceleration, deceleration or arrest, as well as crack bifurcation at the interface are predicted in dependence of the plastic property mismatch of the two solids and the interface properties. The outcome of the simulations is in very good agreement with trends of published experimental data.     

Effect of friction on subsurface stresses in sliding line contact of multilayered elastic solids

A numerical integral scheme based on Fourier transformation approach is employed to investigate the effect of friction on subsurface stresses arising from the two-dimensional sliding contact of two multilayered elastic solids. The analysis incorporates bonded and unbonded interface boundary conditions between the coating layers. Two line contact problems are presented. The first one is the contact problem between a rigid cylinder and a two-layer half space and the second one is the indentation of a multilayered elastic half-space by a flat rigid punch. The effects of the surface coating on the contact pressure distribution and subsurface stress field are presented and discussed.     

小腿残肢与接受腔界面的摩擦行为分析及测试

利用断层扫描数据,图像处理和反求技术建立了骨骼、软组织以及假肢接受腔的三维有限元模型,分别施加Heel Strike、Foot Flat、Mid-Stance、Heel Off和Toe Off五个典型步态时相的载荷,计算残肢皮肤在一个步态周期内的应力、滑移量变化,确定最大临界参数;借助UMT-II多功能摩擦磨损试验机模拟残肢皮肤和接受腔摩擦界面,研究皮肤在步态周期临界参数下的摩擦行为变化.有限元结果表明:一个步态周期内五个典型时相,软组织表面最大正应力和剪切力均发生在髌韧带处,在Heel off步态时相,最大临界正应力为384.3 k Pa,剪应力为102.1 k Pa;随着正压力的减小,摩擦力减小,摩擦系数增大;当正应力下降到39.5 k Pa临界值时,接触面发生相对滑动,摩擦系数达到最大值.摩擦学试验结果表明:法向载荷分别为17 N和7 N时,残肢皮肤均处于黏着状态;皮肤弹性变形越大,黏着程度越大.     

氟橡胶O型圈低压气体密封黏滞摩擦特性实验

开展了低气体密封压力条件下的氟橡胶O型圈低压黏滞摩擦特性实验研究.对氟橡胶O型圈与不锈钢2Cr13摩擦副的摩擦力的时变曲线进行测量,重点分析了密封压力、O型圈压缩率和滑动速度对O型圈黏滞摩擦力的影响规律.结果表明:随压缩率的增加,氟橡胶O型圈摩擦力增加,并呈现出明显的回弹特征,即滞后摩擦力增加且释放时间增加;随着密封压力的增加,其所受最大静摩擦力、滑动摩擦力与滞后摩擦力均呈现非线性变化,且存在一极大值;减小往复运动速度,摩擦力数值增大,但摩擦力极值对应的密封压力值变小.     

Modelling the traction of a prototype track based on soil–rubber friction and adhesion

An experimental track layer tractor, based on an Allis Chalmers 8070 tractor (141 kW) was tested on bitumen covered concrete and on cultivated sandy loam at 7.8%; 13% and 21% soil water content. The two articulated beam-type tracks (500 mm wide × 2000 mm soil contact length) were constructed out of 500 mm long and 70 mm wide rubber covered steel track elements, carried by five steel cables (36 mm diameter). The tracks resisted inward deflection but allowed outward articulation between two smooth rear driving and two smooth front pneumatic truck tires (1060 mm diameter) per track. The contact pressure and the tangential force on an instrumented track element, as well as the total torque input to one track, were simultaneously recorded during the drawbar pull/slip tests.

Different possible pressure distribution profiles under the tracks were considered and compared to the recorded data. Two possible traction models are proposed, one constant pressure model for minimal inward track deflection, and a deformable track model with inward deflection and a higher contact pressure at both the front free-wheeling and rear driving tires. For both models, the traction force was generated mainly by rubber/soil friction and adhesion and limited soil shear. A close agreement between the measured and predicted contact pressures and traction force for individual track elements, based on the deformable track model, was observed. The recorded and calculated coefficient of traction based on the summation of the force for the series of track elements were comparable, but were considerably lower than the predicted values, probably due to internal track friction rather than soil sinkage. The tractive efficiency for both a hard or soft surface was also unacceptably low, probably caused by internal track friction.     

弹性波与单侧界裂纹相互作用问题的边界元法

措助边界元法设计了一种迭代修正方法来求解单侧界面裂纹模型与弹性波的相互作用问题,作为对算法的检验,用这种方法我们具体地分析了平面简谐弹性波对一个则界面裂纹的入射,给出了裂纹面的接触形态及应力场。     

拉扭耦合作用下柱形纤维与基底的界面黏附性能研究

受壁虎刚毛可逆黏附性能的启发,本文建立了单根弹性圆柱纤维与刚性基底黏附接触的理论和数值模型,同时考虑了拉伸和扭转载荷的耦合作用及纤维半径对界面黏附性能的影响.研究发现耦合载荷作用下柱形纤维同样存在一个临界半径,当纤维半径小于该临界尺寸时,界面应力达到均匀的理论强度分布,接触边界应力集中消失,出现缺陷不敏感现象;当纤维半径大于该临界尺寸时,界面以裂纹扩展而失效.在耦合载荷作用下纤维的临界半径小于纯拉伸而大于纯扭转时的临界尺寸,且该临界半径随着施加扭转载荷的增大而减小.表明在纯拉伸载荷下使界面黏附强度达到最优的柱形纤维,在拉伸和扭转载荷耦合作用下,由于界面失效形式的转变使界面易发生脱黏,并且界面脱黏时的拉脱力随着扭转载荷的增大而减小,理论和数值结果一致.本文结果进一步应用揭示了壁虎可以通过调控施加在其最小黏附单元上的载荷形式实现纯拉伸载荷下强黏附及耦合载荷下易脱黏的力学机制.