更全的杂志信息网

基于介观力学信息的颗粒材料损伤--愈合与塑性宏观表征1)

更新时间:2009-03-28

引言

颗粒材料广泛的存在于自然界与实际工程应用中,如土体,土石混合体,地质体以及混凝土结构中.在介观尺度,颗粒材料由大量不同形状、尺度的离散固体颗粒及空隙组成,为高度非均质、非连续介质.颗粒材料非线性、耗散等力学行为与其离散介观结构随力学响应的演变过程有着直接的联系.联接宏--介观力学行为模拟的多尺度计算分析方法,被用以研究颗粒材料的复杂力学行为.

说完,她跑回座位,趴在桌子上“嘤嘤”地哭起来。意料之中,张绍凡哭过后就去班主任那告状了,在校我挨了老师一顿批,回家又挨了老爸一顿揍。

在介观尺度上颗粒材料的每个离散颗粒具有独立的旋转自由度与平动自由度,在接触点处能同时传递力和力偶.因而在多尺度分析中,它在宏观尺度上模型化为Cosserat连续体,而不是Cauchy连续体[1-7].非均质材料的介--宏观计算均匀化模拟基于Hill定理,据此可以给定表征元(representation volume element,RVE)的静力学或运动学容许边界条件,使得Hill-Mendal宏--介观能量等价条件得到满足[8-11].

经典 Cosserat连续体的 Hill定理表明 [4,9],宏观Cosserat连续体仅能下传与施加均一性宏观应变场于表征元边界.作用于表征元均一宏观应变场的合理性基于尺度分离概念,考虑到表征元尺寸并非远小于宏观连续体有限元分析的特征尺寸,特别是在高宏观应变梯度(如应变局部化)区域尺度分离条件不再成立.为施加非均一宏观应变场于表征元边界,要求在宏观尺度采用梯度增强Cosserat连续体模型[10],并发展了相应的二阶计算均匀化方法[11]和在宏--介观尺度的混合元--离散元嵌套算法[12].

然而,基于离散颗粒系统--连续体模型的颗粒材料计算多尺度方法的表征元介观尺度分析并不直接提供和上传表征相应积分点处损伤、塑性等的内状态变量至宏观连续体,而这正是从事颗粒材料结构工程破坏分析的工程技术人员所关心和需要的数值模拟结果.本文主要工作是在颗粒材料计算均匀化多尺度分析框架下,建议一个基于介观力学的宏观连续体中材料损伤--塑性表征方法.

塑性的基本特征是连续体变形的不可恢复性以及在塑性变形过程中所做的塑性功的不可逆性.经典连续损伤力学中材料损伤定义为以弹性模量张量表示的弹性刚度降低[13-15].

“到了2000年左右,新课程改革了……改革的显著之处在于将‘双基目标’改为‘三维目标’……于是,我努力将自己的教学目标调整为‘三维目标’.可是,从此我发现,写教案的时候,我已经不会写教学目标了.因为我发现每节课都有特定的基本知识、基本技能,却很难区分出每节课的思想方法.”

材料塑性与损伤描述了不可逆的能量耗散过程.在连续介质力学唯象方法中,塑性与损伤理论研究基于热动力学框架与引入以唯象形式定义的内变量以定量描述材料内部微结构演变和所伴随的不可逆能量耗散[16-20].

损伤力学的近期研究表明,包括颗粒材料在内的许多工程材料在一定条件下具备愈合与部分恢复它们的刚度的潜能[21-25].需要拓展经典连续损伤力学理论,以计及具有非耗散本质的愈合过程[26].大部分模拟愈合过程的本构模型基于唯象理论[26-33].Ju等[31-32]首先对岩土材料提出了一个新的弹塑性--损伤--愈合唯象耦合模型.

其中是下面矩阵的转置

由于颗粒材料一般为各向异性材料,其损伤--愈合效应也一般为各向异性,甚至在一个方向损伤,而在另一个方向可能呈现材料愈合.区分梯度增强Cosserat连续体中局部材料点处的损伤或愈合行为的准则为d=0,即定义d>0表示材料损伤,d<0表示材料愈合.

本文从颗粒材料介观尺度的离散颗粒集合体力学响应分析出发,在二阶计算均匀化框架内对具有离散颗粒介观结构的表征元建立了等价连续体的增量非线性本构方程;把连续体本构关系中所定义的塑性--损伤--愈合行为以离散颗粒集合体中颗粒间相对耗散摩擦滑移、直接相邻颗粒间接触取向的改变、接触的丧失与再生等介观力学信息表示.

