小Tabor数分形粗糙表面之间的黏着弹性接触分析

基于粗糙表面的分形描述和适用于小Tabor数微突体的黏着弹性接触理论,采用积分方法建立了小Tabor数分形粗糙表面之间的黏着弹性接触模型,获得弹性接触条件下的真实接触面积和载荷表达式,在此基础上采用单因素分析法分析表面粗糙度和材料性质对分离力的影响.结果表明,当分形维数增加时,粗糙表面单位面积上的微突体数目增加且高度减小,从而导致两表面间的分离力增大;当分形粗糙度参数增大时微突体高度增加,从而导致分离力减小;当材料弹性模量增加时黏着作用减弱,从而减小了分离力,而表面黏着能的增加会使分离力急剧增大.     

边界膜剪切弹性模量对于微接触性能的影响

建立接触模型,理论分析了微接触中边界膜剪切弹性模量对于接触性能的影响.接触区由两平行平面形成,属一维接触.上接触表面为粗糙表面,具有矩形微凸体.下接触表面为光滑平面.两接触表面均处理成刚性表面.微接触区中充满流体.它分成两个子区,在微接触的出口区由于极小的接触间隙充满边界膜,在微接触的入口区由于接触间隙较大充满流体膜.边界膜和流体膜行为决定整个微接触性能.当膜厚较大时,这里边界膜可看成纳米级薄膜.由于上接触表面处有限的剪应力承受能力,边界膜可于上接触表面滑移.设下接触表面处剪应力承受能力很大而边界膜在下接触表面不滑移.由于边界膜-接触表面间相互作用,边界膜黏度、密度和剪切弹性模量均沿膜厚变化,在理论分析中使用它们的等效值,这些值与边界膜厚度有关.流体膜在两个接触表面均不发生滑移,分析中不考虑流体膜剪切弹性模量.流体膜采用传统分析法.给出了理论分析和若干变工况参数下的计算结果.     

平行运动的粗糙表面润滑承载性能研究

分析了微米级润滑膜条件下,光滑表面在静止粗糙表面上平行运动时的润滑状态及其承载机制,利用Reynolds流体润滑方程分析粗糙度对油膜压力、载荷及摩擦系数的影响,采用有限差分法计算在正弦和随机粗糙峰条件下油膜的压力分布曲面图,通过改变正弦粗糙度的峰高和波长分析油膜承载能力和摩擦系数随粗糙度变化的规律,同时分析了最小油膜厚度对润滑状态的影响.结果表明:两光滑平行运动的平面无法承载,而粗糙表面微粗糙峰的收敛楔形部分可以形成流体动压润滑膜并承受一定载荷;在给定最小油膜厚度的条件下,随着正弦波峰值增加,承载能力达到最大值后缓慢降低,摩擦系数达到最小值后缓慢增大;除了粗糙峰波长很小时摩擦系数很大以外,波长对摩擦系数的影响很小,而承载能力随波长以二次曲线变化并出现最大值;在给定粗糙度幅值条件下,当最小膜厚在1～100 μm时,随着最小油膜厚度的增加,承载能力减小,摩擦系数逐渐增大.     

介观双粗糙弹塑性流体动力润滑界面法向接触刚度模型

弹性流体动力润滑状态通常出现在机械高副零部件的点/线接触部位,如齿轮、轴承和蜗轮蜗杆等. 宏观上点/线接触在介观层面表现为两粗糙表面的接触,在微观层面上则又表现为微凸体间的接触. 由于在中/重载荷作用下,粗糙表面上的微凸体发生接触后会产生弹塑性/塑性变形,从而使得两粗糙表面的弹流润滑接触转变为弹塑性流体动力润滑接触. 此外,界面的接触刚度决定了机械装备的整机刚度. 为了精确获得弹性流体动力润滑状态下界面法向接触刚度及其主要影响因素,基于界面的法向接触刚度由固体接触刚度和润滑油膜刚度两部分构成的思想,根据固体弹塑性理论和流体动力学理论,分别对界面间微凸体侧接触及部分膜流体动力润滑进行分析,从微观入手揭示双粗糙表面弹塑性流体动力润滑接触机理,进而建立考虑微凸体侧接触弹塑性变形的流体动力润滑界面法向接触刚度模型. 通过仿真分析,揭示了法向载荷、卷吸速度、表面粗糙度及润滑介质特性等因素对润滑界面法向接触刚度的影响规律. 研究表明:在相同速度、粗糙度及润滑油黏度的工况下,固体接触刚度和油膜接触刚度均随着法向接触载荷的增加呈非线性增大;在相同载荷、速度及润滑油黏度的工况下,接触表面粗糙度越大,表面形貌对于润滑状态的影响较强,固体接触刚度占界面总刚度的主要部分,界面主要由固体承载;在相同载荷、粗糙度及润滑油黏度工况下,随着卷吸速度的增大,固体接触刚度逐渐减小,油膜刚度占界面总刚度的主要部分;在相同载荷、粗糙度及速度工况下,随着润滑油黏度的增大,油膜刚度基本保持不变,固体接触刚度基本不受润滑油黏度的影响. 通过理论建模准确获得单位面积弹塑性流体动力润滑结合面法向接触刚度,对改善机械装备动态性能、提高机械装备的可靠性具有重要的理论和实际意义.      

