大直径SHPB装置的数值模拟及实验误差分析

介绍了用于混凝土冲击试验的Φ100SHPB装置。利用数值模拟和应变测试相结合的方法讨论了应力波在该装置中的传播特性。二维效应导致应力波在大直径SHPB装置传播过程的弥散现象。同时在实验的大部分时间内大尺寸试样的应力(应变)存在不均匀性。因此大直径SHPB装置进行混凝土冲击压缩应力-应变曲线的实验研究会产生误差。利用数值模拟分析了误差影响程度,同时验证了用混凝土试样应变和透射杆应力得到应力-应变曲线初始段的可行性。     

用HJC本构模型模拟混凝土SHPB实验

运用Holmquist-Johnson-Cook (HJC) 本构模型对混凝土的SHPB实验进行了数值模拟。解决了罚函数算法中罚因子合理数值的选取问题。利用模拟结果按SHPB两波法重构了试样的应力应变曲线。分析了混凝土材料的SHPB实验得到应力应变曲线的有效段范围和各段的力学规律。通过比较实际混凝土材料SHPB实验和数值模拟得到的应力应变曲线,发现两者体现的力学行为很相似,即HJC模型是描述该类材料的一种合理本构模型。模拟了试样不同平行度公差下的SHPB实验,发现在一定应变率范围内其影响程度远大于试样应力(应变)不均匀性。     

Al2O3陶瓷材料应变率相关的动态本构关系研究

采用改进的SHPB实验方法对Al2O3陶瓷的动态力学性能进行了研究,得到了材料在较高应变率范围内的动态应力应变曲线。结果表明,Al2O3陶瓷为弹脆性材料,其动态应力应变呈非线性关系,在较高的应变率范围内,陶瓷材料的动态应力应变关系是应变率相关的;材料的初始弹性模量、破坏应力、破坏应变值随应变率的增大而增大。基于损伤力学的基本理论,给出了Al2O3陶瓷的一维损伤型线性弹脆性本构模型。根据SHPB实验结果确定模型中的参数,得到了Al2O3陶瓷应变率相关的损伤型动态本构方程。     

高应变率下多孔未极化PZT95/5铁电陶瓷的非线性力学行为

采用添加造孔剂的方法制备了4种不同孔隙率的未极化PZT95/5铁电陶瓷。采用基于超高速相机与数字图像相关性方法的试样全场应变测量技术以及分离式霍普金森压杆（SHPB）技术,对多孔未极化PZT95/5铁电陶瓷进行高应变率单轴压缩实验研究。全场应变测量结果显示:轴向应变仅在试样中部分布较均匀,将该区域的平均应变作为应力-应变关系中的试样应变测量值较为合理,而由SHPB原理计算的试样应变值明显偏大,需要摒弃或修正传统的SHPB数据处理方法。通过波形整形技术实现了恒应变率加载,弱化了径向惯性效应的影响,揭示出多孔未极化PZT95/5铁电陶瓷的压缩强度具有显著的应变率效应。通过分析试样轴向应变和径向应变随着加载应力的变化,阐明多孔未极化PZT95/5铁电陶瓷的非线性变形行为的物理机制是畴变和相变共同作用,并发现畴变临界应力和相变临界应力都随着应变率升高而增大。保持加载应变率不变,讨论了孔隙率对多孔未极化PZT95/5铁电陶瓷动态力学行为的影响,发现随着孔隙率的升高,动态压缩强度呈非线性衰减,而畴变临界应力和相变临界应力则基本呈线性衰减。     

PBX炸药本构关系的实验研究

在SHPB装置上开展了某典型PBX炸药的单轴压缩、间接拉伸实验。利用入射波整形技术,结合 半导体应变片和石英晶体联合测试的方法,实现了PBX炸药的应力平衡和恒应变率加载,得到了不同应变率 下某PBX炸药的压缩应力应变曲线,初步建立了该炸药的修正Sargin唯象模型;开展了3种PBX炸药的动 态巴西实验,通过高速摄影获得了炸药的破坏过程,结合数字相关技术,获得了试样裂纹附近的应变场分布, 初步建立了描述3种炸药动态拉伸行为的修正Johnson-Cook模型。模型曲线与实验结果符合较好。     

