舱内爆炸载荷及舱室板架结构的失效模式分析

通过对典型半穿甲导弹打靶实验中舰艇结构破坏模式的观察，结合数值模拟，分析了舱内爆炸载荷的特征以及舱内爆炸下舱室板架结构的失效模式。结果表明，舱内爆炸下，舱室板架结构承受的冲击载荷及失效模式与敞开环境爆炸下加筋板结构承受的冲击载荷及失效模式有较大区别，其动态响应难以用敞开环境爆炸下加筋板结构的动态响应描述；舱内爆炸载荷主要有壁面反射冲击波、角隅汇聚冲击波以及准静态气体压力，其中两壁面和三壁面角隅汇聚冲击波的强度分别为相同部位壁面反射冲击波强度的5倍和12倍以上；舱室板架结构主要有4种失效模式，其中模式Ⅲ、Ⅳ较常发生；舱室板架结构加强筋布置在迎爆面将使板架中部的局部破坏程度增加，但有利于削弱角隅汇聚冲击波强度，减小板架沿角隅部位的撕裂破坏。     

原型和模型不同材料时加筋板冲击动态响应的相似预报方法

:根据 定律提出了加筋板结构在水下爆炸冲击波作用下模型和原型动态响应相似的条件。在此 基础上,考虑到加筋板模型制造加工过程中的工艺和成本等因素,提出了利用普通钢加筋板模型预报船体钢 加筋板原型在水下爆炸冲击波作用下动态响应的相似预报方法。该方法考虑了流固耦合和材料应变率强化 效应等因素的影响,可用于指导利用常用的普通钢模型试验预报实船局部板架在水下爆炸冲击波作用下的 动态响应。     

浅层水中爆炸冲击波对混凝土墩斜碰撞作用试验研究

通过含铝炸药JHL-3的实爆试验,得到混凝土墩在单个装药浅层水中爆炸、两个装药浅层水中对称及不对称设置同步起爆爆炸作用下、混凝土墩迎爆面上的冲击波压力响应数据;获得了迎爆面中心反射压力峰值计算模型;分析了两个装药浅层水中爆炸冲击波对混凝土墩绕射及透射作用效应。     

水下成组药包爆炸作用下冰盖动态响应数值模拟

利用大型有限元软件ANSYS/LS-DYNA建立了成组药包水下爆炸冰盖动态响应模型,通过数值计算,得到了冲击波峰值压力变化规律,峰值压力与理论计算结果基本吻合。本文分析了冰盖在成组药包水下爆炸冲击载荷下的应力分布及垂直位移响应特征。结果表明,冰盖迎爆面为压缩破坏,背爆面为拉伸破坏,两者数值均十分接近,其中迎爆面最大压力值可达18.12MPa,冰盖近爆炸点最大垂直位移为1.211cm,两药包连线中点最小位移为0.15cm,达到脆性冰盖形成贯通裂隙的基本条件,从而确定了冰盖在水下成组药包大间距布设条件下,其动态破坏形式是以产生裂隙为主要特征。     

接触爆炸作用下高聚物复合板毁伤特性分析

高聚物材料具有成型快和膨胀性能好的特点，该材料与碎石和钢筋的复合结构应用于地基处理和城市道路脱空除险加固，具有明显优势。本文中，设计并制作了高聚物碎石板和钢筋高聚物板，开展了接触爆炸冲击下的试验研究，通过毁伤尺寸和所测冲击波数据探讨了2种板的毁伤特性。基于ANSYS/AUTODYN非线性显式有限元程序，建立了试验中毁伤更严重的钢筋高聚物板的接触爆炸全耦合模型，并通过与试验结果的对比，验证了所建耦合模型的准确性和适用性。参数化分析了钢筋高聚物板对炸药量和板厚的敏感性，并利用多参数非线性回归分析方法，提出了钢筋高聚物板迎爆面和背爆面破坏直径的预测公式。结果表明：接触爆炸作用下，高聚物碎石板的毁伤模式以接触部位的局部震塌冲切破坏为主，除此之还有一些毁伤裂纹；钢筋高聚物板的破坏模式主要是迎爆面爆坑毁伤、背爆面剥落损伤和中心冲切贯穿破坏。高聚物碎石板和钢筋高聚物板对爆炸冲击波都具有良好的衰减作用，有望应用到抗爆炸冲击防护领域。     

水下冲击波载荷作用下气背固支圆板动态毁伤实验

