轴线平行圆柱体滚动接触时的轴向移动阻力

应用摩擦学中的预位移原理，求解轴线平行的两粗糙圆柱体滚动接触时的轴向移动阻力，分析了轴移速比、法向接触载荷、圆柱体表面粗糙层参数对轴向移动阻力的影响，并给出了计算轴向移动阻力的简化公式．     

基于Hertz 理论圆柱和平面之间的滑动接触分析

当圆柱与平面之间产生摩擦滑动接触时,基于Hertz理论,推导了平面内应力分量的计算公式,分析了应力的分布,描述了应力分布的特点.在接触表面上,存在最大拉应力和最大压应力;最大拉应力出现在拖动边的边缘,而最大压应力发生在导向边的区域.因此,如果接触平面因拉或压应力引起失效,首先会出现在接触表面;第一型裂纹的产生及扩展首先会发生在拖动边的边缘.剪应力,主剪应力及Von-M ises等效应力(当摩擦系数较小时)的最大值均会出现在接触体内,存在于导向边.因此,塑性滑移开始首先会发生在接触体内,然后,才扩展到接触表面.这些结论为工程应用研究提供了依据.     

非稳态纯滚动接触弹塑性分析

建立了二维弹塑性非稳态循环纯滚动接触有限元模型.材料本构采用一种较好的循环塑性模型,并通过材料用户子程序在通用有限元软件ABAQUS中自定义该本构模型.通过在弹塑性无限半空间表面上重复移动随时间按简谐规律变化的赫兹法向载荷来模拟非稳态循环纯滚动接触过程.通过数值模拟,得到接触表面附近的残余累积变形、应变和残余应力.不同的最大赫兹接触压力对残余应力和残余应变影响较大.在简谐变化的法向接触载荷作用下接触表面的变形呈波浪形,随着滚动次数的增加,该波状表面沿载荷移动相反方向逐渐移动,但移动速率要衰减.波状表面波谷处的残余应力、应变和变形大于波峰处.随滚动次数的增加,残余应力增大但很快趋于稳定,残余应变也增大但增大速率衰减.     

考虑摩擦因素的结合面切向接触刚度分形预估模型及其仿真分析

为建立更完善和精确的结合面接触刚度模型,本文根据分形理论和摩擦学原理,从微观角度建立了考虑摩擦因素的结合面切向接触刚度分形预估模型.通过数值仿真分析研究了接触载荷、分形维数、摩擦系数和接触面积等因素对结合面切向接触刚度的影响.分析结果表明:结合面切向接触刚度随法向载荷和分形维数的增加而增大,而随分形尺度参数的增大而减小;摩擦系数对结合面切向接触刚度的影响较大,不同实际接触面积下的切向刚度相差较大;当分形维数较小时,摩擦系数对结合面切向刚度的影响将降低.这些研究对于进一步开展结合面的动力学特性研究具有重要意义.     

锐缘硬质轮在金属表面的滚动与滑动

在工程应用的很多摩擦对偶中，往往有一方是锐形的，因此，研究锐形滑块的滚动与滑动具有重要的实用意义。以锐缘硬质轮在挤压铝和软钢两种平面试样表面的滚动和滑动为研究对象，考察了两种摩擦过程中材料的受力状况和变形状况，结果表明，要在同种金属试样表面造成相同宽度的槽痕，滚动法所需要的垂直负荷比滑动法所需要的大，然而滚动法所需要的切向牵引力却比滑动法所需要的小．对两种摩擦过程中界面上的应力与材料的强度性能之关系的分析讨论指出，滚动试验中接触界面上的正应力相应于变形材料的屈服应力，因而滚动试验可以作为材料硬度的一种测量方法；滑动试验中接触界面上的剪应力相应于变形材料的流动剪应力；但滚动试验中界面剪应力与材料的流动剪应力之间，以及滑动试验中界面正应力与材料的屈服应力之间却都没有直线关系．     

钢轨短波长波浪形磨损的安定性分析

