混凝土本构模型及其动态有限元算法研究

混凝土的力学性质及其本构研究是当前力学领域的一个重要课题。混凝土的实验数据表明其力学性质表现出很强的各向异性及并且在高加载速率的情况下，混凝土的性质物美价廉 时时明显不同。这就要求我们采用应变率相关的弹粘弹性本构模型来描述混凝土的力学行为。本语文从Ｏｔｔｏｓｅｎ的四参数混凝土模型出发，引进损伤和应变率的影响，参考关联塑性的理论，改进了混凝土的本构模型。为了把混凝土的数值模拟推向实用，本文概括混凝土     

基于统一强度理论的岩石弹塑性本构模型及其数值实现

基于弹塑性力学理论，以统一强度准则为屈服准则，建立了考虑硬化/软化行为和应变率效应的岩石弹塑性本构模型；采用Fortran语言通过LS-DYNA的用户自定义材料接口(Umat)对该弹塑性本构模型进行编程，并把该程序生成求解器以达到对该模型进行应用的目的；通过岩石的单轴压缩实验和SHPB实验对所建的弹塑性本构模型进行验证，结果表明，该弹塑性本构模型能够反映岩石在准静态和动态下的力学行为。     

描述大应变率范围下材料响应的粘塑性本构模型

以位错动力学理论中的Ｏｒｗａｎ和Ｇｉｌｍａｎ关系为基础建立描述率相关材料非弹性响应的基本方程，选择材料准静态实验的单轴响应作为强化演化的规律，并考虑应变率敏感程度随变形产生变化的特性，建立了适用于大应变率范围内率相关材料的统一型粘塑性本构模型。对铝１１００－０在应变率范围１０－５～１０４ｓ－１内产生的有限塑性应变的单轴响应进行了理论预测，与Ｋｈａｎ和Ｈｕａｎｇ［１］的实验数据及模型预测结果进行了比较，结果表明本文模型具有较高的预测精度，在高应变率和较大应变下不容忽视率敏感参数随变形的变化。     

Comp.B复合炸药动态压缩力学性能和本构关系的研究

本文利用自制的含能材料变温动态压缩实验装置，采用准静态应变速率（１０－４／ｓ）和中等应变速率（３／ｓ），对国产复合炸药Ｃｏｍｐ．Ｂ进行了动态压缩实验。测试了Ｃｏｍｐ．Ｂ在不同温度、不同应变速率条件下的杨氏模量、断裂强度等力学性能参数。实验结果证明，Ｃｏｍｐ．Ｂ具有明显的应变率相关和温度效应。采用Ｊｏｈｎｓｏｎ提出的温度、应变率相关的本构模型，根据实验数据拟合了Ｃｏｍｐ．Ｂ材料的本构关系，分析表明该本构模型可以很好地描述材料的应变率和温度效应。这些基础研究为含能材料动态力学性能的研究和炸药早爆机理的理论分析提供了依据。     

热粘塑性体的积分—微分型本构关系

应用ИЛЪЮШИН关于应力是五维偏应变空间变形历史的泛函的概念和Ｖａｌａｎｉｓ有关内时理论的描述，本文提出，对热粘塑性体，应力可设为应变、应变率和温度历史泛；并应用Ｍｉｌｌｅｒ和其它一些作者有关内变量演化方程的描述，由此建立了热粘塑性体的积分－微分型本构方程这一积分－微分型本构关系大体和Ｍｉｌｌｅｒ微分型模型等价。对１０２０钢的单轴本构响应进行了数值模拟，和Ｔａｎａｋａ与Ｍｉｌｌｅｒ的分析及一些实     

传爆药静态压缩力学性能及损伤特性研究

为研究聚黑-14C（JH-14C）传爆药静态压缩力学性能及损伤特性,开展准静态压缩实验,获得了不同应变率下的应力-应变曲线,建立了描述不同应变率下JH-14C力学行为的非线性本构模型;利用扫描电镜（SEM）对回收试样进行细微形貌观测,获得了准静态压缩JH-14C损伤特性的表征。结果表明:JH-14C压缩强度随应变率的升高而提高;实验与计算结果对照验证了本构模型的有效性;准静态压缩实验中,JH-14C主要损伤模式为脱湿和穿晶断裂。     

