沉积粉尘的激波点火

采用瞬态阴影技术及红外光电传感器技术实验研究了沉积玉米粉的激波点火过程，并对此进行了理论分析。实验与理论分析结果表明，激波掠射沉积粉尘床后，粉尘颗粒先上扬到一定高度后才点火，颗粒的点火延迟时间与激波波前马赫数、气相氧气含量等因素有关。另外，沉积粉尘的激波点火延迟时间比相同条件下的悬浮粉尘激波点火延迟时间长。     

激波点燃粉尘的数值模拟研究

建立了用于模拟入射激波后可燃粉尘颗粒点火的一维非定常两相化学反应流模型，该模型考虑了气固两相间的相互作用、粉尘颗粒的加速、加热和化学反应。粉尘颗粒着火前的化学反应用发生在颗粒外表面和内孔表面的非均相反应描述，颗粒内部的温度变化用一含有化学反应源项的非稳态热传导方程来描述，以颗粒外表面温度的突跃上升作为可燃粉尘颗粒点燃的着火条件。我们用该模型和ＰＳＩＣ方法，对由中等强度激波从纯气相传入煤粉－氧气混合物而引起的非定常两相流动现象，包括气固两相间的相互作用、粉尘颗粒的加速、加热以及点火过程进行了数值研究，计算了对应于不同载荷比、马赫数为４～５的入射激波后煤尘颗粒的点火延迟时间，分析了由于可燃粉尘颗粒的存在，入射激波及波后气固两相流动参数的变化规律。数值计算结果与实验数据符合较好。文中建立的模型和所用的基于ＰＳＩＣ算法的数值方法，用最自然的方式描述气固两相流动，即用连续流模型（欧拉方程）描述输运相（气相）的流动，用轨道颗粒模型（拉格朗日方程）描述分散相（颗粒相）的运动。用这种方法模拟含尘介质中激波后颗粒的点火是很有效的，它可以清楚地确定哪一个颗粒群最先着火，它的初始位置以及在整个点火延迟时间内     

铝粉尘激波点火的实验研究

介绍了用铝粉氧化反应所生成的中间产物的特征光谱来测量铝粉尘受激波点火延迟的新的测量方法．测量了三种不同形状和粒度的铝粉尘激波点火延迟。由实验证实，当环境温度在Ａｌ２Ｏ３的熔点左右时，铝粒都可以被点火；由于机制的不同，点火延迟相差很大；点火延迟与铝粉颗粒的比表面积和活性铝的含量有关，而与环境的氧含量基本无关；但当氧含量小于１％时，铝粉尘不能被点火。     

激波与沉积可燃粉尘相互作用的实验研究

在新建立的激波管中实现了激波对沉积可燃粉尘的卷扬与点火。利用自行设计的瞬态阴影摄影系统和高空间分辨率的新型红外光敏二极管成功拍摄了激波后燃烧粉尘云的阴影图像和测量了激波后可燃粉尘的点火诱导时间。通过对实验结果的较深入的分析，揭示了激波与可燃粉尘床相互作用的基本特征     

点火准则和稀释气体对乙烯点火延时的影响

利用矩形截面激波管研究点火准则和稀释气体对乙烯点火延时的影响。采用压电传感器记录测点压力时间历程，采用光谱仪和光电倍增管记录自发光强时间历程，以压力、总自发光强与·OH和·CH自由基特定能级发射光强等信号判定是否发生自点火，给出自点火过程的时间起始点和终止点，得到了不同点火准则和稀释气体对应的乙烯/氧气/氮气和乙烯/氧气/氩气点火延时。结果表明:相同工况的乙烯点火延时测量数据相对误差约为15%，数据验证了本文实验和测量方法可靠性。针对当量比为1.0、压力为0.2 MPa，得到了温度范围为905～1 489 K，稀释气体的摩尔分数为75%氮气和75%氩气时的乙烯点火延时，给出点火延时和温度拟合的Arrhenius型表达式。不同点火准则会影响所测点火延时数据，但多次测量结果确定的点火延时和温度变化规律近似相同。不同稀释气体对激波管自点火流场的影响表现为和流场均匀性以及混合物比热相关。相同工况的乙烯/氧气/氮气点火延时大于乙烯/氧气/氩气点火延时。高温区和低温区的乙烯/氧气/氩气点火延时与温度的拟合关系不同，转折温度约为1 121 K。     

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程，研究了形 成爆轰波的点火特性. 对理想化学当量比氢气/空气混合气体，在初始压强20kPa的条件下， 马赫数2.6-2.8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的，第 2个点火区是由入射激波在抛物面上发生马赫反射引起的. 这种条件下流场中会出现爆燃转 爆轰，起爆点分别分布在管道壁面、抛物反射面和第2点火区附近. 起爆机理分别为激波管 道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相 互作用. 不同的点火和起爆过程导致了不同的流场波系结构，同时影响了爆轰波传播的波动 力学过程.     

