时间推进方法在叶轮机械内部流场计算中的进展

时间推进方法目前广泛应用于叶轮机械内部流场的计算,这种方法可以划分成显式时间推进方法和隐式时间推进方法．本文简要介绍了时间推进方法用于叶轮机械内部流场计算的起源及这种方法的优点和缺点,重点记述了国内外叶轮机界应用显式时间推进方法和隐式时间推进方法计算叶轮机内部流场的最新进展,并简要介绍世界上一些著名的研究机构开发的时间推进方法计算程序所采用的数值计算方法      

含分流叶片的离心压缩机级内三维流场数值分析

使用Fine/Turbo三维粘性计算程序对含分流叶片半开式叶轮的离心压缩机级内三维粘性流场进行数值分析.该程序应用Jameson的中心差分格式结合Yang&Shih k-ε湍流模型,使用时间推进法求解雷诺平均N-S方程,用四阶显式Runge-Kutta法进行时间推进.计算中CFL数取2.5,残差收敛至10-4数量级、计算域进出口质量流量误差小于0.05%时认为收敛.并为进一步的产品优化设计及改进研究打下基础.计算结果表明该模级的数值模拟曲线与实验值吻合良好;叶轮出口截面上两个通道的速度分布极不均匀,右通道的速度比左通道的变化剧烈.为了加快收敛,计算时采用了隐式残差平均法及完全多重网格技术.     

蜗壳内紊流流场的有限元法计算

用有限元法计算径流式叶轮机械蜗壳的紊流时均流场.有关紊流模型采用K-ε两方程模型,用关于压力p的罚函数方法求解.所得结果可供分析蜗壳流场用.文中方法也可用于计算其他形状的二维通道流动.     

离心风机子午通道内湍流场数值模拟

由进风口－叶轮－无叶扩压器－蜗壳等部件组成的离心风机通道内流分析是非常复杂的，目前还只能是分别计算各部件内的流场，但必须考虑部件间的相互影响。本文采用轴对称Ｎ－Ｓ方程，根据三维叶轮通道计算给出的叶片力分布，求解了考虑叶片力的进风口－叶轮－无叶扩压器组成的子午通道问题，所得结果可用来给出三维叶轮通道计算的进口条件，并可用于优化设计进风口及叶轮前、后盘形状。该方法已得到实践检验。     

一组低转速径流式压气机叶轮三元流动参数分布的计算及与实验数据的对比

本文采用构造叶片面的上下游延伸流面,在接长了的叶轮通道内,用准正交面法计算径流式压气机叶轮三元流场。用于计算S.Mizuki低转速叶轮二个流量系数情况下的流场,对无分离的B型叶轮,计算结果与实验数据符合得相当好;对分离严重的A型叶轮,计算结果亦可对分离区的产生提供合理的解释。本方法可应用于径流式压气机叶轮变工况流场分析。     

水泵叶片表面三维边界层计算

本文用有限差分法计算混流式可逆水力机械水泵工况叶片表面的三维边界层。水泵叶轮中主流区的三维势流由直接边界元法计算。对于叶片面附近的粘性流动。用三维半正交坐标系中的边界层方程表示。为了提高计算精度采用贴体坐标技术生成边界层区域的计算网格。并利用Cebeci等变换函数及Keller差分格式离散方程。用分块解法求解。计算叶轮叶片表面的压力分布与相应试验结果进行了对比。     

弹塑性有限元法在机械工程中的应用——介绍轮盘的超速预应力处理计算

本文叙述了旋转机械零部件的超速预应力工艺及其计算方法。文中采用了自编弹塑性有限元分析程序,对某高速燃汽轮机的试验叶轮分别用变刚度法以及初应力法进行了计算,得到了满意的结果。采用本文推荐的方法及程序,对具有解析解的标准叶轮,进行了比较分析,其最大百分相对误差不超过3％,因此所荐方法及程序是超速预应力工艺计算中较为有效的分析工具。     

级环境下叶片扩压器流场的实验与数值研究

为了研究离心压缩机级环境下的非定常干扰的基本流动现象,并验证多级叶轮机械的CFD软件的分析能力, 对一大尺度离心压缩机的叶片扩压器流场进行了实验测量和数值计算. 实验采用了固定热线、相位锁定------系综平均技术,用常温热线风速仪对叶轮后的叶片扩压器通道内不同周向、径向和轴向位置处的非定常速度进行了测量,同时提出了非定常强度的概念,以定量考核非定常的影响.实验结果表明, 叶片扩压器内的非定常流动非常复杂,其时间周期并非叶轮叶片通过时间,随着与离心叶轮之间的距离增大,非定常扰动逐渐减弱,但一直延续到叶片扩压器的出口.另外,对该实验压缩机级开展了两个不同的数值计算,并与实验数据进行了比较：定常数值计算软件采用了作者发展的确定应力模型,非定常数值计算是用商业软件NUMECA实现的,计算采用了滑移界面技术. 两个计算结果与实验在扩压器的进口截面处吻合得很好.     

计算离心泵叶轮流场的扫描流束有限元方法

本文提出一种求解离心式叶轮流场的数值方法，将流动求解区域离散为有限个由流线构成其边界的单元，采用伽辽金法建立的单元方程在一条流束上集合为方程组，流线上的节点坐标亦作为未知量包含在有限元方程中，通过扫描计算，逐步解得流线位置及流动参数。本文应用叶轮的通流理论流动模型，采用扫描流速有限元方法对离心泵叶轮流场进行了计算，并与有关文献作了比较。     

各向加权异性张力有限元法及其应用

本文利用各向异性张力有限元的特点设计了灵活而有效的各向加权异性有限元算法，采用这种方法，计算了离心叶轮内部准三元流场。     

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

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

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

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

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

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