针对轮-轨滚动接触的短波长波浪形磨损现象,采用有限元法分析了三维实体模型的接触状态,通过计算分析了高频力作用下接触表面的塑性变形过程.结果表明:在一定的运动条件下,由于重复滚压作用,接触表面发生硬化并达到安定极限状态,生成有规律的短波长变形;钢轨表面塑性变形受枕木间距的影响;就具有随动硬化特性的钢轨材料而言,当摩擦系数μ<0.3时,屈服现象发生在材料表层下方;随着摩擦系数的增大,接触表面的切向力增大,安定极限的临界接触压力Po降低,屈服点移向接触表面,材料失效加快.     

风沙环境下钢结构涂层的冲蚀磨损力学性能研究

针对钢结构涂层耐久性受风沙侵蚀和劣化问题,在研究风沙环境特征和钢结构涂层力学性能的基础上,利用接触理论和LS-DYNA有限元分析程序分析了涂层受冲蚀时的最大接触动力、最大接触半径、最大接触应力,并分析了应力场分布规律和屈服极限风沙流速度。分析结果表明：最大接触动力、最大接触半径、最大法向动应力随冲蚀速度和角度的增加而增大,当冲蚀速度、角度分别为35m/s、90°时,其达到最大值4.09×10-2N、4.24×10-5m、1.08×107Pa;最大切向动应力随冲蚀角度的增加而减小,当冲蚀速度为35m/s 时,冲蚀角度为5°和90°时对应的应力分别为9.36×107Pa 和0;得到了涂层表面和内部应力场,在对称轴0.56z a=处,剪应力最大值0.277 pτ=1max 0,剪应力最大值点即是材料接触流动的起始点,因此可预计塑性流动将在涂层表面下方开始发生;当垂直冲蚀时涂层材料不屈服的极限射流速度为18m/s。     

作用在岩体裂隙网络中的渗透力分析

从岩体裂隙网络渗流的特点出发, 以单裂隙渗透力分析为基础, 分析了二维及三维情况下岩体裂隙网络渗流对岩体裂隙壁施加的两种作用力: 垂直于裂隙壁使裂隙产生扩容的法向渗透静水压力以及平行于裂隙壁和裂隙水流方向一致的切向动水压力, 推导出二维及三维情况下裂隙单元因这两种作用力而产生的等效结点力, 并应用算例定量分析了岩体裂隙网络渗透静水压力和动水压力共同作用对岩体应力的影响, 结果显示: (a)渗透力作用下裂隙上部岩体压应力减小, 而裂隙下部岩体压应力增大, 最大压应力增大 10 .5 3%; (b)渗透力作用下裂隙岩体拉应力增大, 最大拉应力增大 9.0 9%; (c)裂隙渗透力使岩体剪应力增大, 最大值达 2 3.75 %。     

干气密封滑动摩擦界面切向接触刚度分形模型

为揭示干气密封滑动摩擦界面高频微幅自激摩擦振动规律，用分形参数表征摩擦界面形貌特性，根据重新建立的微凸体接触变形方式，以及对非协调弹性体在切向力作用下初始滑动问题的研究，建立了干气密封滑动摩擦界面切向接触刚度分形模型. 通过数值对切向接触刚度的影响因素进行了分析，研究结果表明:切向接触刚度随分形维数、真实接触面积和材料特性系数的增大而增大；切向接触刚度随特征尺度、摩擦系数的增大逐渐减小. 相比于分形维数、特征尺度和材料特性系数对切向接触刚度的影响，摩擦系数的影响相对较小. 这些研究结果为进一步研究干气密封高频微幅自激摩擦振动奠定了基础.      

钢丝绳捻制成形接触问题的有限元分析

基于虚功原理建立了钢丝绳捻制成形过程中考虑接触摩擦时的刚度方程,将小球分裂算法用于接触搜索,Augmented Lagrange法用于计算接触力．以钢丝绳一次捻制成形过程为例,分析了摩擦系数和自扭转系数对接触应力和加工应力域变的影响．研究结果表明。摩擦系数对剪应力的影响比对等效应力和等效塑性应变的影响大;考虑摩擦时钢线表面的轴向残余应力计算结果比无摩擦条件下的计算结果更接近实验结果．     

2D FE–DEM analysis of contact stress and tractive performance of a tire driven on dry sand