进一步,在颗粒材料热动力学框架[38]内定义了基于表征元介观力学信息的塑性应变以及各向异性损伤、愈合因子张量和综合了各向异性损伤、愈合效应的净损伤因子张量.进一步,以塑性、损伤耗散能与非耗散的愈合能等具有标量特征的内状态量表征并定量比较材料的塑性--损伤--愈合效应.应变局部化数值算例表明了基于介观力学信息的宏观梯度Cosserat连续体中塑性--损伤--愈合表征方法的有效性.

1 凝聚于表征元边界的离散颗粒集合体非线性增量本构方程

一对典型接触颗粒,在时间段[ti−1,ti],颗粒m的位移增量为颗粒n的位移增量为如图1所示.

按照加固施工图的要求,在需要加固的人行道托架梁体挡砟墙外侧定位螺栓位置(利用托架角钢上的螺栓孔在梁体上定位),用钢筋探测仪探测附近钢筋,确定孔位与梁体内钢筋不冲突。在梁体相应位置钻孔,钻孔深度140mm、钻孔直径28mm,注意过程不要碰到梁体内钢筋。钻孔完毕后,先用毛刷深入孔中清渣,然后将梁体表面的浮尘除掉,最后用气管清除孔中灰尘3遍,要求孔内干燥、干净无物。

Stephan等[29]的一项研究中对肿瘤组织和炎性病灶的达峰时间和流入时间进行测量,炎性病灶的达峰时间(tpeak 22 s)低于正常胰腺组织(tpeak 30 s),然而胰腺癌的达峰时间更长(tpeak 57 s)。流入时间为对比剂达到最大强度10 %的时间(tarr),在肿瘤病灶(tarr 26.2 s)中显著高于正常胰腺组织(tarr 14 s)和局灶性炎性病变(tarr 17 s)。在普通多普勒勒超声检查中发现的可疑病变,可以通过测量到达时间和达峰时间是否延迟,在为小胰腺癌的组织病理学的改变提供依据。

  

图1 典型的接触颗粒对mnFig.1 A typical contact between particle m and particle n

颗粒n为参考颗粒,按照右手螺旋法则定义与其接触颗粒m在接触点c处的局部坐标nt.颗粒n和颗粒m之间的相对接触位移增量为

 
 

考虑颗粒接触间相对切向滑动和滚动以及相对转动对接触力的影响,颗粒m作用于颗粒n的接触力增量在局部坐标下可表示为

 

其中为接触力矢量增量,分别为接触点处,法向、切向力增量与转动力矩增量;Kl为颗粒间接触点刚度阵;塑性力增量为其中

 

热力学框架下,与损伤--愈合因子张量增量Δd共轭的热力学广义力为[35]

 

其中Tg为局部到全局坐标转换矩阵,为参考全局坐标系的刚度阵,为全局坐标系下塑性力增量.

颗粒n的所有直接相邻接触颗粒作用于颗粒n的接触力求和表示为

 

其中Nm为颗粒n的直接相邻颗粒数目.整个表征元内离散颗粒集合体的内力与外力增量关系可以写成

 

其中 KII,KIB,KBI,KBB是按照表征元内部颗粒与边界颗粒划分的表征元刚度阵的分块子矩阵,ΔuI,ΔuB为内部颗粒和边界颗粒位移及转角增量,为内部颗粒与边界颗粒的塑性力分量增量,为作用于内部颗粒与边界颗粒形心处的外力增量.在收敛状态下可以忽略,通过表征元的边界接触点作用在边界颗粒上面,因此式(7)改写为

改革森林旅游行政管理体制,建立以资源资产产权为中心的管理和经营新体制,实现旅游管理的协调统一和监督全覆盖[2]。

 

其中式(9)可以写成

在对《老子》有了越来越多的研究之后,王庆节提供了一个新的视域。他从“物”(存在体)、“物论”切入,通过对老子哲学与海德格尔哲学的比较研究,提出了这样一个问题:“存在着一个老子的‘物论’吗?”③据王庆节的分析,“物论”是存在的,其论说自有合理的一面。这里暂且将“物论”搁置,只看“物”,《老子》的确十分强调“物”。王弼本《老子》中“物”出现了37次,笔者认为其“物”指存在体,但是其外延超过了许慎、顾野王所划定的范围。

 

