等离子熔射粉末颗粒飞行过程格子Boltzmann法仿真

为考察等离子熔射过程中粉末的飞行过程,本文在已开发的正六边形7-b it格子Bo ltzm ann(LB)方法等离子射流的温度场和速度场的计算模型基础上,采用单个颗粒加速方程,建立了一个随机算法,实现了对粉末颗粒在射流场中运动过程的仿真;计算结果通过动画演示了粉末飞行的全过程,表明初始位置越靠近射流场出口中心的粉末颗粒加速越充分,并且在射流场一定的情况下,减小粉末颗粒直径可以提高粉末速度,但会降低粉末利用率。     

膛内等离子体点火及燃烧增强过程数值模拟

在初始点火过程中 ,假设等离子体为充分发展射流流动 ,因而可采用积分近似模型进行描述。而在随后的膛内燃烧过程中通过以经验公式给出等离子体混合长度的方法来处理等离子体的运动 ,并给出了合适的相间阻力及传热公式。通过在计算网格设计中引入了自适应方法提高了数值解精度。计算膛压曲线与实验符合良好。相对常规内弹道情形 ,等离子体能量的注入使在不同时刻膛内的气、固两相速度、压力和膛底温度等出现提高 ,这是等离子体注入后对内弹道的增强过程。计算结果为固体工质电热化学炮的弹道设计提供了理论依据。     

出口条件对圆形湍流射流远场区自保持性的影响

报道出口条件对圆形湍流射流自保持性影响的实验研究结果. 对来自渐缩和长管两种不同结构喷嘴的射流,在相同雷诺数条件下,沿轴线进行了速度测量; 研究的统计量包括平均速度、湍流强度、高阶矩、能谱和积分尺度. 实验结果表明,渐缩喷嘴射流比长管射流发展得更快、更容易达到自保持状态. 通过对比发现,在两射流的速度(温度)场中,平均速度(温度)、湍流强度、偏斜因子和平坦度因子都存在明显的异同. 同时发现两射流的积分尺度随轴向距离的增加都成线性增长,且在渐缩喷嘴射流中增长得更快. 通过对比两射流的边界层厚度、径向与轴向湍流强度的比值、湍动能能谱图并结合前人的研究结果,对两射流湍流场所表现出的不同的统计学行为给出了合理的解释.      

纳米W粉冲击烧结的分子动力学模拟

粉末冲击烧结是制备高品质W的一种有效方法，而分子动力学方法在尺度极小、过程迅速的数值模拟上有着独特的优势。因此运用分子动力学方法，结合W的嵌入原子势，对常温下的纳米W粉末的冲击烧结过程进行模拟，得到颗粒微观压实过程图、体系速度分布云图、p-U_p、T-U_p、T-p曲线以及径向分布函数。研究了不同颗粒速度及产生的射流对纳米W粉末冲击烧结影响，分析了微观冲击烧结机理。结果表明，低速冲击条件下（500 m/s以下），纳米颗粒无法压实。高速条件下（1 000 m/s及以上），颗粒能获得致密化很高的压实。颗粒间的相互挤压造成的高应力使颗粒表面的原子发生流动变形，原子向颗粒间空隙流动，形成压实。颗粒间产生的射流以及高速冲击导致的颗粒熔化，均促进烧结获得致密度更高的烧结体。     

聚能射流形成数值模拟

基于网格自适应方法,ALE方法和SPH方法模拟了聚能射流形成过程.分析了射流头部速度随时间的变化关系,初始时刻位于药形罩不同位置的金属形成射流后的速度分布情况以及射流形成后沿射流方向射流整体速度分布情况.对比了三种方法用于模拟聚能射流形成问题的特点以及与实验的近似程度.数值模拟结果与射流形成过程与机理一致,可以体现射流形成特点,且SPH方法计算结果与实验数据最为接近.     

托卡马克中凯尔文-亥姆霍尔茨不稳定性与密度分布及速度分布的关系

近年来,在托卡马克装置实验中观察到等离子体环的密度及环向速度沿半径都有变化。当反平行中性束注入时,密度分布对应空洞型温度(即中心温度低而外侧高),等离子体的稳定性变坏。本文对密变分布和速度分布对等离子体环的凯尔文-亥姆霍尔茨稳定性的影响分别进行了研究,得出,在托卡马克实验参数范围内,密度沿半径的分布会使凯尔文-亥姆霍尔茨稳定性变坏,从而为反平行注入对等离子体约束有明显破坏提供了一种可能的解释。     

消息与动态(九)

一种计算等离子喷涂颗粒吸收热量的方法计算等离子喷涂过程中颗粒吸收热量的方法一般直接引用能量和质量转移联立方程的解.可是被喷涂的颗粒材料在等离子射流中温度变化十分激烈,变化范围太宽,以致 ...     

双环状电极大气压等离子体射流的时间分辨诊断

大气压等离子体射流在众多领域中发挥着至关重要的作用,其中,拥有双环状电极结构的射流装置表现十分突出,此类装置具有结构简单、放电稳定的优势,并且产生的射流长度较长、温度较低.本工作利用ICCD相机采集技术,对高频高压交流电源驱动的双环状电极等离子体射流进行了纳秒级时间分辨诊断.研究发现,在一个完整的放电周期内,存在两个明显的放电阶段,分别位于外加电压的正、负半周期的峰值附近,并且在正半周期的放电阶段中,亮度、射流长度、放电发展速度均明显强于负半周期,值得注意的是,仅在负半周期的放电阶段中,才能观察到等离子体离开管口自主向前传播的现象.本工作对于了解等离子体动力学过程,揭示等离子体的行为规律及优化等离子体设备等诸多方面有积极地推动作用.     

消息与动态

一种计算等离子喷涂颗粒吸收热量的方法计算等离子喷涂过程中颗粒吸收热量的方法一般直接引用能量和质量转移联立方程的解.可是被喷涂的颗粒材料在等离子射流中温度变化十分激烈,变化范围太宽,以致     

消息与动态(九)

一种计算等离子喷涂颗粒吸收热量的方法计算等离子喷涂过程中颗粒吸收热量的方法一般直接引用能量和质量转移联立方程的解.可是被喷涂的颗粒材料在等离子射流中温度变化十分激烈,变化范围太宽,以致 ...     

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

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

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

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