混凝土SHPB试验数值模拟研究

基于混凝土HJC动态本构模型,分析了该模型方程及参数获得方法,以混凝土材料的基本力学参数确定了HJC模型中各参数值,并运用LS-DYNA软件完成混凝土SHPB试验的三维数值模拟.对比分析数值模拟结果重构所得试样的应力应变曲线与实际试验结果的相似性,研究表明,数值模拟结果可再现试验中混凝土的动态力学行为与破坏特征,而且,数值模拟可很好的反映应力波在波导杆与试件中传播过程,符合一维弹性波理论规律.     

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

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

混凝土三维细观模型的建模方法与力学特性分析

根据混凝土材料的细观组成和结构特点，基于三维Voronoi图形提出了一种简单高效的混凝土细观模型生成方法，利用塑性损伤模型对该细观模型进行了单、多轴应力状态下的准静态分析以及SHPB动态有限元分析。结果表明，数值模拟得到的应力应变曲线和破坏模式与实验结果基本吻合，本文中提出的混凝土三维细观模型可较好地模拟混凝土的静、动态力学特性，为进一步从细观力学角度研究混凝土损伤演化规律和破坏机理提供了模型基础。     

吸能包装模型结构的冲击响应

利用空气炮冲击实验对吸能包装结构的跌落过程进行模拟，进行了缩比模型的正撞和30°斜撞实验，针对模型实验进行了数值分析，获得了吸能包装结构模型在撞击过程中的应力分布和塑性变形，并将计算情况与实验结果进行了分析。结果表明:在撞击中吸能包装结构主要通过缓冲木材的塑性变形及外钢壳屈曲产生的塑性铰吸收能量，塑性变形主要集中于撞击端，而远离撞击端未见塑性变形；计算中木材本构参数采用顺纹方向压缩应力应变曲线具有一定的有效性。     

单轴循环载荷下镁合金宏观塑性本构关系研究

基于微观织构演化,本文对镁合金的宏观塑性本构模型进行了改进,得到了室温下镁合金AZ31B和AZ61A单轴循环拉压载荷作用下的应力应变曲线。模型中区分了滑移、孪晶与去孪三种不同微观塑性变形机制,并且研究了它们对材料宏观循环塑性行为的影响。通过引入初始非零背应力来模拟循环应力应变曲线中拉压不对称性,给出了材料在单轴循环拉压载荷下背应力及屈服面的演化规律,并运用线性预测塑性回拉技术与隐式积分法求解相应的非线性方程组,编写了与.ABAQUS接口的UMAT用户子程序。预测与实验结果的比较显示该模型能较好地模拟镁合金材料在单轴循环拉压载荷作用下的力学行为,具有较高的准确性与良好的实用性。     

不同应变率下聚乙烯材料的压缩力学性能

为了研究聚乙烯材料在不同应变率下的压缩力学性能，通过准静态实验和动态实验获得聚乙烯材料不同应变率下的应力应变曲线，分析发现:聚乙烯的弹性模量和屈服强度随应变率增大而增大，具有明显的黏弹塑性；聚乙烯材料进入塑性阶段，其应力应变曲线在不同应变率下具有相近的变化趋势，即塑性切向模量近似相同。根据聚乙烯材料的压缩力学性能，建立了弹性区、屈服点和塑性区的分段本构模型。该模型的屈服点和塑性段与实验结果吻合较好，由于弹性段采用线弹性模型，与实验结果存在一定偏差，可近似描述材料的弹性行为。     

45钢的J-C损伤失效参量研究

为了在结构碰撞效应的有限元分析中描述材料行为,通过开展45钢在不同应力状态和温度下的准静态材料力学性能实验及拉伸SHB实验,考察了应力状态三轴度、温度和应变率对材料失效应变的影响。由实验数据得到了Johnson-Cook失效模型参量,并通过出现失效的Taylor撞击实验和数值模拟进行了一定的验证,表明模型描述与实验结果的趋势一致。     

考虑塑性应变记忆恢复的循环棘轮本构模型研究

在统一粘塑性循环本构理论框架下,以Ohno-Abdel-Karim非线性随动硬化模型为基础,建立了一个循环本构模型。模型通过引入塑性应变幅值记忆效应,并在塑性应变记忆项中加入恢复系数,提高了对循环硬化材料单轴棘轮行为的预言能力。将模型应用于316L不锈钢单轴棘轮行为的描述中,模拟不同平均应力、应力幅值下的棘轮应变,均与实验数据吻合较好,证明本文改进的本构模型能合理地描述循环硬化材料的单轴棘轮行为。     

冲击载荷下脆性空心颗粒破碎机理

为考察脆性空心颗粒在冲击载荷作用下的应变率效应和破碎行为的细观机理,以粉煤灰漂珠为研究对象,基于低速冲击实验和有限元数值模拟,对比了典型空心颗粒材料在不同加载速率下的力学响应特性和细观压溃行为,阐释了材料宏观应变率效应产生的细观机理,获得以下结果。（1）在0.001～300 s?1应变率范围,漂珠颗粒的破碎率和Hardin破碎势平均提升了约21%和10%～30%,材料比吸能提升了50%～125%,比吸能的额外增加主要与动态颗粒滑移产生的摩擦耗能相关。颗粒平均尺寸较大的试样体现出更强的应变率效应。（2）初始压溃阶段的应力应变响应特征的数值模拟结果与实验结果较吻合,低速冲击下动态二次压溃现象产生的细观机理为动态颗粒滑移和压紧行为对加载速率的依赖性。(3) 数值模拟表明,冲击加载下产生相同应变时颗粒的损伤程度和范围大于准静态加载,这与实验所得破碎势随应变率增加的结果一致。对比低速冲击实验的相对破碎势分析和细观数值模拟结果可知,脆性颗粒堆积材料在动态冲击下表现出的宏观应变率效应主要归因于颗粒压溃行为的率敏感性和动态加载下颗粒破碎能量利用率的降低。     

一个综合模糊裂纹和损伤的混凝土应变软化本构模型

本文研究就变软化材料的本构关系，提出了一个考虑损伤的粘塑性模型，损伤不仅影响材料的临界应力，而且影响材料的粘塑性，为模拟材料的应变软化行为，假设受损混凝土的破坏局部区域由模糊裂纹和损伤所统治，软化模量和局部区域尺度参量依赖于模糊裂纹扩展时释放的断裂能的参变量，用文中提出的模型计算了混凝土单轴压缩时不同应变速率下的瞬时应力应变响应以及等应力长期作用下的徐变，均得到很的结果。     

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

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

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

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