金属圆管内爆轰波相互作用效应的数值模拟研究

采用动力学有限元程序,针对大变形带来的网格扭曲造成计算的锁死现象进行了ALE技术处理,结合高速摄影和X光照相实验结果,对金属圆管内爆轰波的相互作用效应进行数值模拟,给出了爆轰波对碰作用下金属圆管运动计算结果,复现了金属圆管对碰区鼓包现象,计算结果与实验吻合较好。可以预见,基于该方法能够初步分析影响金属圆管运动规律的诸多因素。     

空心弹高速入水机理及特性数值模拟研究

为分析空心弹高速入水的机理及其特性,基于雷诺时均Navier-Stokes方程、VOF(volume of fluid)多相流模型、Realizable k-ε湍流模型,引入Schnerr-Sauer空化模型和重叠网格技术对空心弹高速入水进行数值模拟研究,获得了通孔孔径和头部形状对空心弹的空化特性、空泡形态和入水运动特性的影响规律。研究显示数值计算的空泡形态和入水速度、位移曲线与实验结果吻合较好,验证了数值模拟方法的可行性。结果表明：当通孔孔径不同时,通孔孔径越大,空化现象越明显,通孔射流越长,但对空泡半径的影响不大;通孔孔径越小,空泡闭合时间越早,与水面碰撞产生的阻力系数峰值越高,空心弹入水稳定后其阻力系数也越大;无量纲直径在0.575～0.600之间时,空心弹的运动最为稳定。当头部锥角不同时,头部锥角越大,空泡直径越大,空化现象出现得越晚,但空化生成的速度更快;随着头部锥角的增大,阻力系数变大,空心弹的速度衰减变快,相同时间运动的距离较短;头部锥角越大,俯仰角的变化越小,空心弹的运动越稳定。     

基于FTM的Kelvin-Helmholtz不稳定性衍化数值分析

采用FTM研究了两相流体界面中的K-H不稳定性,分析了初始振幅、波数、表面张力、密度比、速度梯度、韦伯数We对K-H不稳定性发展的影响。结果表明:增大初始扰动对K-H不稳定性具有促进作用,扰动越强,K-H波浪翻卷高度越高;波数的增加能加快K-H不稳定性的衍化过程,但对K-H波浪高度的发展影响较小;减小表面张力或密度比都能使不稳定性的发展速度增加;初始速度梯度对不稳定性的发展有较大影响,梯度越大,K-H波浪翻卷高度越高;We对K-H波浪的对称性影响较大,We越小,对称程度越高。     

爆轰波对碰区产物驱动金属圆管的研究

利用闪光X光照相、光学高速分幅照相实验,测得了金属圆管在两个爆轰波对碰后产物驱动外壳的膨胀、变形过程,该过程表明爆轰波对碰作用使得对碰处金属圆管相对贯穿断裂时间明显提前,对碰处金属圆管径向膨胀速度较非对碰处明显增加。利用拉格朗日二维流体动力学程序TTD2C,数值模拟了对碰驱动过程,计算结果与实验结果符合很好,且用Taylor断裂判据得到的金属圆管相对贯穿断裂时间也明显提前。     

表面粗糙度对6061铝合金薄壁管冲击膨胀断裂性能的影响

采用分离式霍普金森杆实验技术,对表面加工后不同粗糙度的6061铝合金薄壁圆柱管进行动态膨胀断裂冻结回收实验,并对薄壁金属圆柱管动态膨胀断裂过程中裂纹萌生、扩展情况以及最终断裂模式等进行了研究。结果表明:相同冲击压力条件下,薄壁金属圆柱管表面粗糙度越大,材料越容易发生膨胀破裂;裂纹萌生于外壁面,由外向内扩展,并且裂纹的扩展主要受裂纹处应力状态的影响;薄壁金属圆柱管的断裂模式由拉伸和剪切断裂机制起主导作用,其断口为拉剪混合型断口。     

重力和接触角对表面张力贮箱内液体流动的影响

采用流体体积函数方法数值模拟板式表面张力贮箱内液体的流动过程. 主要考虑了不同重力加速度和接触角等因素的影响. 发现在接触角等于10° 的前提下, 重力加速度小于10-2g₀, 液体沿着外侧导流叶片爬升至贮箱顶部; 当重力加速度小于10-3g₀ 时, 重力加速度的减小对液体爬升速度的影响较小. 在爬升过程中, 导流叶片附近液体的速度呈逐渐减小的趋势, 且重力加速度越小, 初始速度越大. 接触角对液体的爬升高度也会产生影响,接触角越小, 液体爬升至顶部的时间越短; 当接触角大于45° 时, 外侧导流叶片附近的液面不能爬升至顶端. 另外, 接触角越大, 液面高度的斜率越小,导流叶片附近液体的初始速度也越小.      