粗糙度纹理对有限长线接触混合润滑影响

采用统一Reynolds方程建立有限长线接触混合润滑模型,研究横向、纵向和二维规则表面粗糙度的波长、幅值及工况变化对润滑影响.结果表明:波长、幅值与工况对三种表面粗糙度接触副的润滑影响类似;随着载荷增大,平均膜厚降低,摩擦系数、接触载荷比与接触面积比均增大;随着转速升高,平均膜厚增大,摩擦系数、接触载荷比与面积比均降低,其中摩擦系数随转速进一步增大而小幅升高.在润滑状态转换区域润滑特征参数变化显著,而其他润滑区域变化平缓.沿卷吸速度方向的压力与膜厚波动分布存在相位差,垂直方向则同相位;相同的工况和粗糙度参数时,纵向粗糙度分布更有利于接触润滑.     

齿向修形对滤波减速器润滑性能的影响分析

综合考虑了滤波减速器齿向修形参数、真实齿面粗糙度和瞬态效应等因素,建立了轮齿混合润滑数学模型,数值计算了不同修形参数值对应不同啮合点的最大压力和中心膜厚,分析了齿面粗糙度和转速对润滑性能的影响.结果表明:修形参数r和Ry均存在一个优化范围,使得轮齿表面最大油膜压力显著降低,边缘效应弱化,而中心膜厚则随着r和Ry的增大而逐渐增大;未修形轮齿边缘油膜压力受粗糙度的影响而急剧增大,边缘效应更加显著,修形后轮齿的边缘效应得到了明显改善,因此,轮齿修形也因粗糙表面的存在而显得更加重要;随着转速逐渐降低,轮齿表面的平均油膜厚度逐渐变小,接触比逐渐增大,轮齿表面由弹流润滑逐渐转为混合润滑,最后演变为边界润滑.     

基于平均膜厚和压力流量因子的机械密封泄漏分形模型

研究接触式机械密封端面泄漏模型建立问题.采用分形参数表征机械密封端面形貌,通过引入压力流量因子来反映实际粗糙表面对泄漏通道的影响,推导出了压力流量因子的分形表达式,建立了基于平均膜厚和压力流量因子的泄漏分形模型.通过理论计算对机械密封泄漏率的影响因素进行了分析,并在自制的试验装置上对2套B104a-70型机械密封进行了试验,试验密封流体为20℃清水,压力为0.5 MPa,转速为2 900 r/min,弹簧比压分别为0.15和0.30 MPa.研究结果表明:泄漏率随着弹簧比压的增大略有下降,随着密封流体压力及转速的增大而增大,且端面越粗糙增大的幅度越大;当端面较粗糙时,泄漏率随着端面分形维数的增大或特征尺度系数的减小而迅速减小,而当端面较光滑时,泄漏率的变化很小;泄漏率的理论计算值与试验值吻合较好,特别是在进入正常磨损阶段后相差很小.     

机械密封织构化表面粗糙度效应的有限元模型与摩擦学特性分析

为研究粗糙度效应对机械密封织构化端面承载能力、摩擦学特性以及密封性能的影响，基于有限元法建立了机械密封的混合润滑模型，针对确定性非高斯随机分布粗糙表面，考虑润滑液膜的宏微观空化作用和粗糙峰的接触，研究了圆孔型织构化机械密封的摩擦学性能和密封性能. 结果表明:平衡状态下处于全膜润滑区间的非高斯型粗糙织构表面，在研究范围内均方根值越大的表面其减摩效果越好，且泄漏率越大，从混合润滑区到全膜润滑区转变的速度也随之增大；非高斯表面的偏态值和峰态值对密封表面的承载能力，润滑性能以及密封性能有影响，但无明显规律；在混合润滑区，圆孔织构具有增摩效果，而在全膜润滑区表现出减摩效果.      

考虑界面粗糙度动态变化的点接触弹流润滑特性研究

针对点接触弹流润滑的粗糙度效应,建立了考虑表面粗糙度动态变化的点接触弹流润滑模型,实现了油膜厚度和压力分布的快速求解. 对点接触弹流润滑下的粗糙表面弹性变形进行了定性和定量研究,同时分析了表面均方根粗糙度、载荷、相对运动速度和滑滚比对最小膜厚和最大压力的影响,以及表面形貌动态变化对膜厚比的影响. 结果表明:形貌变化改变了弹流油膜和压力分布特性,相对于光滑表面,表面粗糙度总体上提高了最大接触压力、降低了最小膜厚,在轻载工况下表面粗糙度对油膜厚度的削弱更加显著,而不同速度下粗糙度的影响程度基本相同,呈现线性变化趋势,膜厚比随载荷增大呈现先增后减的变化趋势,并在530 MPa左右达到峰值.      

多跨换热器传热管在静水中的阻尼特性

本文对多跨换热器传热管在静水中的结构阻尼与粘滞阻尼作了分析。通过求解带N—1个同心支承板的管子运动方程,导出计算挤压膜阻尼的公式. 为了验证理论分析结果,对四跨传热管在静水中的挤压膜阻尼进行了实验研究.得出挤压膜阻尼与支承板数、支承板厚度、管子—支承板管孔偏心率、以及管子—支承板间隙之间的关系.在分析实验结果的基础上,得出了计算传热管挤压膜阻尼的改进公式.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.     

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).