瓦斯/空气预混气体爆燃流场测试系统多参数耦合同步控制实验方法

为了更精确地获得爆炸激波管内瓦斯/空气预混气体爆燃过程中，激波形成过程、压力和火焰传播速度以及火焰与惰性阻燃剂相互作用的流场演化图像。通过分析激波管测试系统中多个目标的时间响应特征及控制方式，利用超高速相机、光电倍增管、时间延时器、固态继电器、电荷放大器和数据采集系统等设备，设计实验方案，分别对激波管中瓦斯/空气预混气体爆燃高压点火系统的响应时间和惰性介质阻燃剂喷射系统的响应时间进行测试。实验结果表明电火花点火的响应时间为微秒量级，而阻燃剂喷射系统的响应时间为毫秒量级，以响应时间为依据，通过设置精确的延迟时间实现多目标同步控制，为完成激波管内瓦斯/空气预混气体爆燃过程的微观流场显示奠定基础。     

爆轰波管中铝粉尘爆轰的数值模拟

用两相流模型对爆轰波管中的铝粉尘的爆轰波进行了研究。模型考虑了气体和颗粒两相间速度和温度的不同及由于管壁引起的对流热传导和粘性引起的耗散,考虑了由于铝颗粒表面粗糙使得表面积增加的因素。铝颗粒的点火判据使用了新的判据,即铝颗粒在激波后的气流中温度达到铝的熔点且铝全部熔化即被点火。数值模拟了内径为15.2cm的爆轰波管中铝粉尘中爆轰波的传播和发展,得到了爆轰波速度及铝颗粒点火距离,还得到了爆轰流场中物理量的分布。从前导激波面到CJ面处,两相间的速度和温度有明显的差别。还考虑了粒子由于粗糙引起的表面积增加对爆轰波的影响,这个因素对铝颗粒的点火距离的影响较大,对这里计算的铝粉尘爆轰波的速度基本没有影响。结果表明,两相流模型可以较好地描述铝粉尘的爆轰过程,得到具有很粗糙表面、平均粒子直径为3.4m的铝粉尘浓度为304g/m3时爆轰波的速度为1.63km/s,点火距离为3mm,与实验值符合较好。     

粉尘火焰加速现象的实验研究

粉尘火焰的发生、加速及由爆燃向爆轰转捩的机理是个至今尚未弄清的问题。需要解决的技术关键之一是在实验室实现弱点火条件下的粉尘火焰加速直至达到爆轰状态。着手发展了一种球形喷粉扬尘装置，令产生的扬尘湍流在水平实验管中形成空间均匀分布和维持秒级悬浮的粉尘云状态。采用以上扬尘装置的水平实验管，在６ｇ黑火药的六点平面点火条件下获得了微细铝粉火焰经５ｍ长的传播过程加速至１０００ｍ／ｓ的实验结果。给出了扬尘湍流强度、粉尘粒度与浓度、点火能量及方式等因素对粉尘火焰加速过程中所起作用，及变截面效应（由小变大）对粉尘火焰减速的影响。     

激波诱导的燃烧粉尘云边界层的结构

足够强的激波扫过铺有可燃粉尘的界面时，波后形成燃烧的粉尘云边界层。为揭示其内部结构，本文对该现象进行了实验和理论两个方面的研究。理论计算表明燃烧粉云边界层可分为三个区域：诱导区、反应区和扩散区，诱导区的粉尘浓度最高，计算获得的粉尘云轮廓和点火延迟与实验结果基本吻合．     

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

正~~     

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

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

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

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

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。