基于气炮实验的PTFE/Al复合材料冲击反应阈值

基于16 mm口径气炮撞击实验,对铝颗粒增强的聚四氟乙烯(polytetrafluoroethylene)/Al(PTFE/Al)冲击反应复合材料的冲击反应阈值开展了研究。为研究不同撞击加载条件下应变率和碰撞应力对PTFE/Al冲击诱发反应的影响,实验中采用铝、钢和低密度聚乙烯(low density polyethylene, LDPE)这3种不同材料的靶板及不同长度的试样,进行不同加载条件下的测试分析。实验结果显示,PTFE/Al材料的冲击诱发同时受到碰撞压力和加载应变率的影响。同时,通过对试样撞靶过程进行数值模拟,并与实验和理论结果进行对比。基于实验数据,拟合出PTFE/Al材料冲击反应的的预测曲线。     

超高速撞击厚靶过程的能量分配研究

超高速撞击过程的能量分配研究,对于解决动能撞击、发展导弹拦截技术、判定空间飞行器被撞事件及评估碰撞破坏程度具有重要的理论意义。本文在总结前人关于超高速撞击过程能量分配的基础上,将超高速撞击厚靶过程中弹丸的动能分配归纳为靶板的变形能、弹丸与靶板作用过程应力波传播使靶板内能的增加、撞击产生碎片的崩溅能和产生电磁辐射的辐射能,并结合理论推导、实验和数值模拟对撞击速度为2.61km/s且正碰撞2A12铝靶的能量分配进行了定量计算。研究结果表明：无论在弹坑的形貌、尺寸还是辐射温度等方面,实验测量结果、理论推导结果与数值模拟的结果均基本吻合。该研究成果在解决行驶中的车辆碰撞问题以及航空飞行器遭遇鸟撞等领域亦有重要的参考价值。     

冰弹丸高速发射方法研究

冰弹丸的高速发射技术是认识冰弹丸高速撞击靶板损伤效应的瓶颈技术。利用二级轻气炮作为发射装置,采用液氮杜瓦制冷系统实现温度控制,选择铝质膜片对气压进行控制,利用高速摄像机对冰弹丸的完整性进行观测,并利用丝网测速仪对弹丸速度进行测量,实现了对冰弹丸的高速发射。通过实验验证了该方法的可行性,得到了完整冰弹丸在1516m/s下撞击铝合金靶板的实验结果。     

材料动态屈服应力的实验测定

本文提出了用平头柱形弹丸撞击刚性靶后的残余总长度为实测特征量的一个测定的新方法,在碰撞速度大于或小于弹丸材料塑性波速的两种情况下,给出了推算公式,得到材料动态屈服应力的稳定数据。     

飞机加强蒙皮在12.7 mm弹丸撞击下的变形与破坏

为了研究飞机蒙皮在12.7 mm标准机枪弹丸射击下的损伤,对3 mm厚LY-12 CZ材料的单蒙皮及其加筋板进行了模拟弹击试验。通过试验研究,建立了一个由高速气炮、弹体与弹托分离机构、连续位移激光测速装置和弹丸回收装置组成的系统并被有效地用于弹丸正撞击试验。通过对四边固支的3 mm厚蒙皮用12.7 mm直径弹丸进行速度约60~300 m/s的正撞击试验,结果表明,靶板从微小损伤到完全击穿;弹击造成的变形区有效直径随弹丸速度的增大呈幂指数趋势下降;弹击引起的变形深度随弹丸撞击速度的增加呈直线下降;靶板上的应变随弹丸速度的增加逐渐降低。弹丸剩余速度随弹丸撞击速度的增加呈直线上升。最后利用DYNA3D程序对单蒙皮及其加筋板进行了弹击数值模拟,模拟结果与弹击试验结果较吻合。     

靶板厚度对卵形弹丸垂直贯穿中等厚度混凝土靶的影响