其中可拆分成边界颗粒形心上的外力与外力矩两部分,ΔuB对应拆分成两部分,则拆分成两部分,式(10)可以写成

 

边界颗粒形心处力矩与边界接触点处的外力与力矩的关系可以写成

 

Formulation design key points for baby sunscreen cosmetics 6 67

 

其中对于第i个边界颗粒定义为

 

式中分别是第i个边界颗粒形心和此颗粒与边界接触点的位置向量.把式(12)代入式(11)得到

 

式中注意边界颗粒形心位移增量与边界颗粒和表征元边界接触点位移增量之间关系式(15)改写成

 

式中中的第i个子向量为通过表征元边界接触点作用于表征元第i个周边颗粒的外力与外力矩,且可表示为

 

其中,均为矩阵.

 

其中Nc为边界颗粒总数,分别是的子矩阵.

2 基于介观力学信息的宏观梯度Cosserat连续体非线性增量本构关系

颗粒材料梯度增强Cosserat连续体中局部样条点(即有限元网格积分点)处的应力等于与该点相关联的表征元的平均应力.它们可通过如下表征元的边界积分,边界积分可进一步离散到表征元中所有边界颗粒与表征元边界接触点处的力学与几何量,表示为

 

其中为Cauchy应力,为内扭矩,为应力矩的对称部分,为偶应力;VS为每个积分点处表征元体积与边界,xtm分别为表征元边界上任何一点的位置向量、面力与面力偶向量,e为三阶置换张量.

根据Hill定理、满足Hill-Mandel能量条件而应施加于表征元边界离散颗粒处的边界条件为[11]

 

其中为与表征元相关联的积分点处功共轭的宏观梯度Cosserat连续体应变度量.

利用式(18)和式(19)以及式(24)和式(25),式(20)~式(23)的增量形式可写成

 

首先看模型6,工作稳定度对社会距离的影响虽然是正向的,但是很小,并不具有统计显著性,因此没有把握将结论推论到研究总体,假设四没有得到统计分析上的支持。其原因是这种工作稳定度和主观社会经济地位有着比较大的相关性,因此其中很大一部分被主观社会经济地位解释掉了。

 

式(26)可以写成

首先对工艺过程进行HAZOP分析,确定是否采用SIS作为SIF的保护层,然后利用LOPA分析方法确定各SIF回路的SIL等级,最后采用紧急切断阀与储罐液位仪表相结合构建安全联锁保护系统,实现罐区的安全监测与自动控制。

 

其中为当前时刻包含16个弹性子矩阵[11]的弹性模量张量,为塑性应变增量,为对应的塑性应力增量,其计算形式如下所示

和蔼可亲的老师便于学生亲近,能促进师生的互动,情感的交流。老师知识的深度和广度直接影响着学生学习的兴趣。在语文教学实践中,我深深体会到“书到用时方恨少”这句话的含义。如果老师学识渊博,教学中能化抽象为具体,学生就会学的轻松,感受到学习的乐趣。教学语言清晰流畅,幽默有趣,让学生听得懂,听的有趣,学生学好就有了保证。教学方法灵活,让学生的学习处于新奇的状态,有利于激发学生的求知欲。因此,教师要不断修身、养德、广才,不断探索好的教学方法,以适应教学需求。

试验原料为四川会理地区红土镍矿,主要成分见表1,矿中铁以三氧化二铁、针铁矿等形式存在。将红土镍矿与硫酸铵进行混合焙烧,焙烧后铁以可溶性硫酸铁盐形式存在于熟料中。

 

式(27)为根据平均场理论与多尺度力学理论中Hill定理所导出的基于表征元介观结构与力学信息的联系宏观梯度Cosserat连续体局部材料点应力、应变、模量张量及塑性应变等物理量的非线性增量本构方程.它是下节中基于表征元介观信息的宏观连续体塑性及净损伤表征的出发点.

3 基于介观结构的损伤--愈合因子张量与损伤--愈合--塑性过程的热动力学耗散能

弹性损伤--愈合能密度增量为

在连续损伤力学框架内,损伤--愈合因子张量d可表示为

 

其中I为单位矩阵,是初始状态下的弹性模量张量,ti时刻的弹性模量张量,它们可由式(28)确定.经典时间段[ti−1,ti]内损伤--愈合因子张量增量为

 