Normal and tangential stresses acting over a contact interface of a tire driven on dry sand were investigated to expand the applicability of our model incorporating 2D FE–DEM with proportional–integral–derivative (PID) control. A simple averaging method for contact reaction was introduced: computational segments were defined over the lower half part of the tire circumference that translates without rotation with the tire; then the contact stresses were calculated segment by segment. For the analysis, it was assumed that the tire was in rigid contact mode and that it would travel on the model sand terrain in stationary condition. The integration of normal and tangential contact stresses with respect to the angle of rotation was then applied to calculate the vertical contact load, gross tractive effort, net traction, and running resistance of the tire by parametric (or semi-empirical) analysis. The result of tractive performance obtained through the parametric analysis was found to be similar to the result of tractive performance obtained directly using FE–DEM analysis. A forward shift of the consistent angle of rotation for maximum normal contact stress and that for maximum tangential contact stress with the increase of slip from 22% was also observed in the FE–DEM result.     

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

利用断层扫描数据,图像处理和反求技术建立了骨骼、软组织以及假肢接受腔的三维有限元模型,分别施加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时,残肢皮肤均处于黏着状态;皮肤弹性变形越大,黏着程度越大.     

A photoelastic study of contact between a cylinder and a half-space

Three-dimensional photoelasticity was employed to study a cylinder in contact with a half-space. Both bodies were modeled in epoxy resin. Three loading cases were examined, namely, the cylinder lying on its side subject to a load normal to the plane, the cylinder on its side subject to both normal and tangential loads and the cylinder standing on its end and subject to a normal compressive load, i.e., as a circular punch. The cylinders and the half-space, which was represented by a large block, were stress frozen with a known coefficient of friction and using relatively small loads so that the strain levels were low. After slicing the cylinders, which resulted in lower fringe orders than could be readily analyzed manually, an automated system based on phase stepping was used to record and process the data. Distributions of maximum shear stress and Cartesian shear stress were obtained for a large area of the slice. Stress separation was performed, using the shear difference method, to obtain the Cartesian stress components in the plane of symmetry of the half-space. These results provide confirmation, by experiment, of the theoretical and numerical models of this type of contact obtained by other investigators.     

A frictional receding contact plane problem between a functionally graded layer and a homogeneous substrate

This paper investigates the plane problem of a frictional receding contact formed between an elastic functionally graded layer and a homogeneous half space, when they are pressed against each other. The graded layer is assumed to be an isotropic nonhomogeneous medium with an exponentially varying shear modulus and a constant Poisson’s ratio. A segment of the top surface of the graded layer is subject to both normal and tangential traction while rest of the surface is devoid of traction. The entire contact zone thus formed between the layer and the homogeneous medium can transmit both normal and tangential traction. It is assumed that the contact region is under sliding contact conditions with the Coulomb’s law used to relate the tangential traction to the normal component. Employing Fourier integral transforms and applying the necessary boundary conditions, the plane elasticity equations are reduced to a singular integral equation in which the unknowns are the contact pressure and the receding contact lengths. Ensuring mechanical equilibrium is an indispensable requirement warranted by the physics of the problem and therefore the global force and moment equilibrium conditions for the layer are supplemented to solve the problem. The Gauss–Chebyshev quadrature-collocation method is adopted to convert the singular integral equation to a set of overdetermined algebraic equations. This system is solved using a least squares method coupled with a novel iterative procedure to ensure that the force and moment equilibrium conditions are satisfied simultaneously. The main objective of this paper is to study the effect of friction coefficient and nonhomogeneity factor on the contact pressure distribution and the size of the contact region.     

Fast evaluation of friction forces in traction-drive contacts including thermal effects

Summary  In automotive traction drives, power is transmitted by friction forces. The friction forces result from the shear stresses developed in lubricated and highly loaded contacts between rolling bodies. Due to the kinematics of a traction drive, shear velocities occur in both the rolling direction and perpendicular to it. Due to these shear velocities and by normal pressure, the lubricant is forced to build up shear stresses. The increase of the shear stresses may be modelled by a nonlinear viscous element. The describing differential equations are coupled by the equivalent shear stress, which defines the nonlinear behaviour of the element. A fast method is described to evaluate the coupled differential equations. By using a known analytical approximation for the equivalent shear stress, the differential equations are decoupled and can be solved analytically. In an iterative procedure the equivalent shear stress is updated, and the complete solution is found. The iterative method is extended to account for thermal effects in the contact. Received 17 June 1999; accepted for publication 26 October 1999     