扫路机气力系统工作效率研究

本文旨在探讨提高吸式扫路机气力输送系统的效率,采用计算流体力学与离散单元法对圆柱颗粒在气体中的运动进行数值模拟,分别研究了颗粒密度、颗粒长径比、颗粒初始平面与吸入口的距离、吸入速度、颗粒所处圆周对于颗粒吸入效率的影响。研究结果表明：颗粒密度越大,颗粒的吸入效率越低;颗粒长径比对吸入效率的影响较小;颗粒初始平面与吸入口的距离与颗粒吸入效率为负相关;吸入速度增大有利于提高颗粒吸入效率,但是当吸入速度增大至一定程度时,其对吸入效率影响减小;颗粒所处圆周与吸入口圆周径向距离越小,颗粒的吸入效率越高。     

航行体出水破冰的多场耦合效应与相似律

航行体出水破冰中的耦合效应及载荷特征, 是出水冰结构安全性评估的重要依据. 针对航行体出水破冰问题, 通过量纲分析, 获得了影响航行体动载荷及头部应力的主控参数和相似律. 基于LS-DYNA流固耦合计算方法, 得到了航行体在不同冲击速度、冰层厚度、冰层大小条件下的载荷特性. 计算结果表明, 航行体速度越大, 不同冰层对其过载和头部应力的影响差别越大, 这主要是因为航行体速度越大, 通过水介质对不同冰层的前期破坏程度不同. 对于无限大冰层, 当其厚度大于3倍航行体直径时, 航行体穿冰后期呈现稳定侵彻现象, 航行体的过载和头部应力只与航行体的速度和冰的动力学性能相关; 而对于薄冰, 航行体速度越大, 其头部应力反而越小, 这是因为航行体初速度越大, 其通过水的运动对冰的前期冲击破坏越严重, 冰层易开裂上鼓, 所以造成航行体头部应力较小. 对于径向尺寸为6倍航行体直径的碎冰, 当其厚度大于5倍航行体直径时, 碎冰对航行体运动特性的影响和无限大冰层几乎相同; 而当其厚度小于3倍的航行体直径时, 只有在初速度较低时, 碎冰的尺寸效应才可以忽略. 此外, 对比碎冰和无限冰层对航行体运动的影响可以看出, 越厚的冰受前期水的冲击破坏越小, 碎冰和无限冰层的影响规律基本一致; 而较薄的冰在前期水的冲击下破坏严重, 碎冰和无限冰层对航行体运动的影响都较小; 只有中等厚度的冰, 在较高冲击速度下碎冰和无限冰层才表现出径向尺寸效应相关的破坏程度, 如无量纲厚度为3的两种冰在航行体较高初速度40 m/s的条件下前期破坏差别较大, 导致后期对航行体运动特性的影响具有显著差异.      

爆轰波对碰驱动下金属圆管膨胀变形特性研究

采用高速分幅摄影和激光速度干涉仪（velocity interferometer system for any reflector, VISAR）联合测试技术,对HR-2钢管在爆轰波对碰驱动下的膨胀变形特性进行了研究。实验得到了圆管膨胀变形的过程图像以及圆管特定部位径向膨胀的位移-时间曲线和速度-时间曲线。采用LSDYNA-3D程序对HR-2钢管的膨胀变形过程进行了数值模拟,计算结果与实验结果吻合较好。     

三维气泡与自由表面相互作用的直接数值模拟

采用VOF中的PLIC界面重构方法数值模拟了三维气泡在液流中上升并与自由表面相互作用的运动.分别考察了不同初始高度,有无来流及有无再生气泡对气泡上升高度、上升速度、压力及与自由表面相互作用等的影响.结果表明:气泡初始位置越低,顶端上升的高度越大,自由面隆起的范围更广.越靠近自由表面,底部射流横向发展越窄,而向上的压力梯度,气泡上升速度,底部射流上升高度越大,反之则反;但如果底部射流均在接近自由表面以前已横向发展充分,则差别不大.气泡外形、上升高度、破裂时间以及上升速度与来流无关.产生再生气泡后,原生气泡与再生气泡相吸,相互加速对方的上升;自由表面抬升的高度增幅较大,范围拓宽,上升速度也大大增加,且再生气泡越多,自由表面隆起的范围越大.     

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

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

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

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