为研究卵形弹丸贯穿中等厚度混凝土靶体的贯穿规律,开展直径60 mm尖卵形弹丸贯穿不同厚度混凝土靶体的侵彻实验,获得了不同撞击速度的弹丸贯穿不同厚度混凝土靶体的剩余速度规律。结合无网格SPH方法、RHT混凝土本构以及状态方程,对贯穿实验进行数值模拟,对不同工况下的弹丸过载规律以及靶体的损伤过程的分析发现:弹丸贯穿中等厚度混凝土靶体的贯穿过程分为开坑阶段、隧道稳定侵彻阶段以及靶背影响出靶阶段,在相同初始撞击速度下的靶背影响区的厚度随着靶体厚度的增加而增大。实验结果与数值模拟结果对比,表明模型能够有效模拟弹丸贯穿混凝介质问题,研究结果可为贯穿机理的研究提供参考。     

评价炸药安全性能的苏珊试验

苏珊试验是测量炸药撞击感度的一种射弹试验。将炸药装入壳体中,用空气炮或火炮以30～300m/s,甚至更高的速度发射,弹丸约经3.7m,与钢靶相撞。本文简介了试验的基本方法,用氦氖激光器、接收转换器、测时仪构成的光电系统测量弹丸的飞行速度;用红宝石序列脉冲激光器、纹影仪以及等待型转镜扫描相机构成的高速摄影系统拍摄弹丸的着靶姿态以及撞靶后挤压变形直至点火爆炸的全过程;用传感器、电荷放大器、瞬态波形存贮器构成的数据采集系统测量炸药爆炸后形成的空气冲击波超压。试验得到弹丸撞靶速度与相对释放能关系的苏珊感度曲线。本文给出了JO-9159、JOB-9003以及TNT/RDX 40/60等几种炸药的试验结果,对其安全性能作出了初步评价。     

动载荷下固体材料的剪切断裂判据

本文通过对钢质圆柱体弹丸高速撞击平面铝靶板实验的二维数值模拟计算,提出了以最大工程剪切应变作为剪切断裂判据的建议。     

单个弹丸撞击316L不锈钢引起的变形场

运用金属材料表面纳米化试验机对单个弹丸撞击316L不锈钢表面进行了撞击实验;采用激光共聚焦显微镜观察了弹坑的三维形貌,测量不同振动频率下弹坑的直径及离面位移;采用云纹干涉法对弹坑周围的面内应变场进行测量,并分析振动频率及撞击方式对弹坑尺寸、塑性应变大小以及塑性应变区范围的影响;采用有限元方法对单个弹丸垂直撞击试件表面的应变场进行数值模拟,与实验结果进行比较,分析了弹坑周围残余应力的分布。结果表明:随振动频率的增加,弹坑直径和离面位移都增加,频率在50~55Hz,弹坑直径有突变,离面位移和振动频率呈线性关系;振动频率越大,塑性应变越大,塑性应变分布范围均大于弹坑直径的2倍;同一振动频率下弹丸垂直撞击比倾斜撞击的塑性应变大,而塑性应变分布范围相差不大;面内残余应变场的数值模拟结果和实验结果吻合较好,最大误差小于10%。     

单层5A06铝合金板高速撞击实验研究

采用非火药驱动二级轻气炮发射球形弹丸,对单层5A06铝合金板进行高速撞击实验研究,从而模拟空间碎片对航天器防护结构的高速撞击作用。实验得到了该铝合金板在不同的速度区间的损伤模式。结果表明,弹丸撞击速度一定时,弹坑深度和弹坑直径均与弹丸直径呈线性关系。当撞击速度在4km/s至5km/s时,靶板上的弹坑深度和弹坑直径随撞击速度的增大而减小,在其它速度范围内,弹坑深度和弹坑直径随撞击速度的增大而增大。通过固定弹丸直径,变化撞击速度,寻找临界撞击速度的方法获得了该铝合金板在弹丸撞击速度为1.0km/s至4.2km/s时的撞击极限曲线,并将实验弹坑深度与由Cour-Palais方程得到的预测弹坑深度进行了比较,实验弹坑深度大于预测值。     

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

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

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

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