Contact forces in anisotropic frictional granular materials

The probability density functions (PDFs) of contact forces in anisotropic, cohesionless and frictional granular materials are studied numerically and theoretically. Using discrete element simulations of biaxial deformation of a large two-dimensional assembly consisting of 200,000 disks, it is observed that the PDFs for the normal and tangential components of the contact forces depend significantly on contact orientation. The PDFs exhibit exponential decay and the PDF for the tangential component of the contact forces is not always symmetrical with respect to zero tangential force. The shape of the PDF for the normal component of the contact forces changes with shear strain. A qualitative explanation for this change is given that is related to the biaxial deformation mechanism in which the disrupted contacts are predominantly oriented in the direction of the minor principle stress.A maximum entropy method is employed to study these PDFs theoretically, using a prescribed stress tensor as constraint. It is found that the theoretical results correspond qualitatively to many of the results obtained from the discrete element simulations. Discrepancies between theory and simulations are attributed to the fact that the kinematics have not been taken into account in the theory.     

高速列车车轮踏面剥离引起的轮轨冲击力学响应有限元模拟

车轮踏面剥离是轨道车辆车轮非圆化损伤的常见形式之一。轮轨滚动接触过程中,车轮踏面剥离会循环冲击钢轨,诱发异常大的轮轨动态相互作用,严重影响高速列车运行平稳性和安全性。基于三维轮轨滚动接触有限元模型,模拟了高速列车车轮踏面剥离引起的轮轨冲击力学响应,分析了轮轨冲击过程中的轮轨接触力/压力、接触斑及黏/滑特性、钢轨表面节点速度分布和应力/应变状态等响应特征,讨论了列车速度、剥离长度和剥离深度等关键参数对轮轨冲击响应的影响。结果发现,车轮踏面剥离引起的轮轨动态垂向接触力随列车速度的提高呈现出先增大后减小的变化趋势,并在列车速度为300 km/h出现最大值,约为轮轨准静态垂向接触力的1.35倍;随着剥离长度的增大,轮轨动态接触力、轮/轨von Mises应力和等效塑性应变均显著增大;随着剥离深度的增大,仅车轮von Mises应力和等效塑性应变显著增大。     

基于Iwan模型的接合面切向响应建模

针对工程中常见预紧力作用下的搭接接头,研究其在小幅切向位移激励时的切向位移响应问题,为此提出一种新的基于实际表面形貌和材料性能参数的滑移力密度分布函数.应用该分布函数得到搭接接头切向响应本构模型,并获得单位加载周期内的迟滞曲线和能量耗散值, 通过与已出版的实验结果相对比,发现得到的模拟值与实验结果吻合, 证明该模型的合理性.在此基础上利用该分布函数研究了接合面切向位移与切向力、切向接触刚度及能量耗散之间的关系,结果表明: 建立的模型能很好地描述接合面间切向力与切向位移之间的关系,临界滑移力函数开始迅速上升, 到达最大值后迅速收敛到零;切线力与切向位移之间表现出非线性特性, 随着切向位移的增大,切向接触刚度表现出"软化"现象;初始切向刚度与法向载荷、粗糙度参数及塑性指数有关, 对于确定的接触表面,法向力越大, 初始切向刚度越大; 初始切向刚度同样也随着塑性指数的增大而增大.      

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

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

Multiaxial fatigue life prediction of kiln roller under axis line deflection

This paper investigates the multiaxial fatigue life of the roller in rolling contact with wheels with respect to axis line deflection. The multiaxial fatigue criteria proposed by Wang and Brown, together with the ralnflow counting method and Miner- Palmgren＇s rule, are applied to the cumulative damage estimation and life prediction. As the axis line deflection of overlong kilns generally results in asymmetric load distribution on each roller, the load ratio is introduced to describe the deflection for quantitative stress analyses. The stress analyses are performed within the finite element code ANSYS. The tangential friction stress is calculated in terms ：of the condition of the rolling contact area. By taking one roller as an example, the plotted fatigue life versus load ratio curve discovers how the axis line deflection affects the fatigue life. This study is significant to prevent the fatigue failure of the roller and can provide basis to adjust and optimize the axis line of the rotary kiln.