分别为切向塑性接触力增量与塑性转动力矩增量,为切向相对滑动位移增量的耗散部分,为切向相对滚动位移增量的耗散部分为接触颗粒n,m间的相对转动位移增量的耗散部分,分别为接触颗粒间切向滑动摩擦力、切向滚动摩擦力、转动力矩刚度系数.把式(4)转换到全局坐标,可以写成

 

文献[35,38]推导了基于颗粒材料离散介观信息的等价经典Cosserat连续体Voronoi胞元的非线性弹性本构关系,在热动力学框架下定义了损伤--愈合因子张量,给出了损伤--愈合和塑性的耗散能量密度的表达形式.在此基础上,本文在颗粒材料二阶计算均匀化框架下,上节导出了梯度增强Cosserat连续体非线性弹性本构关系.

 

需要指出的是,在损伤--塑性耦合分析的宏观唯象理论框架内,损伤现象一般假定为各向同性[34].然而,对于包括颗粒材料在内的具有非均质和复杂介观结构的材料,材料的损伤与愈合在本质上为各向异性.各向异性来源于颗粒材料局部初始介观结构的各向异性以及由于颗粒材料力学响应导致的局部介观结构演变的各向异性[35-38].

根据式(27),宏观塑性应变增量可以表示为

 

塑性耗散能密度增量可表示为

 

4 数值算例

平板压缩算例,如图2(a)所示.忽略重力效应,因为满足对称性条件,数值模拟仅对1/4平板进行,如图2(b)所示,其中L=30m,被划分成20×20有限元网格.方板与刚性板之间假设为理性粘结,使得方板上边界节点水平和垂直方向自由度分别施加于有限元分析的边界条件如图1(b)所示.

  

图2 方板压缩问题Fig.2 Compression problem of a square panel

有限元网格中,每一个积分点处配置一个初始构形相同的表征元.3种具有相同规则介观结构,但具有不同尺寸的表征元,分别命名为RVE40,RVE60,RVE84,如图3所示,被用于考察表征元尺寸对宏观塑性--损伤--愈合效应表征的影响.每一个表征元尺寸是l×l,含nr个颗粒,每个颗粒半径为0.02m.颗粒材料性质数据参见文献[39].

作为一种先天畸形,盘状半月板损伤可使膝关节正常解剖结构发生明显改变,进而造成胫骨关节不协调、关节软骨磨损加快等现象的产生[6]。目前认为该疾病发生率与种族有较大相关性,且对于损伤后出现疼痛、绞锁、关节功能障碍等症状者需及时给予手术治疗,以帮助其缓解各相关症状,并促进膝关节功能的恢复[7,8]。

文献[12]详细讨论了二阶计算均匀化方法下所构造混合有限元的收敛性以及不需要宏观唯象本构与破坏模型,可以模拟颗粒材料失效与捕捉应变局部化现象.

  

图3 相同介观结构但不同尺寸的表征元Fig.3 Three samples of RVE with different RVE sizes but the same regular configuratio

本文以每一个积分点处累积的塑性及损伤耗散能密度及愈合能密度作为内状态变量表征宏观塑性、损伤与愈合,并定量比较它们对颗粒材料破坏的效应.

图4~图7所示分别为采用表征元 RVE40,RVE60,RVE84的方板应变局部化问题模拟加载结束时位移载荷曲线与耗散能密度分布云图.颗粒材料等价梯度增强Cosserat连续体宏观损伤--愈合现象的介观机理为表征元内颗粒接触丧失与再生和直接相邻颗粒间接触取向的改变,以及相伴随的表征元体积增大或减小,这些因素的综合将导致连续体局部材料点处弹性刚度的降低(损伤)或增大(愈合),塑性耗散的介观机理为表征元内颗粒在接触点处的耗散性相对滑移、滚动.颗粒材料等价梯度增强Cosserat连续体软化结束时即加载结束时,方板损伤--愈合以及塑性耦合过程的热力学耗散主要集中在应变局部化带.对比图5(c)、图5(d)至图7(c)、图7(d),可以看到,表征元对应的等价梯度增强Cosserat连续体材料点处,表征颗粒间耗散性滑移、滚动的塑性演化首先发生,表征颗粒间接触丧失的等价梯度增强Cosserat连续体局部点的体积膨胀与净损伤演化随后发生.

  

图4 位移载荷曲线Fig.4 Curves of load-displacement

  

