45钢动态塑性本构参量与验证

运用静态实验机和SHPB装置,对45钢在常温～750C、应变率为10-4～103 s-1下的力学行为进行了研究,拟合得到了Johnson-Cook本构模型参量。开展了Taylor圆柱撞击的火炮实验,运用LS-DYNA进行了相应的数值模拟,通过二者的比较对本构模型参量进行了验证,表明所得模型参量可以较好地描述材料在高速变形下的大应变力学行为。     

高温下混凝土的本构模拟及破坏分析

针对已建立的高温下混凝土中化学-热-水力-力学耦合过程分析的分级数学模型,发展了混凝土的化学-热-水力-力学(CTHM)耦合本构模型。在已有的Willam-Warnke弹塑性屈服准则基础上发展了考虑脱水和脱盐引起的材料损伤及化学塑性软化、塑性应变硬化/软化和吸力硬化的广义Willam-Warnke本构模型,模拟高温下混凝土的材料非线性行为。为保证全局守恒方程的Newton迭代过程的二阶收敛率,导出了非线性化学-热-水力-力学(CTHM)耦合本构模型的一致性切线模量矩阵。数值结果显示了本文所发展的化学-热-水力-力学(CTHM)耦合本构模型在模拟高温下混凝土中复杂破坏过程的能力和有效性。     

5A06铝合金的动态本构关系实验

运用材料试验机和分离式霍普金森压杆装置(SHPB)对3种不同加工及热处理状态的5A06铝合金在常温~500 C、应变率为10-4~103 s-1 下的力学行为进行了实验研究。基于Johnson-Cook (JC) 本构模型,通过实验数据拟合得到了每种状态下材料的本构模型参量。对Johnson-Cook本构模型中的应变率强化项作了修正,修正后的Johnson-Cook本构模型与实验数据基本吻合,从而确立了3种状态下5A06铝合金的动态本构关系。     

宽应变率范围下2A16-T4铝合金动态力学性能

为了研究2A16-T4铝合金的动态力学性能,利用电子万能试验机、高速液压伺服试验机及霍普金森压杆(SHPB)装置进行常温下准静态、中应变率和高应变率的动态力学性能实验,得到不同应变率下的应力应变曲线,基于修正的Johnson-Cook本构模型对它进行拟合,并分析材料中应变率力学特性对模型应变率敏感参量的影响。结果表明:2A16-T4铝合金在应变率10-4~102 s-1范围内应变率敏感性较弱,而在102~103 s-1范围内应变率敏感性较强,且应变率强化效应随塑性应变的增大而减小;同时,在10-4~103 s-1范围内具有较强的应变硬化效应,且应变硬化效应随应变率的增大而减小;此外,修正Johnson-Cook本构模型的拟合结果与实验结果吻合很好,能够很好表征材料的动态力学行为,且考虑材料中应变率力学特性可提高本构模型参量的准确性。     

高温下混凝土化学塑性-损伤耦合本构模拟及破坏分析

提出了一个用于模拟高温下混凝土中化学-热-湿-力学耦合行为的化学塑性-损伤耦合本构模型.发展了用于积分率形式的耦合本构方程的三步算子分裂算法.导出了化学-热-湿-力学(CTHM)耦合本构模型的一致性切线模量矩阵,以保证对于全局耦合控制方程Newton迭代过程的二阶收敛率.算例数值结果显示了本文发展的化学塑性-损伤耦合本构模型在重现火灾和热辐射条件下的混凝土中化学-热-湿-力学耦合行为的有效性.     

基于Khan-Huang模型高强度钢DP600率相关特性实验与本构模型研究

对汽车用高强钢DP600在10-4s-1~103s-1应变率范围内进行力学拉伸实验,在此基础上,采用唯象的方法对Khan-Huang本构模型进行修正建立DP600考虑率敏感效应的本构模型,对高应变率下试样的拉伸过程进行数值模拟验证模型有效性.结果表明,高强钢常温下具有明显的应变率敏感特性,高应变率下材料的屈服应力接近低应变率下的两倍.对Khan-Huang本构模型进行修正后可以比较准确地描述材料在不同应变率下力学行为,且参数容易确定,便于在有限元软件中实现接口,因此该模型可进一步应用于汽车的碰撞安全数值仿真中.     

强动载荷下金属材料塑性变形本构模型评述

爆炸、高速冲击等强动载荷作用下金属材料塑性流动的本构模型研究是冲击动力学领域备受关注的课题,近几十年来的研究成果极大地推动了冲击波物理、冲击动力学、材料科学等领域的发展. 本文对金属材料塑性变形动力学数值计算中重要的、得到普遍应用的本构模型进行了回顾和评述,进一步明确其适用范围、优点和不足,为冲击载荷下材料响应的数值模拟对本构模型的选取提供参考. 其中最新的、适用范围最广的PTW本构模型应给以特别地关注,其研究方法及应用前景都值得我们借鉴和期待.      

考虑温度影响的混凝土弹塑性本构模型及其应用

基于弹塑性理论的框架,引入热力学中的自由能概念,建立了考虑温度影响的混凝土弹塑性本构模型。假设混凝土为弹塑性解耦材料,通过采用包含温度的自由能函数建立了弹性应力-应变关系的表达式;基于Drucker-Prager强度准则构造屈服函数,以塑性剪应变作为硬化参量,通过在硬化规律中引入温度变量来反映温度对混凝土塑性性质的影响,得到了塑性应力-应变关系的表达式。基于ABAQUS有限元软件平台,应用半隐式回映算法编制了上述弹塑性本构模型的UMAT子程序,完成了该模型的程序实现。通过对C20及C40混凝土单元在不同温度条件下单轴压缩与双轴压缩的数值模拟,并将模拟结果与已有文献试验结果进行了对比,二者符合较好,验证了该模型的正确性。应用该模型模拟混凝土材料的力学特性,开展了不同温度下集中荷载作用的混凝土简支梁及混凝土轴心受压柱力学响应的数值模拟。计算结果表明,随着温度的升高,相同荷载作用下混凝土构件的力学响应逐渐增大,反映了混凝土材料的力学性能随着温度的升高逐渐退化,这符合不同温度条件下混凝土构件力学响应的基本规律。     

SS304不锈钢高温非比例多轴循环变形行为的粘塑性本构描述

在350℃和700℃下SS304不锈钢的非比例多轴循环变形行为进行了系统的实验研究。在此基础上,在统一粘塑性本构理论的框架下改进和发展了一个新的循环本构模型。该模型给出了新的背应力演化方程,引入了非比例度参量,考虑了温度效应和最大塑性应变幅值记忆效应,能够对材料的高温非比例多轴应变循环变形行为和棘轮行为进行统一描述。模型的模拟结果与实验结果比较表明:该模型对SS304不锈钢高温非比例多轴循环变形行为的描述比较合理。     

混凝土材料冲击特性的研究

基于混凝土材料强冲击加载下的试验研究，提出了两种损伤型动态本构模型： 损伤型黏弹性本构模型和损伤与塑性耦合的本构模型. 通过模型计算结果与冲击试验结果的 比较可发现，随着冲击速度的提高，混凝土材料内部产生了显著的塑性变形，由此 损伤型黏弹性本构模型的应用就存在一些不足. 而损伤与塑性耦合的本构模型由于考虑了裂 纹扩展引起的材料强度和刚度的弱化，以及微空洞缺陷塌陷引起的塑性变形，因而能更好地 用于模拟强冲击载荷作用下混凝土材料的冲击响应特性.     

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

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

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

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