基于裂纹扩展脆性岩石动力压缩力学特性研究

动态压缩荷载作用下，脆性岩石内部动态细观裂纹扩展特性，对岩石宏观动态力学特性有着重要的影响。然而，对岩石内部动态细观裂纹扩展与宏观动态力学特性的关系研究较少。基于准静态裂纹扩展作用下的应力-应变本构模型、准静态与动态裂纹扩展断裂韧度关系、裂纹速率与应变率关系模型及应变率与动态断裂韧度关系，提出了一种基于细观力学的动态应力-应变本构模型。其中裂纹速率与应变率关系，是根据裂纹长度与应变关系的时间导数推出；应变率与动态断裂韧度关系，是根据推出的裂纹速率及应变率关系，与裂纹速率及断裂韧度关系相结合而得到。研究了应变率对应力-应变本构关系及动态压缩强度影响。并通过试验结果验证了模型的合理性。讨论了岩石初始损伤、围压、模型中参数m、ε₀和R对应力-应变关系、动态压缩强度和动态弹性模量的影响。研究结果可为动态压缩荷载作用下深部地下工程脆性围岩稳定性分析提供了一定的理论支持。     

一种低应变率下钢纤维混凝土本构关系

采用岩石力学试验机,开展了同一低应变率下,钢纤维体积含量分别为0%、0.75%、1.5%、3%、4.5%、6%的混凝土材料的准静态力学试验,并获得了材料的应力-应变曲线,由试验结果可知:钢纤维混凝土材料的峰值应力、峰值应变随着纤维含量的增加而提高,材料达到峰值应力后呈现应力的应变软化特征。基于试验应力-应变曲线,提出了一种合理的钢纤维混凝土准静态本构关系,通过拟合试验应力-应变曲线,获得了相应的材料参数,拟合结果表明采用该形式的本构关系方程能很好地拟合试验曲线。     

热粘塑性体的积分－微分型本构关系

应用　　　关于应力是五维偏应变空间变形历史的泛函的概念和Ｖａｌａｎｉｓ有关内＊时理论的描述，本文提出，对热粘塑性体，应力可设为应变、应变率和温度历史的泛函；并应用Ｍｉｌｌｅｒ和其它一些作者有关内变量演化方程的描述，由此建立了热粘塑性体的积分－微分到本构方程．这一积分－微分型本构关系大体和Ｍｉｌｌｅｒ微分型模型等价．对１０２０钢的单轴本构响应进行了数值模拟，和Ｔａｎａｋａ与Ｍｉｌｌｅｒ的分析及一些实验结果符合较好．     

NEPE推进剂低高应变率下改进的黏-超弹本构模型

为研究低高应变率条件下NEPE推进剂的力学特性,通过电子万能试验机和分离式霍普金森杆装置,对NEPE推进剂进行了准静态和冲击实验,得到了不同应变率下（1.667×10?4～4 500 s?1）的应力-应变曲线。实验结果表明NEPE推进剂具有明显的非线性弹性和应变率敏感性,随着应变率的增加,材料的强度、屈服应力和弹性模量显著增加,与低应变率相比,高应变率条件下材料的应变率敏感性更高。在高速冲击下材料内部瞬间产生大量热量无法及时散发出去,使得材料内部温度升高,导致材料出现软化效应,力学性能降低。本文建立了一个非线性黏超弹本构模型,其中采用Rivlin应变能函数来描述稳态超弹响应部分,采用积分型本构模型来描述材料的动态黏弹性响应部分,考虑到松弛时间具有应变率相关性,本文采用了一个率相关松弛函数来替代传统的Prony级数形式。使用极慢速压缩实验数据对本构模型中的超弹部分进行拟合获得超弹参数,然后用准静态和动态实验数据对本构模型进行拟合得出其他参数。不同应变率下的预测曲线与实验曲线具有较好的重合度,证明了该模型可以很好地描述低高应变率下NEPE推进剂的力学特性。     

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

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

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

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

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

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