图5 RVE40加载结束耗散能与愈合能密度分布(N/m)Fig.5 Contours of density of dissipative and healing energy distribution in the square panel at the end of the load history obtained by RVE40(N/m)

  

图5 RVE40加载结束耗散能与愈合能密度分布(N/m)(续)Fig.5 Contours of density of dissipative and healing energy distribution in the square panel at the end of the load history obtained by RVE40(N/m)(continued)

  

图6 RVE60加载结束耗散能与愈合能密度分布(N/m)Fig.6 Contours of density of dissipative and healing energy distribution in the square panel at the end of the load history obtained by RVE60(N/m)

  

图6 RVE60加载结束耗散能与愈合能密度分布(N/m)(续)Fig.6 Contours of density of dissipative and healing energy distribution in the square panel at the end of the load history obtained by RVE60(N/m)(continued)

  

图7 RVE84加载结束耗散能与愈合能密度分布(N/m)Fig.7 Contours of density of dissipative and healing energy distribution in the square panel at the end of the load history obtained by RVE84(N/m)

  

图7 RVE84加载结束耗散能与愈合能密度分布(N/m)(续)Fig.7 Contours of density of dissipative and healing energy distribution in the square panel at the end of the load history obtained by RVE84(N/m)(continued)

图8所示为采用RVE40,RVE60,RVE84等3种具有相同介观结构而不同尺寸表征元样本的方板全域随加载过程的塑性总耗散能、净损伤总耗散能与塑性--净损伤总耗散能的演变曲线.可以看到,当采用RVE40样本时,导致结构失效的总净损伤耗散能远小于总塑性耗散能;当采用RVE60样本时,总净损伤耗散能尚小于总塑性耗散能,但已较接近;而在采用RVE84样本的情况下,总净损伤耗散能与总塑性耗散能基本持平.从这一方面看,表征元的大小对表征结果是有影响的.

另一方面,图9所示为分别采用RVE40,RVE60和RVE84三种样本情况下方板的塑性--净损伤总耗散能随加载过程的演变曲线比较.从图9可以看出,3条曲线比较接近,特别是在软化后直至加载结束这一后软化阶段的总耗散能曲线重合很好.从这一方面看,表征结果并不病态地依赖于表征元尺寸.图4显示,分别采用RVE40,RVE60和RVE84三种样本情况下方板的外力耗散功十分接近,表明本文多尺度表征工作所基于的协同(concurrent)二阶计算均匀化过程的颗粒材料结构多尺度模拟[11−12]结果的可靠性,也即多尺度模拟结果不应病态地依赖于表征元的(窗口)尺寸.图9所给出的多尺度表征结果与图4显示的多尺度模拟结果比较,验证了本文工作在数值误差意义上近似满足了结构的外力耗散功模拟结果应等于结构的内耗散能表征结果这一基本要求,表明了本文提出的基于一致二阶计算均匀化过程的颗粒材料塑性--损伤--愈合多尺度表征方法的有效性.

  

图8 由所有积分点处耗散能密度累计的方板塑性、净损伤和总耗散能随加载过程的演化Fig.8 Evolutions of accumulated net damage,plasticity and total dissipative energies added up from those associated to all integration points of the square panel with increasing prescribed vertical displacement

  

图9 由所有积分点处耗散能密度累计的方板总耗散能随加载过程的演化Fig.9 Evolution of accumulated total dissipative energies added up from those associated to all integration points of the square panel with increasing prescribed vertical displacement

5 结论

颗粒材料在宏观和介观尺度分别模型化为梯度增强Cosserat连续体和离散颗粒集合体.本文基于颗粒材料力学行为的协同二阶计算均匀化模拟,对具有离散颗粒介观结构的表征元建立了等价梯度增强Cosserat连续体的增量非线性本构关系.

通过2组测站可以看出,顶板锚索的受力变化规律基本相同,呈波动状态;但两帮的锚杆受力状况差别较大,左帮锚杆卸载严重,右帮锚杆平均受力大于左帮,没有明显卸载现象。这是因为巷道掘进方向与原岩主应力方向之间的夹角和煤层倾角,使得两帮煤体的受力状态不同,从而导致致节理剪胀破坏的形态不同,锚杆受力规律也不同。左帮煤体卸荷严重,平行于巷帮平面剪胀变形大,导致锚杆载荷下降较大。因此,对于节理化巷道,还应进一步研究节理方向、开挖方向、应力分布特征等对巷道变形破坏和支护体的影响。

