含能微球敏化乳化炸药水下爆炸能量输出特性研究

将内部含有烷烃的含能微球引入乳化基质,得到一种新型乳化炸药。采用水下爆炸实验探究微球质量分数对乳化炸药水下爆炸性能的影响,得到含能微球质量分数为0.2%～7%的乳化炸药水下爆炸冲击波压力-时程曲线。依据压力结果,通过公式计算和分析得到炸药的水下冲击波峰值压力、比气泡能、比冲击波能以及比爆炸能等水下爆炸参数。实验结果表明:含能微球质量分数0.2%的乳化炸药的峰值压力最大,并且随着微球质量分数增大而下降;乳化炸药的比气泡能随着含能微球质量分数的增大先上升再下降,微球质量分数为4%的比气泡能最大;乳化炸药的爆速、比冲击波能以及比爆炸能均随着含能微球质量分数的增大而下降。     

爆炸塔内120kgTNT当量爆炸载荷的数值分析

爆炸塔设计中常采用等当量集中装药估算爆炸塔内部载荷。为研究该方法的合理性,本文采用有限体积方法离散求解积分型Euler方程组;利用瞬时爆轰模型描述炸药爆炸初场,分别计算了乳化炸药和TNT集中装药爆炸对爆炸塔的内部加载情况;给出了塔体及封头内壁面9个特征点的压力时间曲线。结果表明:可以采用集中装药模拟乳化炸药对爆心截面处的超压作用,但不可以采用集中装药模拟乳化炸药对封头中心处的超压作用。因为在封头中心处乳化炸药爆炸会聚产生的超压极值约为TNT集中装药的2倍;塔体柱段的爆心截面处与封头中心处为超压作用最强的部位,两处冲量大致相当;压力最大值出现在封头中心处,由反射冲击波会聚产生。建议塔体建筑设计时对封头中心处和塔体爆心截面处采取相应的防护措施。     

HMX基多组分PBX结构和性能的模拟研究

用分子动力学(MD)方法，对HMX基含少量TATB、F2311(粘结剂)和石蜡(钝感剂)的4组分PBX的结构和性能进行模拟研究。为细致考察各组分对主体炸药的作用，对2组分体系(HMX/TATB, HMX/F2311和HMX/石蜡)也进行类似的MD模拟。为深入揭示钝感机理和钝感剂的作用，还对HMX和石蜡超分子体系的相互作用进行量子化学第一性原理DFT计算。此外，对纯HMX及以它为基的多组分PBX的爆热和爆速进行了理论估算。结果表明，各PBX的弹性较纯HMX的有所改善，以粘结剂组分对主体炸药的力学性能影响最大。各组分的加入均或多或少地降低主体炸药的爆热和爆速。钝感剂与HMX的相互作用很弱，PBX的钝感性不是由电子结构因素所造成。多组分PBX的理论配方设计需综合考虑各种复杂因素。     

一类爆轰合成用乳化炸药的爆轰反应特征

结合热分解特性、爆轰产物状态及数值计算探讨所制备的含Fe、Mn(Zn)元素乳化炸药的爆轰反 应特征。对该系列乳化炸药及其爆轰产物分别进行DSC和DTG实验、XRD检测和TEM 表征,通过比较不 同类型炸药及其爆轰产物的热分解特征、爆轰产物的成分和形貌、晶型畸变度等研究该类炸药爆轰反应特征; 并通过数值计算研究和佐证了不同温度和压强状态下,爆轰固体产物可能存在状态。结果表明:硝酸铵有利 于炸药爆轰时形成相对比较均匀和平稳的爆轰反应结构,而这种结构有利于可重复性合成比较单一和均匀 的爆轰产物。高密度炸药的爆轰产物比低密度炸药的相对纯净。低密度炸药不完全爆轰,导致了爆轰产物杂 质较多;数值计算结果表明,爆轰固体产物分布相图可以辅助分析爆轰固体产物可能存在状态及Mn(Zn)的 爆轰合成机制。     

时序控制预裂爆破参数优化及应用

针对白鹤滩水电站地下厂房的爆破开挖,为降低爆破对岩体的损伤,采用数值模拟与现场试验研究小孔径时序控制预裂爆破中,后爆孔孔间延期时间及后爆孔孔距对时序控制预裂爆破成缝的影响,从而获取合理的延期时间和最佳的后爆孔孔距。研究结果表明,孔径为42mm的时序控制预裂爆破,先爆孔孔距为35cm时,合理的起爆延期时间为75~100μs;综合考虑炸药爆炸能量利用率和成缝效果,最佳的后爆孔孔距为70cm。科学地采用时序控制预裂爆破不仅可以减小钻孔工作量与炸药用量,还能降低对岩体的损伤,可为地下厂房爆破开挖安全、高效的进行提供一条有效的途径。     

