考虑微裂纹相互作用的的岩石微观力学弹塑性损伤模型研究

本文建立基于微裂纹扩展的岩石弹塑性损伤微观力学模型。用自洽方法考虑裂隙间相互影响,压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的微观损伤,基于应变能密度准则用Newton迭代法求复合型断裂的翼裂纹扩展长度,并采用微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映塑性内变量对硬化函数的影响;建立岩石的弹塑性损伤本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性损伤模型的程序。从围压和微裂隙长度等因素分析弹塑性损伤模型的岩石的损伤和宏观塑性特性。     

压缩载荷作用下岩石的细观损伤和塑性研究

建立岩石微裂纹扩展的细观力学模型,研究了岩石的细观损伤和塑性性质.压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的细观损伤,采用应变能密度准则求解复合型断裂的翼裂纹扩展长度,微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映了塑性内变量对硬化函数的影响;建立岩石模型的本构关系及其数值算法,并用回映隐式积分算法编制了模型的本构程序.分析弹塑性损伤模型的围压对岩石损伤的影响,并从围压和短微裂隙长度等因素分析模型的岩石的损伤和宏观塑性特性.     

弹丸侵彻混凝土的SPH算法

给出了弹丸侵彻混凝土的数值计算,其中弹体作为刚体处理并划分成Lagrangian标准有限元网格,而混凝土划分成光滑粒子并经历大应变、高应变率和高压作用。为了描述混凝土的非线性变形及断裂特性,在计算中引入了Holmquist Johnson Cook本构模型及损伤模型。计算结果与实验结果的对比表明,利用光滑粒子流体动力学方法对混凝土材料进行大应变、高应变率的变形计算是有效的,并可避免网格重分或网格消蚀。此外,将光滑粒子流体动力学方法和有限元方法结合可以保持计算过程中材料界面的清晰。     

弹塑性损伤本构关系的显式积分算法研究

首先引入弹塑性损伤本构关系,分别从材料软化与残余应变两个方面,描述伪脆性材料的非线性行为.针对结构动力分析中的强非线性问题,给出了弹塑性损伤本构关系的显式积分算法.算法中引入算子分解的思想,将弹塑性本构关系分成塑性与损伤两个模块.首先求解塑性模块,根据有效应空间塑性演化公式,采用前进欧拉算法,直接构造塑性演化的预测值,并且根据屈服函数的漂移构造了误差限公式,作为衡量显式算法精度的指标.将塑性模块求解的结果代入损伤模块,可以方便地求得损伤变量的演化,并最终得到更新后的应力.整个求解过程不需要迭代,可最大程度的算法稳定性.将论文建立的本构关系显式算法与结构分析显式算法结合,构造了结构显式分析方法,并模拟了两个经典算例,算例结果验证了论文方法的有效性.     

纯钽动态本构关系的实验研究

利用材料试验机和分离式霍布金森压杆实验装置测量了纯钽在较宽温度(-100～550℃)和应变率(10-5～103s-1)范围内的应力 应变曲线,该曲线表明钽的流动应力对温度和应变率都很敏感。利用所测得的应力 应变曲线拟合了钽在不同应变率和不同温度下的Johnson Cook与Zerilli Armstrong本构方程,Zeril li Armstrong本构方程与实验数据基本吻合,Johnson Cook本构方程与实验数据吻合得不太好。对Johnson Cook本构方程中的应变率强化项作了修正,修正后的Johnson Cook本构方程与实验数据吻合比修正前的好。     

混凝土黏塑性动力损伤本构关系

从静力弹塑性损伤本构关系的基本框架出发, 综合考虑塑性应变与损伤 演化的率敏感性, 建立了能够较为全面地描述混凝土在动力加载条件下非线性性能的混凝土 黏塑性动力损伤本构模型. 为了考虑塑性应变的率敏感性, 基于Perzyna理论推导了有效应 力空间黏塑性力学基本公式, 采用改进的Perzyna型动力演化方程, 将损伤静力演化方程推 广到动力加载情形. 基于并联弹簧模型, 从概率论的角度推导给出了一维损伤静力演化方程, 并基于能量等效应变的基本概念将其推广到多维损伤演化. 利用数值模拟, 计算得到了 混凝土在不同应变率下的应力应变全曲线, 同时得到了一维动力提高因子和二维动力强度包 络图, 数值结果与试验结果的对比表明了该模型的有效性.     

基于对数应力率大应变本构关系的参数标定研究

在所有率型弹塑性本构模型中,只有对数应力率对应的本构模型能够满足自适应准则.基于对数应力率,采用实心圆轴扭转实验,对大应变弹塑性本构模型中的参数标定问题进行了讨论.推导出了考虑Swift效应时端部自由实心圆轴扭转变形的变形率、对数旋率、Kirchhoff应力及Kirchhoff应力的对数应力率.对于等向强化大应变弹塑性本构关系,给出了由实心圆轴扭转实验标定的、基于Kirchhhoff应力对数应力率的本构关系中塑性刚度函数的表达式.分析了扭转圆轴的Swift效应对塑性刚度函数的影响.结果表明,实心圆轴扭转的轴向伸长变形和径向变形对基于对数应力率大应变本构关系中的塑性刚度函数都有影响.当不考虑Swift效应时,所得塑性刚度函数表达式与不考虑Swift效应时基于Jaumann应力率的塑性刚度函数表达式相同.     

聊聊动态塑性和黏塑性

固体力学研究者致力于具有应力-应变本构关系（以下简称为形变型本构关系）的变形体的力学响应研究，而流体力学研究者致力于具有应力-应变率本构关系（以下简称为流动型本构关系）的流动体的力学响应研究。当涉及结构和材料的动态塑性时，到底应该用“塑性变形”还是“塑性流动”来表示？本文从宏观塑性本构理论和微观位错动力学机理两个角度，分别讨论并指出塑性本构关系属于流动型黏塑性率相关本构关系，且同时适用于加载和卸载；因而不应该用应力-应变图来描述塑性加-卸载过程。弹塑性本构关系则是一种形变型和流动型本构关系的耦合。     

内爆炸加载条件下圆筒的膨胀、破裂规律研究

以经典热粘塑性本构关系 (John Cook本构关系和Powerlaw本构关系 )为基础 ,建立了更适当表征材料特征的新的、解耦的本构关系。采用SHPB(分离式霍普金森杆 )技术测定了常用弹箭材料 35CrMnSiA的不同应变率下的应力、应变关系 ,并采用拟合的方法确定了本构方程中的材料常数。与高速摄影技术测得的径向位移函数及圆筒材料的损伤演化方程相结合 ,建立了控制内爆炸加载圆筒膨胀 ,直至破裂的完备方程组 ,完成了内爆炸加载圆筒问题的数值模拟 ,计算结果与以往有关箭弹材料圆筒膨胀的实验结果符合较好。     

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

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

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.     

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

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

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

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