为了研究水下近爆载荷作用下舰艇水下结构的动态变形及失效毁伤模式，利用水下爆炸冲击波等效加载装置结合高速摄影技术，对两种厚度的气背固支5A06铝合金圆板进行了水下冲击波加载实验。得到了气背固支圆板塑性大变形、中心拉伸撕裂和边界剪切破坏3种典型失效模式的动态响应历程。比较分析了冲击波强度、冲击因子、损伤参数和响应参数4种毁伤判据对该类靶板毁伤模式的判别能力。实验结果表明:考虑了结构因素的损伤参数和响应参数能够更为全面的判别结构的失效毁伤情况。     

预制孔靶板在爆炸冲击波载荷作用下的动态响应

对四边固定带预制孔靶板在爆炸冲击波载荷作用下的动态响应问题进行了实验研究。通过设计 爆炸冲击波载荷对2种材料不同孔数、不同孔径的靶板的冲击工况,获得了各靶板的中心点挠度。通过分析, 得到了预制孔靶板中心点挠度与预制孔数、孔径的相关规律,结合无孔平板在爆炸冲击波载荷下的挠度公式, 建立了广义的靶板中心点挠度公式,可用于计算不同孔靶板的中心点挠度。研究结果可作为平板材料在破片 冲击波联合毁伤研究的基础,为杀爆战斗部威力评估提供一些依据。     

长持时爆炸冲击波荷载作用下梁板组合结构的动力响应

为研究钢筋混凝土梁板组合结构在长持时远爆冲击波荷载作用下的动力响应及毁伤形态,通过实验获得了梁板组合结构的破坏形态和背爆面中心点位移变化。利用有限元软件对钢筋混凝土梁板组合结构的动态响应过程进行数值模拟研究,模拟得到的结构破坏现象与实验吻合较好。在此基础上,分析了梁板组合结构在相同冲量、不同峰值爆炸荷载作用下组合结构的动态响应和破坏过程,并结合挠跨比与破坏形态划分破坏模式。研究结果表明,相同冲量作用下,随着爆炸荷载峰值强度增加,梁板组合构件的破坏程度逐渐增加,破坏模式从弯曲破坏向弯剪联合破坏转换,最后呈现冲切破坏模式;组合构件中板部分发生破坏的时间早于交叉梁部分、破坏程度大于交叉梁。     

水下爆炸对沉箱重力式码头毁伤效应

通过开展不同爆距下水下爆炸对沉箱重力式码头模型毁伤效应试验，对水下荷载进行了采集分析，对毁伤因素、毁伤模式和毁伤机理开展了研究，初步讨论了爆距的影响。结果表明:试验中未形成完整的气泡脉动过程，荷载超压主要出现在冲击波传播阶段，爆炸冲击波、水底反射波和侧壁反射波是主要的毁伤因素；水下爆炸对沉箱重力式码头造成的毁伤程度大、模式多、机理复杂，主要毁伤部位为迎爆面外墙、迎爆侧管沟、封仓板和面板；爆距越近码头毁伤越严重；当爆距过近时，爆炸能量被迎爆面结构变形大量吸收，迎爆面毁伤程度大幅增大，码头其他部位毁伤程度增幅放缓。     

水下近距/接触爆炸加载下圆柱壳结构动态响应行为试验研究

为加深水下近距/接触爆炸加载下圆柱壳结构动态响应行为认识，设计典型圆柱壳结构模型，开展了水下近距/接触爆炸加载下圆柱壳结构动态响应光电联合测试，获得了冲击波、气泡与圆柱壳结构相互作用高速光学物理图像、动态应变、超压载荷、毁伤模式等试验数据。通过高速光学物理图像和三维激光扫描毁伤形态的分析，给出了冲击波、气泡与圆柱壳结构相互作用物理过程及最终毁伤模式；通过动态应变的分析，给出了圆柱壳结构迎爆面和背爆面在加载过程中应变拉伸压缩转变和响应阶段的划分；通过超压载荷的分析，明确了装药爆轰完全性以及接触爆炸加载下结构吸能对超压的影响。研究表明：爆距的变化会显著影响圆柱壳结构的毁伤形态，近距加载下圆柱壳结构主要呈现塑性大变形，接触加载下圆柱壳结构主要呈现撕裂破坏；近距加载下圆柱壳结构迎爆面空化区的形成及溃灭形成的二次加载毁伤效应不容忽视，值得深入研究；研究成果可为水下近距/接触爆炸加载下圆柱壳结构毁伤评估提供参考和依据。     

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

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

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

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