钢筋混凝土拱的水下抗爆性能

为探究水下爆炸荷载作用下钢筋混凝土拱的动力响应特性和破坏特征,制作了两个钢筋混凝土拱试件,并开展了水下爆炸试验。试验分为拱外爆炸和拱内爆炸两组,采用10 g乳化炸药,试验时爆源距结构面最小距离为10 cm（起爆点位于拱结构正上方和正下方）,通过传感器记录爆炸试验中钢筋混凝土拱典型断面处的水压力及加速度时程曲线。基于Arbitrary Lagrange-Euler (ALE)算法,建立了空气-水-炸药-钢筋混凝土拱等多介质动态耦合作用模型,将数值模拟结果与试验结果对比,验证了数值方法的可靠性。采用验证后的数值模型进一步研究了拱外及拱内爆炸荷载作用下钢筋混凝土拱的动力响应差异。结果表明：相同炸药当量下,内部爆炸有更多的能量作用于混凝土拱,使结构的动力响应更强烈;外部爆炸下,拱顶、拱腰处产生较大裂缝;内部爆炸时,迎爆面裂缝数量明显增多,拱肩位置出现裂缝。钢筋混凝土拱形结构抵抗外部爆炸荷载的能力明显强于内部爆炸荷载。     

基于通用炸药状态方程分析飞板运动规律的特征线法

从小扰动波(马赫波)的物理概念出发,导出了不依赖流体状态方程表达形式的平面二维超声速定 常流的特征线方程;重新定义了以流体密度为单自变量的Prantl-Meyer函数,形成了求解平面二维超声速定 常流的封闭方程组。还利用这种通用物态方程的特征线差分解法,针对滑移爆轰驱动飞板运动问题构建了爆 轰产物流场内部和飞板边界特征线差分法格式。对TNT炸药和乳化炸药采用JWL状态方程和多方方程进 行了对比计算。结果表明,炸药爆轰对飞板的驱动能力与状态方程表示的炸药的做功能力是一致的。     

采用粉末冶金法制备Ag_2MoO_4/Ag自润滑涂层及其性能研究

由于Ag_2MoO_4在100℃～600℃温度范围内具有良好的润滑性能以及良好的高温稳定性,本文采用高能球磨和粉末冶金技术制备以Ag_2MoO_4和Ag为复合润滑相的试样,并研究了室温时,润滑相含量对试样摩擦磨损性能的影响。研究表明:随着润滑相含量增加,烧结试样孔隙率明显降低;但当润滑相含量增加到15%以后,孔隙率降低不明显。常温时的摩擦实验表明:在润滑相含量(Ag_2MoO_4-Ag)从5%增加到20%时,试样的摩擦系数随之降低;但是当润滑相含量增加到15%后,润滑相含量的增加对摩擦系数影响不明显。同时,随着润滑相含量增加,试样摩擦磨损形貌更优,表面磨屑、犁沟更少,更有利于润滑膜形成。试样中润滑相(Ag2MoO4-Ag)含量为15%时其摩擦性能最佳。     

水下爆炸冲击波载荷作用下冰层破碎特性及其影响因素

水下爆炸破冰是复杂的爆破工程，为了研究冰层在水下爆炸冲击波载荷作用下的破碎特性及规律，利用几何动力分析软件LS-DYNA对水下爆炸破冰的过程进行数值模拟，并将计算结果与实验结果进行对比，误差在8%以内，验证了数值模型的有效性。根据本文中的建模方法及建立的模型，计算不同的实验工况:实验场地环境不变，调整爆距分析不同爆距下冰层破碎特性；调整药量、爆距和冰厚，通过正交设计方法设计9组实验方案，应用灰色系统理论对3种因素进行分析，建立了各个因素与破冰半径之间的灰色关联度系数及灰色关联度。分析结果表明:药量为100 g，冰厚为29 cm，水深为2.9 m，爆距范围为0.3～1.5 m破冰的半径范围为0～1.1 m，最佳爆距范围在0.3～0.45 m之间；根据以上9种工况的分析可知，药量(100、200、300 g)、爆距(0.3、0.6、0.9 m)和冰厚(24、28、32 cm)对破冰半径的影响的主次关系依次为爆距、药量、冰厚。     

常温下JBO-9021高能钝感炸药直径效应实验

通过光电联合测试方法,测量了常温下直径为10、15和30mm 高能钝感炸药柱JBO-9021炸药的 定态爆速和波形,根据拟合的爆轰波阵面形状,分析了波阵面法向速度Dn与当地曲率 之间的关系。研究结 果表明,JBO-9021炸药的定态爆速随药柱直径的增大而增大,JBO-9021(w(HMX)=15%)相对JB-9014炸 药(w(TATB)=95%),爆轰波形更平坦,爆轰波非理想传播行为得到一定改进。     

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

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

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

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