在颗粒材料热动力学理论框架下,定义了基于介观结构与介观力学响应的塑性应变以及各向异性损伤、愈合因子张量和综合了各向异性损伤、愈合效应的净损伤因子张量,提出了以净损伤和塑性耗散能(标量)表示的多尺度表征方法.应变局部化的数值例题结果显示了所提出表征方法的有效性.

参考文献

1 Alonso-marroquin F.Static equation of the Cosserat continuum derived from intra-granular stress.Granular Matter,2011,13(3)∶189-196

2 Chang CS,Kuhn MR.On virtual work and stress in granular media.International Journal of Solids and Structures,2005,42(13)∶3773-3793

3 D’Addetta GA,Ramm E.A particle center based homogenization strategy for granular assemblies.Engineering Computation,2004,21(2-3-4)∶360-383

4 Li XK Liu QP,Zhang JB.A micro-macro homogenization approach for discrete particle assembly Cosserate continuum modeling of granular materials.International Journal of Solids and Structures,2010,47(2)∶291-303

5 Ehlers W,Ramm E,Diebels S,et al.From particle ensembles to Cosserat continua∶Homogenization of contact forces towards stresses and couple stresses.International Journal of Solids and Structures,2003,40(24)∶6681-6702

6 Kruyt NP.Statics and kinematics of discrete Cosserat-type granular materials.International Journal of Solids and Structures,2003,40(3)∶511-534

7 PasternakE,MuhlhausHB.Generalisedhomogenisationprocedures for granular materials.Journal of Engineering Mathematics,2005,52(1-3)∶199-229

8 Qu J,Cherkaoui M.Fundamentals of Micromechanics of Solids.New Jersey∶Wiley,2006

9 Li XK,Liu QP.A version of Hill’s lemma for Cosserat continuum.Acta Mechanica Sinica,2009,25∶499-506

10 Li XK,Zhang JB,Zhang X.Micro-macro homogenization of gradient-enhanced Cosserat media.European Journal of Mechanics A,2011,30(3)∶362-372

11 Li XK,Zhang X,Zhang JB.A generalized Hill’s lemma and micromechanically based macroscopic constitutive model for heterogeneous granular materials.Computer Methods in Applied Mechanics and Engineering,2010,199(49-52)∶3137-3152

12 Li XK,Liang YB,Quan QL,et al.A mixed finit element procedureofgradientCosseratcontinuumforsecond-ordercomputational homogenization of granular materials.Computational Mechanics,2014,54(5)∶1331-1356

13 Kachanov L.Time rupture under creep conditions. Izvestiya Akademii Nauk SSSR,Otdelenie Tekhnicheskikh Nauk,1958 8(8)∶26-31

14 Lemaitre J,Chaboche JL.Mechanics of Solid Materials.Cambridge∶Cambridge University Press,1990

15 唐炳涛.基于连续体损伤理论的硼钢高温成形极限确定法.力学学报,2016,48(1)∶146-153(Tang Bingtao.Prediction of forming limit of boron steel at elevated temperatrue based on CDM theory.ChineseJournalofTheoreticalandAppliedMechanics,2016,48(1)∶146-153(in Chinese))

16 Chaboche JL.Continuum damage mechanics∶Part I∶General concept,Part II damage growth,crack initiation and crack growth.Journal of Applied Mechanics,1988,55(3)∶59-71

17 Chow CL,Wang J.An anisotropic theory of elasticity for continuum damage mechanics.International Journal of Fracture,1987,33(1)∶3-16

18 Ju JW.On energy-based coupled elastoplastic damage theories∶constitutive modeling and computational aspects.International Journal of Solids and Structures,1989,25(7)∶803-833

19 Simo JC,Ju JW.Strain-based and stress-based continuum damage models,part I∶Formulation.International Journal of Solids and Structures,1987,23∶821-840

20 Ladeveze P,LeDantec E.Damage modelling of the elementary ply for laminated composites.Composites Science and Technology,1992,43(3)∶257-267

21 DarabiMK,Al-RubRKA,LittleDN.Acontinuumdamagemechanics framework for modeling micro-damage healing.International Journal of Solids and Structures,2012,49(3-4)∶492-513

22 Herbst O,Luding S.Modeling particulate self-healing materials and application to uni-axial compression.International Journal of Fracture,2008,154(1-2)∶87-103

