原子尺度摩擦行为的分子动力学模拟

应用大规模分子动力学方法,模拟了具有原子级光滑和原子级粗糙形貌的刚性球形探头与弹性平面基体的干摩擦行为,研究了无/有粘附条件下的载荷与摩擦力、载荷与真实接触面积,以及摩擦力与真实接触面积之间的关系,对纳米尺度下的摩擦行为规律进行了分析。几种系统的真实接触面积－载荷关系都与相应的连续力学接触模型定性的一致,它们分别是Hertz光滑表面接触模型、Greenwood-Williamson粗糙表面接触模型和Maugis-Dugdale粘着接触模型。无论是由光滑表面还是粗糙表面构成的摩擦系统,在无粘附条件下摩擦力与载荷成正比,而摩擦力与真实接触面积之间没有一个简单的关系;在粘附条件下摩擦力与真实接触面积成正比,而摩擦力与载荷之间表现为Maugis-Dugdale模型预测的亚线性关系。我们的研究表明,当表面作用从无粘附到粘附时,控制摩擦力的决定因素从载荷转变为接触面积,摩擦行为从载荷控制摩擦转变为粘着控制摩擦。     

考虑弹塑性变形阶段的干气密封接触模型

为探究干气密封摩擦界面在变形全阶段的接触特性,基于分形接触理论及微观接触力学理论,充分考虑微凸体变形的3个阶段,通过余弦函数构建干气密封全阶段接触模型,并分别与GW模型、KE模型和ZMC模型三种经典接触模型及相关文献的试验数据作对比,验证本文中接触模型的合理性与正确性. 最后对干气密封摩擦界面接触特性的主要影响因素进行探究. 模型研究结果表明:干气密封摩擦界面的接触特性与分形维数、特征尺度及两表面的真实接触面积有关. 接触特性与分形维数和两表面真实接触面积呈正相关,与特征尺度呈负相关. 分形维数越大,接触载荷与接触刚度的数量级就越大,且接触载荷变化范围相对较大. 当特征尺度每次以1个数量级递增时,接触载荷与接触刚度的变化范围较小,都在1个数量级内.      

基于基底长度的粗糙表面分形接触模型的构建与分析

为建立更为精确的粗糙表面接触模型,根据微凸体变形特征、分形理论以及摩擦的作用,从微观角度基于基底长度建立了粗糙表面分形接触模型.通过与其他粗糙表面接触模型和实验数据的比较,验证了本文模型的正确与合理,并由数值仿真分析了分形维数、接触载荷与真实接触面积之间的相互关系.分析结果表明:特征尺度一定时,要维持一定的真实接触面积,分形维数越大,所需要的力也越大;分形维数与特征尺度一定时,随着载荷的增加,接触面积也在增加;特征尺度与接触力一定时,随着分形维数的增大,真实接触面积在减小.该模型的建立为进一步研究粗糙表面的摩擦、磨损与润滑提供了理论依据.     

摩擦过程中真实接触面积的灰色分析

灰色系统理论是近几年以随机过程为研究对象而形成和发展起来的，在农业和自动控制等方面的研究与应用发展较快．为了使这种理论应用于摩擦学研究领域，从光全反射测量和计算机数值计算两方面，考察了粗糙表面接触时真实接触面积的随机特性，指出摩擦过程中的真实接触面积是一种含有随机成分的变化量，研究目标应当侧重于真实接触面积的变化趋势和范围，而不是它的精确解．在此基础上，首次利用灰色系统理论，分别对静态和动态接触条件下真实接触面积的变化规律进行了灰色分析，建立了精度比较高的预测真实接触面积的灰色模型ＣＭ（１，１）．这种模型对动、静接触面积的计算值与实测值具有良好的一致性，可以用于动态跟踪描述摩擦过程中真实接触面积的变化规律     

基于分形理论的双渐开线齿轮接触应力研究