非均质结构PBX炸药的动态压缩过程模拟

采用离散元方法,构建了考虑PBX炸药晶体颗粒与黏结剂非均匀性的计算模型,通过在炸药试样两端采用相向飞片加载、将原SHPB实验的应力边界替代为速度边界的方法,开展了PBX炸药的动态压缩数值模拟研究,再现了考虑PBX炸药非均质特性的SHPB实验过程,获得了不同加-卸载路径下炸药应力应变曲线以及相应的损伤破坏图像。模拟结果表明:在动态压缩过程中,虽然PBX炸药处于整体应力平衡,由于PBX高度非均质性,内部应力分布并不均匀,晶体间应力以应力桥形式传递;应力曲线为试样整体平均,试样内局部所受真实应力可能高于曲线峰值应力;出现卸载回滞曲线的试样发生局部损伤破坏,而完全软化曲线的试样发生整体失稳破坏。     

用于高强度材料的SHPB实验添加垫块法

提出了用于高强度材料的改进的SHPB实验方法添加垫块法,运用数值模拟方法,利用有限元程序LS-DYNA3D分析了添加垫块实验方法的合理性和可行性。根据一维应力波理论,给出了数据处理的修正方法。作为应用实例,采用改进的实验方法对高强度的Al2O3陶瓷材料的动态力学性能进行了研究,得到了比常规方法较高的应变率及应力应变范围的动态应力应变曲线,表明Al2O3陶瓷为应变率相关的非线性弹脆性材料。结果表明,添加垫块实验方法可有效地防止实验中压杆端面的变形,提高试件的应力应变及应变率水平。添加垫块实验方法为在SHPB装置上实现高强度材料的动态实验提供了一种方便实用的途径。     

混凝土SHPB束杆装置冲击实验的数值仿真

利用Ansys/Ls-dyna有限元软件模拟不同打击气压下混凝土SHPB束杆实验,对实验装置模型进行了简化;采用梯形应力脉冲载荷,模拟得到了实验中试件的应力-应变曲线。结果显示:模拟结果与实验结果吻合较好,验证了该方案的可行性;同时试样应力峰值随打击气压增大而增大。对不同幅值梯形应力脉冲加载下的SHPB束杆装置和单杆装置实验进行数值仿真,研究表明:束杆装置中试样应力峰值比单杆装置小一些;随应力脉冲载荷幅值增大,束杆装置中试样应力峰值增长速度要高于单杆装置。     

花岗岩动态力学性能的实验研究

介绍了利用Φ74mm SHPB试验装置进行花岗岩的动态压缩试验和动态劈裂拉伸试验的结果。通过试验获得了花岗岩在100s-1～102s-1应变率加载时的动态力学参数。同时,对试件内的动态应力分布进行了数值模拟,验证了动态试验的有效性。试验结果表明,花岗岩的动态压缩和动态劈裂拉伸的力学性能表现出显著的应变率效应。随着应变率的增加,其冲击压缩破坏应力、冲击压缩破坏应变、弹性模量、冲击压缩应变能密度和动态劈裂拉伸破坏强度均有一定程度的增加。     

改进SPH方法在陶瓷材料层裂数值模拟中的应用

为了探讨铝飞片撞击陶瓷材料时的层裂现象,采用改进SPH方法模拟应力波在陶瓷材料中的传播。结果表明,当离散粒子分布不均匀时,数值模拟计算的自由面速度时程曲线与实测曲线吻合良好。对比CSPM方法,改进SPH方法的精度更高。提出适用于数值模拟的陶瓷材料损伤演化方程,对脉冲载荷下陶瓷/钢层合板层裂的破坏过程进行数值模拟,结果表明,由于陶瓷的波阻抗高于钢的,且抗压强度远高于抗拉强度,因此拉应力引起的层裂破坏是主要的。即使在材料内部传播的只是弹性压缩波,当弹性波到达材料界面时,由界面反射引起的卸载波也能导致陶瓷发生层裂破坏。     

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

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

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

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