23 Mergheim J.,Steinmann P.Phenomenological modelling of selfhealing polymers based on integrated healing agents.Computational Mechanics,2013,52(3)∶681-692

24 Aliko-Ben´ıteza A,Doblar´e M,Sanz-Herrera JA.Chemical-diffusive modeling of the self-healing behavior in concrete. International Journal of Solids and Structures,2015,69-70∶392-402

25 陈洁,刘剑兴,姜德义等.围压作用下岩盐应变与损伤恢复试验研究.岩土力学,2016,37(1)∶105-112(Chen Jie,Liu Jianxing,Jiang Deyi,et al.An experimental study of strain and damage recovery of salt rock under confinin pressures.Rock and Soil Mechanics,2016,37(1)∶105-112(in Chinese))

26 BarberoEJ,GrecoF,LonettiP.Continuumdamage-healingmechanics with application to self-healing composites.International Journal of Damage Mechanics,2005,14(1)∶51-81

27 Voyiadjis GZ,Shojaei A,Li G.A thermodynamic consistent damage and healing model for self healing materials.International Journal of Plasticity,2011,27(7)∶1025-1044

28 Voyiadjis GZ,Shojaei A,Li GQ,et al.Continuum damage-healing mechanicswithintroductiontonewhealingvariables.International Journal of Damage Mechanics,2012,21(3)∶391-419

29 VoyiadjisGZand KattanPI.Healingandsuper healingincontinuum damage mechanics.International Journal of Damage Mechanics,2014,23(2)∶245-260

30 张浪平,尹祥础,梁乃刚.地震条件下损伤--愈合模型的初步研究.岩土力学与工程学报,2008,27(2)∶3956-3962(Zhang Langping,Yin Xiangchu,Liang Naigang.Preliminary study on damagehealing model under earthquake.Chinese Journal of Rock Mechanics and Engineering,2008,27(2)∶3956-3962(in Chinese))

31 Ju JW,Yuan KY,Kuo AW.Novel strain energy coupled elastoplastic damage and healing models for geomaterials-Part I∶Formulations.International Journal of Damage Mechanics,2012 21(4)∶525-549

32 Ju JW,Yuan KY.New strain energy based coupled elastoplastic twoparameter damage and healing models for earth moving processes.InternationalJournalofDamageMechanics,2012,21(7)∶989-1019

33 屈家旺,刘泉生,何军等.泥岩弹塑性损伤--愈合模型研究.岩土力学与工程学报,2014,33(1)∶3192-3197(Qu Jiangwang,Liu Quansheng,He Jun,et al.Study of elastoplastic damage-healing model for argillite,Chinese Journal of Rock Mechanics and Engineering,2014,33(1)∶3192-3197(in Chinese))

34 Li XK,Duxbury PD,Lyons L.Coupled creep-elastoplastic-damage analysis for isotropic and anisotropic nonlinear materials.International Journal of Solids and Structure,1994,31(9),1181-1206

35 Li XK,Du YY,Duan QL.Micromechanically informed constitutive model and anisotropic damage characterization of Cosserat continuumfor granular materials.International Journal of Damage Mechanics,2013,22(5),643-682

36 付云伟,倪新华,刘协权等.颗粒缺陷相互作用下复合材料的细观损伤模型.力学学报,2016,48(6)∶1334-1342(Fu Yunwei,Ni Xinhua,Liu Xiequan,et al.Micro-damage model of composite materials with particle and defect interaction.Chinese Journal of Theoretical and Applied Mechanics,2016,48(6)∶1334-1342(in Chinese))

37 郭洪宝,王波,贾普荣等.平纹编织陶瓷基复合材料面内剪切细观损伤行为研究.力学学报,2016,48(2)∶361-368(Guo Hongbao,Wang Bo,Jia Purong,et al.Mesoscopic damage behaviors of plain woven ceramic composite under in-plane shear loading.Chinese Journal of Theoretical and Applied Mechanics,2016,48(2)∶361-368(in Chinese))

38 Li XK,Du YY,Duan QL,et al.Thermodynamic framework for damage-healing-plasticity of granular materials and net damage variable,International Journal of Damage Mechanics,2016,25(2)∶153-177

39 Li XK,Wang ZH,Liang,YB et al.Mixed FEM-crushable DEM nested scheme in second-order computational homogenization for granular materials.International Journal of Geomechanics,2016,16(5)∶C4016004

 
王增会,李锡夔
《力学学报》 2018年第02期
《力学学报》2018年第02期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号