综合考虑接触面粗糙度、材料特性等因素对齿轮接触应力的影响，基于分形理论和经典Hertz接触理论建立双渐开线齿轮分形接触模型. 该模型中，影响载荷和实际接触面积的主要因素包括分形维数、粗糙度幅值和材料特性参数. 理论分析表明:分形维数一定时，真实接触面积随着载荷的增大而增大；载荷一定时，接触面积随着粗糙度幅值的增大而减小；随着材料特性参数值的增加，在一定程度上加强了软材料轮齿承载能力，同时会使得微凸体由弹性变形到塑性变形的临界面积减小. 对比分形接触模型和有限元模型两种计算双渐开线齿轮轮齿接触应力方法，结果证明了分形接触模型计算双渐开线齿轮接触应力的有效性.      

分形粗糙表面弹塑性接触力学模型

为了建立更为精准的粗糙表面接触模型,基于分形理论提出一种修正的弹塑性接触模型,该模型同时考虑了微凸体的弹性、弹塑性、完全塑性三种变形状态;建立了临界第一弹塑性接触面积、临界第二弹塑性接触面积与临界弹性接触面积之间的关系,分别推导出三种变形阶段的接触载荷与接触面积之间的关系式;结合海洋岛屿面积分布密度函数,获得了接触表面上总接触载荷与真实接触面积之间的表达式。计算结果表明:单个微凸体的临界弹性接触面积不但与材料属性和分形参数有关,也与微凸体的尺寸有关;分形参数一定时,随着接触载荷的增大,真实接触面积也在增大,弹性接触面积与真实接触面积的比值逐渐减小。该模型的建立为进一步研究粗糙表面的摩擦、磨损、润滑提供了理论依据。     

气动加载与振动激励耦合试验方法研究

飞机结构在飞行过程中同时承受气动载荷和振动载荷的联合作用,这两种载荷的耦合加载试验对于飞机结构成为一项重要的研究内容,所以有必要对此类试验的可行性及其耦合加载方式进行研究。此次试验以气囊加载静载/常规疲劳载荷状态下试件的振动响应测试为目的,设计符合试验要求的试件和整套试验装置。得到了气囊5种不同加载情况下试件振动响应变化情况,并对此试验结果进行了理论分析,得出以下结论:a)气囊模拟静载/常规疲劳载荷加载不会大幅改变结构本身振动特性,此耦合试验方法所模拟环境比较接近飞机结构真实载荷环境;b)加载气囊的个数、部位及加载力的不同对试件结构的振动响应有一定影响,应增加气囊蓄能器或在试验前进行分析以选择合理的加载点。     

粗糙表面接触分形模型的提出与发展

自表面接触塑性变形模型问世以来，经过近４０年的发展，已经形成以分形几何理论为基础的Ｍ－Ｂ粗糙表面接触分形模型．Ｍ－Ｂ模型以分形参数代替统计学参数表征粗糙表面，推导出了实际接触面积与载荷的关系，以及实际弹性接触面积和实际塑性接触面积的计算公式，指出了影响接触面变形性质的因素与规律．由于分形参数的尺度独立性，可望利用Ｍ－Ｂ模型对接触面积的预测不受测量仪器分辨率和取样长度等因素的影响，故其比基于统计分析的Ｇ－Ｗ接触模型更为合理．尽管如此，Ｍ－Ｂ模型还有待完善，多方面的问题尚待进行深入研究     

基于共轭梯度法和快速傅立叶变换的粗糙表面微动接触数值研究

实际工程表面多为粗糙表面,研究粗糙表面的表面形貌对微动接触中压力和应力的影响具有重要意义。本文研究接触过程中法向载荷保持不变,切向载荷为周期性的交变载荷。首先,建立接触算法和模型,算法核心是利用共轭梯度法CGM(Conjugate gradient method)计算微动接触中的表面压力及切向应力并利用快速傅里叶变换FFT(Fast Fourier transform)加快计算速度。然后,对单峰表面、正弦表面和随机粗糙表面的接触进行数值研究。结果表明,表面幅值对切向载荷-位移曲线以及接触过程中的能量耗散有影响,表面幅值越大,相同切向载荷作用下产生的切向位移越大,能量耗散也越大。     

计及摩擦的球面弹塑性分形模型与应用

