相变材料接触熔化的研究

对有关相变材料接触熔化的研究成果进行了全面的综述, 内容包括:封闭空间的接 触熔化、围绕热源的接触熔化、滑动接触熔化、滚动接触熔化、冰的压力熔化、接 触熔化的量级分析与修正、由液体边界层内的热扩散引起的熔化、不连续接触点上 熔化与润滑、固体内部传热的影响、玻璃的软化与接触熔化.     

等壁温加热微通道内气体流动换热特性的直接蒙特卡罗模拟

用直接模拟蒙特卡罗方法对压力边界条件下气体在微通道内的流动换热特性进行了研究,给出了壁面与来流存在温差时的沿程速度分布特点,以及在可压缩性与换热条件综合作用下的温度分布特点。研究结果表明:微通道内气体可压缩性作用显著,温度分布由可压缩性和换热强度的相对强弱综合决定;壁面与来流存在温差时气体沿程速度分布型线在入口段内上凸;壁温高于来流温度时,气流速度与等温流动工况下的速度的相对大小与气体稀薄性有关。     

用边界积分方程法解具有自由边界的粘滞不可压缩流问题

文中提出用边界积分方程法解具有自由边界的粘滞不可压缩流问题。Navier—Stokes方程中的非线性项是通过线性化近似求解的。作为对这种方法的检验,计算了失重情况下容器内粘滞流体因表面张力变化(由温差引起)而产生的旋流,并将计算结果与用有限差分法所得结果作了比较。此外,还初步计算了在管内压力作用下液滴的形成过程。     

考虑动水压力和流固耦合的库水-坝体-地基系统建模与动力分析

结合动水压力模型和罚函数耦合算法，考虑地基和坝体结构的接触以及边界效应，构建动水压力和流固耦合作用下的库水-坝体-地基地震响应分析模型和方法。通过与试验结果及解析解和实测数据对比，验证了本文模型和方法能准确反映系统地震荷载和分析整体耦合系统的动力响应，引入罚函数处理流固耦合界面能提高计算收敛速度。进一步以某重力坝工程实际为背景，验证本文构建的模型和方法适用于库水-坝体-地基耦合系统动力分析的可行性，并分析了地震作用下地基变形对系统动力响应的影响。     

柔性转子弹性对行波旋转超声电机接触特性的影响

转/定子间动态接触过程决定行波旋转超声电机(TRUM)的输出性能,但是柔性转子弹性在以往界面接触特性研究中常被忽略.本文中基于建立的考虑柔性转子和摩擦层三维弹性的转-定子动力学耦合模型,从接触界面的接触状态、速度、接触应力、摩擦应力驱动和阻碍子区等空间分布角度,研究了柔性转子弹性对TRUM摩擦界面的动态接触行为和摩擦损耗的影响.结果表明：柔性转子在轴向、径向和切向均呈现出一定的振动速度,且在接触界面沿圆周方向与定子各向速度形成一定的相位关系,显著影响界面动态接触特性、摩擦特性和摩擦损耗;柔性转子与定子轴向位移分布的相位一致性表明,定子与摩擦层的接触界面具有类齿轮啮合机制;与以往中间驱动两端阻碍的分布结果不同,考虑转/定子动力学耦合的接触区摩擦应力分布为一端驱动和一端阻碍的新机制.计算结果表明在建模过程中忽略柔性转子弹性得到的界面摩擦损耗要比考虑柔性转子弹性获得的摩擦损耗大得多,主要原因：(1)柔性转子径向刚度小,摩擦力作用下接触区的运动方向与定子相同,从而减小了径向速度差和径向摩擦损耗;(2)柔性转子界面的切向振动特性,改变了接触界面黏-滑子区分布,使得行波远离端的摩擦应力阻碍子区消失...     

温度梯度与缺陷共同作用下的双壁碳纳米管持续驱动器件

纳尺度驱动因其在能量转换、物质输送等众多领域的潜在应用得到广泛关注，但实现长程持续驱动依然在原理上面临重大挑战．本文借助双壁碳纳米管的层间范德华相互作用，以热梯度力为基本动力，利用外管缺陷对层间热梯度力的调控作用，提出了一种新型的纳尺度可控持续驱动原理．其基本机理是：外管在周期性缺陷和温度场共同作用下形成周期性的非对称温度梯度，在内管上产生周期性的非对称热梯度力，使得每个周期（驱动单元）上都存在净驱动力，从而推动内管实现长程持续运动．同时，内管的运动方向或速度可通过外管上的控温区域或温差控制．     

行波型超声波电机摩擦特性的实验研究

超声波电机是通过逆压电效应产生高频振动并通过摩擦实现驱动的电机,其定转子摩擦接触面的高频、微幅的行波振动有其特殊性,跟以往的摩擦接触显著不同,并对电机的换能效率、电机寿命和摩擦材料的选择等都有影响.本文建立了超声波电机定转子接触面摩擦系数测量装置,通过试验比较了摩擦接触面在常规、定子驻波振动以及行波振动下的动态摩擦系数,分析了定子相对转速(速度)、预紧力以及温度对定转子接触面摩擦系数的影响.根据摩擦二项式和电机定转子表面弹塑性接触状态对摩擦系数与预紧力等的对应关系进行了说明,分析了驻波振动和行波振动的超声减摩效果以及温度对摩擦系数的影响,确定了适合行波型超声波电机的摩擦系数计算模型.超声波电机的驱动摩擦特性是1项基础性工作,有助于行波形超声波电机设计、建模和驱动控制等.     

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

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

供油温度对润滑状态影响的试验观察与数值模拟

采用球-盘点接触光干涉润滑油膜测量装置,观察了不同供油温度下的接触区润滑状态.结果表明,随着供油温度的升高,润滑油黏度降低,润滑油成膜能力整体下降,接触区润滑状态在较大速度范围内处于边界润滑和混合润滑区间;温度升高还将导致膜厚曲线斜率增大,壁面温差是其诱因,通过数值模拟证实了壁面温差的影响;速度增加使接触区润滑状态向等黏度-刚性区间转化;温度和载荷使接触区润滑状态趋向于压黏-弹性区间.系统散热也会对油膜厚度产生影响.     

磁驱动等熵压缩和高速飞片的实验技术（续2）

4磁压加载等熵压缩实验的数据处理方法 4．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).