水泥基材料在宽围压范围的本构模型

本文首先进行了水泥石的三轴试验和静水试验,根据试验结果分析了宽围压范围下不同力学机理对水泥石宏观力学性能的影响。提出了基于热力学定律的剪切塑性、孔隙塑性与破坏准则相互作用的弹塑性损伤本构模型。该模型有效的考虑了围压对材料主要力学特征的影响,提出了与损伤准则相关联的塑性硬化函数。数值解答与试验结果比较表明,该本构模型可以很好地描述水泥石在低围压状态下的脆性、中围压下的脆-塑性转换以及高围压下的孔隙塑性为主的力学特征。     

粉砂岩峰后破坏区应力脆性跌落的试验和本构方程研究

对粉砂岩试样进行常规三轴压缩试验,得到了其在不同围压下的应力-应变全程曲线.结果表明,岩样应力-应变曲线峰前(破坏前区)部分的力学性质较为稳定,其力学行为可用经典强度理论进行描述;应力-应变曲线峰后(破坏后区)部分,岩样处于非稳定状态,其力学行为难以用经典强度理论来描述,而脆塑性模型恰好可以描述试验岩样的应力非连续变化特征.岩石的应力脆性跌落是有条件发生的,应力脆性跌落系数是围压的函数,并给出应力脆性跌落系数与围压的关系.对完整岩样而言,使用弹性-线性软化-残余塑性三线性非理想脆塑性模型可以描述岩样的本构特征,而对含有节理岩样而言,采用双线性弹性-线性软化-残余塑性四线型非理想脆塑性模型与试验更为相符.     

考虑围压效应的冻结砂土动态本构模型研究

为描述主动围压作用下冻结砂土的动态力学特性，通过在朱-王-唐模型的非线性体上串联塑性体，建立了能够考虑围压效应的冻结砂土动态损伤本构模型；分析了损伤参数对应力-应变曲线特征、屈服点、峰值应力和峰值应变的影响规律，基于冻结砂土动力学试验数据确定了模型参数；通过将模型和试验数据进行对比，并对不同试验条件下模型的预测误差进行分析，验证了模型的适用性和准确性。结果表明，损伤参数对应力-应变曲线弹性阶段和屈服点无明显影响，而对塑性阶段和破坏阶段的影响较为显著，本构模型预测的应力-应变曲线与试验结果具有较好的一致性。模型能够预测围压引起冻结砂土塑性阶段占比大和屈服点明显的特征，且能够描述围压对冻结砂土动态强度的增强效应；不同负温和主动围压条件下，模型对峰值应力和屈服强度的预测效果优于峰值应变和屈服应变。     

混凝土单轴受拉的非局部本构模型

混凝土受拉本构行为存在很强的局部软化现象,使得单轴受拉试验无法给出应力-应变关系,而只能给出应力-位移关系。本文根据内变量理论和等效应变假设建立了基于真实应变的混凝土单轴受力本构方程,并根据Weibull分布可以描述混凝土等脆性材料断裂过程的试验现象,建立了关于弹性应变的损伤演化规律。然后,通过假设平均应变与真实弹性应变的函数关系,在应力-平均应变的本构关系中采用平均弹性应变以描述其非局部行为,而在材料的损伤演化规律中采用真实弹性应变以描述其局部行为,由此建立了单轴受拉荷载条件下的非局部本构模型。最后,对一个单调受拉试验和一个反复受拉试验的仿真结果表明所提出的非局部本构模型可以准确地模拟试验结果。     

软化Mohr-Coulomb材料强度参数的测定方法

为了得到试件的粘聚力和内摩擦角随轴向塑性压应变变化的曲线提出本方法。试件的弹塑性本构关系遵循相关联的Mohr-Coulomb强度准则;对常规三轴试验,试件受力进入塑性状态后,处在棱椎状屈服面的棱上,加载过程遵循Koiter流动法则。按经典塑性力学理论,推导得到轴向塑性压应变与轴向应力与轴向应变的关系;在常规三轴试验机上获得不同围压下试件的全程应力-应变曲线,进而可得到各自围压下轴向塑性压应变随加载过程的变化曲线;把来自不同围压下对应同一轴向塑性压应变的应力分别代入屈服面方程,即可求得对应的粘聚力和内摩擦角。结果表明,Mohr-Coulomb材料的两个强度参数的变化由轴向塑性压应变确定。轴向塑性压应变可以作为塑性变形的状态参数,它和试件的受力过程可以唯一确定试件的变形过程。     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

基于0.2定围压比状态下的约束混凝土性能研究

通过对混凝土施加0.2定围压比,使单轴受压混凝土裂缝发展受到约束,研究了混凝土在接近无裂纹发展状态下的受压力学性能。本文以在MTS815.02岩石电液伺服系统上对不同强度等级的混凝土试件在定围压比为0.2时的单调、循环加卸载受压常规三轴试验为基础,研究分析了0.2定围压比对混凝土受压应力-应变规律的影响和围压对混凝土裂缝发展的限制效果,并根据0.2定围压比下混凝土材料的特殊性能,提出了无裂纹发展状态下混凝土上升段表达式及参数。分析结果表明,0.2定围压比下,各强度混凝土本构关系均表现出了类似金属材料的本构关系,且低强混凝土试件应力-应变的提高幅度优于中强混凝土试件。综上所述,0.2定围压比改善了混凝土在受压时无征兆失效(脆性破坏)的力学性能。     

混凝土的弹塑性损伤双面本构模型

针对混凝土材料拉压应变空间下损伤机制的不同,结合连续损伤力学和塑性理论建立了一个全新的本构模型。该模型中损伤和塑性变形的演变由应变空间的同一个非弹性曲面来控制,但对拉压应变空间中非弹性曲面的演变分别采用了随动强化法则和各向同性演化规律。计算结果表明,该模型能较好地描述混凝土材料在单轴及多轴单调加载和低周反复荷载下的典型非线性特征。     

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

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

基于元件组合理论的砂岩动态损伤本构模型

采用分离式霍普金森压杆（SHPB）系统,对砂岩进行不同速度下的冲击试验,得到砂岩的应变率效应特征以及典型的动态本构曲线。该曲线分为近似线弹性阶段、塑性阶段、塑性增强阶段和正向卸载阶段。通过组合模型的方法,构建了砂岩含损伤的动态本构模型,借助LS-DYNA软件中的用户材料子程序UMAT接口实现对本构模型的二次开发,并对砂岩在冲击速度为7.5、9.5、11.5和13.5 m/s 4种情况下的SHPB动态冲击压缩试验进行模拟。结果表明:所构建的模型可以很好地描述砂岩的应变率效应和应力-应变曲线弹性段,并且动态峰值强度、最大应变均与试验结果一致,应变率、峰值强度、最大应变与试验结果的相对误差不超过10%。所构建的砂岩动态本构模型能够准确地描述砂岩在冲击作用下的动态力学特性。     

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

正~~     

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

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

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

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

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。