为了研究球面接触微凸体的弹塑性变化规律,在M-B分形模型的基础上考虑基底直径的实际变化情况,引入摩擦力修正因子,并基于修正的球面接触基底直径密度分布函数,最终建立了考虑摩擦因素的球面弹塑性分形模型;并以两球体接触为研究对象,分析了模型中的相关变量对接触强度的影响规律。仿真分析结果表明:当载荷不变时,分形维数的增加会导致真实接触面积先保持不变后增大再减小;当分形维数在1.35～1.7区间时,接触面积上升的速度最快;增大摩擦系数或降低材料特性参数,有利于提高接触承载能力;粗糙度幅值参数对真实接触面积和弹性接触面积占比的影响大略一致;当粗糙度幅值增加时,真实接触面积和弹性接触面积占比都随之降低。     

应用光全反射法及图像处理技术研究真实接触面积

光全反射法是测量真实接触面积的常用手段,然而用其测得的接触图像小,接触斑点的边界模糊。作者应用波动光学的有关理论对光全反射法的基本原理进行了深入的研究,进而提出用图像处理技术识别接触区域的边界。通过钢球与平面接触的试验研究,验证了这种测量方法的准确性和可靠性,并且根据试验结果分析了赫兹接触理论的局限性。文章最后还阐述了对典型工程表面接触面积的测量与结果。     

Hysteretic linear and nonlinear acoustic responses from pressed interfaces

An analysis of the linear and nonlinear acoustic responses from an interface between rough surfaces in elastoplastic contact is presented as a model of the ultrasonic wave interactions with imperfect interfaces and closed cracks. A micromechanical elastoplastic contact model predicts the linear and second order interfacial stiffness from the topographic and mechanical properties of the contacting surfaces during a loading–unloading cycle. The effects of those surface properties on the linear and nonlinear reflection/transmission of elastic longitudinal waves are shown. The second order harmonic amplitudes of reflected/transmitted waves decrease by more than an order of magnitude during the transition from the elastic contact mode to the elastoplastic contact mode. It is observed that under specific loading histories the interface between smooth surfaces generates higher elastoplastic hysteresis in the interfacial stiffness and the acoustic nonlinearity than interfaces between rough surfaces. The results show that when plastic flow in the contacting asperities is significant, the acoustic nonlinearity is insensitive to the asperity peak distribution. A comparison with existing experimental data for the acoustic nonlinearity in the transmitted waves is also given with a discussion on its contact mechanical implication.     

Ultrasonic assessment of rough surface contact between solids from elastoplastic loading-unloading hysteresis cycle

An ultrasonic method to characterize the elastoplastic contact between two rough surfaces is presented. Ultrasonic experiments are performed on three different interfaces formed by aluminum surfaces with different levels of roughness. The frequency-dependent ultrasonic reflection coefficient from the interface is measured during loading and unloading cycles as a function of pressure, from which the ultrasonic interfacial contact stiffness is reconstructed by the least-squares inversion procedure. It is shown that one should distinguish between the ultrasonic (dynamic) interfacial stiffness and static interfacial stiffness for rough surfaces in elastoplastic contact (they are identical for purely elastic contact). It is shown that ultrasonic stiffness is associated with local unloading stiffness. An elastoplastic micromechanical model is used to describe the plasticity-induced hysteresis in the ultrasonically measured interfacial stiffness during loading-unloading cycles. The topographic parameters of the interface contact are reconstructed by matching the model-predicted results with the experimentally determined ultrasonic stiffness. Using these parameters the real area of contact, which is not directly measurable, is predicted during loading-unloading cycles using the model.     

Axi-symmetric bodies of equal material in contact under torsion or shift

Summary Two axi-symmetric bodies are pressed together, so that their axes of symmetry coincide with the contact normal and the normal force is held constant. A small torque about the contact normal or a small tangential force is applied. For bodies of equal material, the normal and tangential stress states are uncoupled, and can solved separately. The surfaces of the bodies are thought as a superposition of infinitesimal rigid flat-ended punches. Consequently, the normal stress distribution can be calculated as a summation of differential flat punch solutions. A formula results, which is identical with the solution of Green and Collins. After application of a torque an annular sliding area forms at the border of the contact area. For reasons of symmetry, the common displacement of the inner stick area must be a rigid body rotation. Similarly to the normal problem, the solution can be thought as a superposition of rigid punch rotations. The tangential solution can be derived analogically, in form of a superposition of rigid punch displacements. The present method also solves the problem of simultanous normal and torsional or tangential loading with complete adhesion. As an example, Steuermann's problem for polynomial surfaces of the formA _2nr²ⁿis solved. The solutions for constant normal forces can be used as basic functions for loading histories with varying normal and tangential forces.     

