常规螺杆泵定子有限元分析

螺杆泵定子橡胶不仅是易损构件，而且它与转子的配合状况对螺杆泵的工作性能影响显著。目前，还没有能够直接对实际工况下的定子橡胶的变形和受力状态进行测试的有效手段，因此对螺杆泵定子进行有限元分析自然是有益的尝试。重点研究定子橡胶衬垫内轮廓线的变形规律，而应力和应变分布情况仅作为辅助性的分析；对螺杆泵工况的考虑方面，重点研究均匀内压作用的情况，再考虑工作压差以及定子与转子的装配等所引起定子橡胶的变形和受力的变化。利用有限元分析软件Abaqus，对螺杆泵定子进行研究，得出了螺杆泵定子在不同工况下的受力状态和变形规律以及定子材料参数对变形规律的影响。利用此规律可对螺杆泵进行优化设计，提高其工作效率。     

等壁厚螺杆泵脱胶机理数值模拟研究

等壁厚螺杆泵作为一种高效的人工举升装置已经被广泛地应用于采油工程中.相比于常规螺杆泵,等壁厚螺杆泵具有定子温胀溶胀均匀,运行稳定,举升效率高等优点.但是在实际工作过程中,等壁厚螺杆泵定子容易脱胶,导致定子失效.论文采用有限元数值模拟方法,建立螺杆泵三维动力学模型.通过模拟相同结构参数常规螺杆泵和等壁厚螺杆泵运行过程中的应力应变场、压力分布规律,对比分析发现由液压压差和定转子过盈接触引起的剪应力是造成等壁厚螺杆泵脱胶的主要因素.在此基础上研究了附胶率以及各模型参数对剪应力的影响规律.研究发现,考虑附胶率后,在井下复杂的工作环境下,剪应力会大幅增加,等壁厚螺杆泵定子容易发生脱胶.这些结果对等壁厚螺杆泵的优化设计具有指导意义.     

SMA复合材料在热载条件下的残余应力分析

本文基于三相复合圆柱模型发展了增量型的分析方法,讨论在ＳＭＡ复合材料中由于ＳＭＡ材料相变以及各相材料热特性随温度变化引起的残余应力。研究基体与过渡恸介面和纤维与过渡界面间的残余应力,同时讨论由于基体相的变化对残余应力的影响。特别研究了涂层和复合材料基体间界面处的残余应力受纤维体积比、涂层厚度、纤维最大相变应以及基体中纤维取向等影响,而且讨论了计及应力对相就运动方程的影响时对ＳＭＡ复合材料相变温度和     

矩形变截面梁的挠曲分析

阐述弹性均质矩形变截面梁挠曲线的近似解法。梁的高度的长度方向呈线性变化，或是二次型抛物线变化，其宽度于线性变化。这种近似解法适用于集中荷载、均布荷载、线性分布荷载以及抛物线型分布荷载的作用情况。     

炭黑填充橡胶材料改进Mooney模型

本文讨论了炭黑填充橡胶材料的唯象本构模型。考虑到Mooney模型无法表征橡胶类材料大变形阶段的力学特性,首先利用实验数据,对Mooney模型进行了分析,讨论了炭黑含量与Mooney模型准确表征橡胶材料应变区间大小的关系,Mooney模型对纯剪切和等比双向拉伸等复杂变形的预测能力,同时也分析了材料参数对Mooney模型的影响。最后在Mooney模型的基础上添加了一个修正项,且改进后的Mooney模型满足Treloar和Ogden六项假设。通过与实验数据对比分析,改进Mooney模型可以较好地描述橡胶材料大变形阶段的应力应变关系,同时提高了预测橡胶材料复杂变形的能力。     

金属橡胶材料基于微弹簧组合变形的细观本构模型

研究了金属橡胶材料细观本构模型。依据金属橡胶材料变形的主要特征,分析了微弹簧沿轴向和径向变形的规律,分别推导了应力应变关系的解析表达式,由两种微弹簧的组和变形构造了代表性体积单元-微元体。为了反映因铺层引起的工艺各向异性,本构方程中引入了一个新的材料参数铺层比例系数iβ,该参数可以较好地表征细观结构中微弹簧取向的分布状况。在大量试验的基础上,提出了确定铺层系数的实验方法,进而建立了包含金属橡胶材料细观结构信息的本构方程。理论与实验结果比较表明,本文建立的本构模型能较好的反映材料的力学行为。     

岩土介质非稳态热固结耦合问题的热源函数法

考虑耦合效应的饱和土体热固结问题控制方程,利用Fourier变换、Laplace变换给出其在变换域上的解,将初始温度场分布视为虚拟的热源或者将热源等价为特定的初始温度分布,利用热源函数法给出瞬时线热源非稳态温度场、应力场和位移场的解析求解方法,通过在时间域和空间域上进行积分,给出有初始温度场分布以及有分布内热源存在且热源强度随时间变化条件下的热固结问题计算方法。对一无限大物体内存在有平面矩形域热源情况下周围介质的温度、孔隙水压力以及位移等的变化特征进行分析。研究表明,热源函数法可有效地求解一系列复杂情况下的热固结问题。     

深水钢悬链线立管疲劳寿命计算方法

基于悬链线理论及立管与海床接触模型,建立了一种估算立管疲劳寿命的方法,模型中考虑了深水浮式平台运动过程中立管悬链线形状的变化及触地区域海床接触刚度大小对疲劳寿命的影响。分析表明:浮式平台的运动会引起立管悬链线形状的改变以及立管内部应力的重新分布,这是引起立管疲劳损坏的主要原因;随着海床土体刚度的增大,立管疲劳寿命逐渐降低,为了获得准确的立管寿命,有必要对海底土体进行更为全面细致的研究。通过与现有分析方法比较,本方法具有计算方便快捷,结果可靠等优点。研究可为钢悬链线立管的设计分析提供一定的参考。     

固体火箭发动机燃烧室光弹性分析

本文运用光弹性法研究了某型固体火箭发动机燃烧室内端燃型药柱和星柱型药柱两种形状的固体复合推进剂的残余应力场,测定了光弹材料的弹性模量、泊松比、线胀系数、材料条纹值和冻结温度,获得了光弹材料在压应力和温度应力下的等差线条纹图,分析了应力分布情况.     

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

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

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

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

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


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).