用光电技术测定爆轰参数

本文介绍了一种用光电技术同时测定爆速和爆温、爆压和爆温的方法。用光电比色技术测定炸药爆轰温度的同时,用两个爆轰面上爆轰的时间差计算爆轰波速度;或用透明介质中冲击波速度来反算炸药中的爆轰压力。方法原理可靠,技术简便。     

A novel design scheme for acoustic cloaking of complex shape based on region partitioning and multi-origin coordinate transformation

Acoustic cloaking is an important application of acoustic metamaterials. This article proposes a novel design scheme for acoustic cloaking based on the region partitioning and multi-origin coordinate transformation. The cloaked region is partitioned into multiple narrow strips. For each strip, a local coordinate system is established with the local origin located at the strip center, and a coordinate transformation in the local coordinate system is conducted to squeeze the material along the strip length direction to form the cloaked region. To facilitate the implementation of the acoustic cloak, the multilayer effective medium is used to approximate the non-uniform anisotropic material parameters. The effectiveness of the proposed coordinate transformation method is verified by comparing the results from our method with those in the literature. Firstly, the results of a circular acoustic cloak in the literature are reproduced by using our finite element (FE) simulations for validation. Then, a comparison is made between the traditional coordinate transformation scheme and our new scheme for simulating an elliptical acoustic cloak. The results indicate that the proposed multi-origin coordinate transformation method has a better cloaking effect on the incident wave along the ellipse minor axis direction than the traditional method. This means that for the same object, an appropriate transformation scheme can be selected for different incident wave directions to achieve the optimal control effect. The validated scheme is further used to design an arch-shaped cloak composed of an upper semicircular area and a lower rectangular area, by combining the traditional single-centered coordinate transformation method for the semicircular area and the proposed multi-origin method for the rectangular area. The results show that the designed cloak can effectively control the wave propagation with significantly reduced acoustic pressure level. This work provides a flexible acoustic cloak design method applicable for arbitrary shapes and different wave incident directions, enriching the theory of acoustic cloaking based on coordinate transformation.     

腔体声传递函数计算方法及实验研究

阐述了低频腔体声场的声传递函数的概念，提出了一种实验求解声传递函数的方法，即在腔体壁面各测点的力与速度、腔内关键位置力与声压之间频响函数测量的基础上，求解声传递函数的线性方程组。由一矩形腔体声传递的求解实例，证明了这种声传递函数求解方法的可行性。     

Experimental Evaluation of the Interfibre Joint Strength of Papermaking Fibres in Terms of Manufacturing Parameters and in Two Different Loading Directions

The strength properties of interfibre bonds play a vital role in the response of fibrous materials, such as paper and paperboard, under mechanical loading. To help tailor the properties of such materials by chemical or mechanical treatments of the fibre wall and fibre surfaces, improved understanding of the microscopic damage and failure mechanisms of interfibre joints is desirable. In this paper, a method for manufacturing and testing of interfibre joint specimens in two principally different modes of loading is presented. The method was applied to investigate the strength of Kraft pulp interfibre joints with different geometries and in two different modes of loading: the conventional shearing mode and also a peeling mode of loading. The method was also used to investigate the influence of drying pressure, defined as the nominal pressure between two Teflon surfaces or between a rubber surface and a Teflon one as well as a simple comparison of a pulp with two different degrees of refining. The results are presented in terms of rupture force and using different methods of normalization such as nominal overlap area, length, and width of the joint region, measured using a microscope. It was shown in this study that normalising the force at rupture by either geometric parameter reduced the scatter of the strength measurements slightly, but, neither were unambiguously more successful than the other. The results of tests done with the peeling type of loading were about 20 % of those done with the conventional shearing type of loading.