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

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

An isotropic unilateral damage model coupled with frictional sliding for quasi-brittle materials

In this paper, we present an original extension of an isotropic damage model for quasi-brittle materials and assess its predictive capabilities. The proposed model accounts not only for unilateral behavior related to the opening and closure of microcracks but also for inelastic strains reflecting the frictional sliding along closed microcracks. More importantly, owing to its careful mathematical formulation with a particular attention paid to the continuous differentiability of the underlined thermodynamic potential, the model ensures the continuity of the inelastic stress–strain response. First applications show that it is able to predict the asymmetric behavior and hysteretic response of microcracked materials such as concrete and some rocks     

Constitutive relationship of brittle rock subjected to dynamic uniaxial tensile loads with microcrack interaction effects

A micromechanical model is proposed to describe both stable and unstable damage evolution in microcrack-weakened brittle rock material subjected to dynamic uniaxial tensile loads. The basic idea of the present model is to classify the constitution relationship of rock material subjected to dynamic uniaxial tensile loads into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress-strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to consistent with the experimental results.     

基于近场动力学理论的准脆性材料的本构力函数的构建

近场动力学理论(PD)是基于非局部思想的连续介质力学新理论,用于研究材料破坏问题。根据准脆性材料破坏的线性和非线性的力学行为,在初始微观弹脆性材料(PMB)的本构力函数中引入了键的损伤模型,将键的断裂过程分成了线性的弹性变形阶段和非线性的损伤变形阶段,以此构建了准脆性材料的本构力函数的基本形式。以典型的准脆性材料为例构建了其本构力函数,通过在压缩载荷下对含预制不同角度单裂纹缺陷的类岩材料的裂纹扩展进行PD数值模拟仿真,裂纹起裂位置和扩展方向与试样试验结果在一定程度上保持了一致,证明了该基于近场动力学理论的典型准脆性材料的本构力函数可用于该类材料的破坏分析。     

A micromechanics-based thermodynamic formulation of isotropic damage with unilateral and friction effects

This paper is devoted to micromechanical modeling of isotropic damage in brittle materials. The damaged materials will be considered as heterogeneous media composed of solid matrix weakened by isotropically distributed microcracks. The original contribution of the present work is to provide a proper micromechanical thermodynamic formulation for damage-friction modeling in brittle materials with the help of Eshelby’s solution to matrix-inclusion problems. The elastic and plastic strain energy involving unilateral effects will be fully determined. The condition of microcrack opening–closure transition will be determined in both strain-based and stress-based forms. The effect of spatial distribution of microcracks will also be taken into account. Further, the damage evolution law is formulated in a sound thermodynamic framework and inherently coupled with frictional sliding. As a first phase of validation, the proposed micromechanical model is finally applied to reproduce basic mechanical responses of ordinary concrete in compression tests.     

带微裂纹物体的有效断裂韧性

按照等效介质的思想，引进有效表面能密度的概念，建立了带微裂纹物体有效断裂韧性的公式．具体计算了微裂纹群分别平行和垂直于宏观裂纹两种情况的减韧比．表明微裂纹群在产生应力屏蔽（或反屏蔽）效应的同时，也降低了材料的有效断裂韧性，减小了对宏观裂纹的扩展阻力．     

Micromechanical analysis of damage in saturated quasi brittle materials

In this paper, we propose a micromechanical analysis of damage and related inelastic deformation in saturated porous quasi brittle materials. The materials are weakened by randomly distributed microcracks and saturated by interstitial fluid with drained and undrained conditions. The emphasis is put on the closed cracks under compression-dominated stresses. The material damage is related to the frictional sliding on crack surface and described by a local scalar variable. The effective properties of the materials are determined using a linear homogenization approach, based on the extension of Eshelby’s inclusion solution to penny shaped cracks. The inelastic behavior induced by microcracks is described in the framework of the irreversible thermodynamics. As an original contribution, the potential energy of the saturated materials weakened by closed frictional microcracks is determined and formulated as a sum of an elastic part and a plastic part, the latter entirely induced by frictional sliding of microcracks. The influence of fluid pressure is accounted for in the friction criterion through the concept of local effective stress at microcracks. We show that the Biot’s effective stress controls the evolution of total strain while the local Terzaghi’s effective stress controls the evolution of plastic strain. Further, the frictional sliding between crack lips generates volumetric dilatancy and reduction in fluid pressure. Applications of the proposed model to typical brittle rocks are presented with comparisons between numerical results and experimental data in both drained and undrained triaxial tests.     

A damage model with evolving nonlocal interactions

We present a damage model for softening materials with evolving nonlocal interactions. The thermodynamic implications and the material stability issue are addressed. The proposed nonlocal averaging scheme provides the obtained constitutive models with an evolving nonlocal interaction which is activated only when damage occurs. In the analysis of structures made of quasi-brittle materials, this feature helps not only to overcome some issues with the incorrect initiation of damage but also to better control the evolving size of the active fracture process zone. This is an essential feature that is usually not considered in depth in many existing nonlocal approaches to the continuum modelling of quasi-brittle fracture. Numerical examples are given to demonstrate features of the proposed modelling approach.     

Quasi-micromechanical constitutive theory for brittle damaged materials under tension

One of fundamental but difficult problems in damage mechanics is the formulation of the effective constitutive relation of microcrack-weakened brittle or quasi-brittle materials under complex loading, especially when microcrack interaction is taken into account. The combination of phenomenological and micromechanical damage mechanics is a promising approach to constructing an applicable damage model with a firm physical foundation. In this paper, a quasi-micromechanical model is presented for simulating the constitutive response of microcrack-weakened materials under complex loading. The microcracking damage is characterized in terms of the orientation domain of microcrack growth (DMG) as well as a scalar microcrack density parameter. The DMG describes the complex damage and its evolution associated with microcrack growth, while the scalar microcrack density factor defining the isotropic magnitude of damage yields an easy calculation of the effects of microcrack interaction on effective elastic moduli. Project supported by the National Natural Science Foundation of China (19891180).     

On the stored and dissipated energies in heterogeneous rate-independent materials

The present paper is a part of a work that aims at building a dissipative model of microcrack friction in quasi-brittle energetic materials. The latter is viewed as an assembly of elementary cells containing the most salient features of the microstructure heterogeneity. It is intended here to build an analytical model describing the mechanical and energetic response of such an elementary cell under confined tension. This is achieved by applying a previously published theory that allows for the determination of the amount of dissipated and stored energies in heterogeneous dissipative structures containing microcracks and other dissipative components.     

Tensile failure of two-dimensional quasi-brittle foams

Stress redistribution caused by damage onset and the subsequent local softening plays an important role in determining the ultimate tensile strength of a cellular structure. The formation of damage process zones with struts dissipating a finite amount of fracture energy will require the macroscopic stress to be increased in order to continue structural damage. The goal of this paper is to investigate the influence of the fracture energy of the solid on the tensile fracture strength and the strain to fracture in quasi-brittle two-dimensional foams using a microstructural model. We analyze the mesoscopic damage and failure mechanisms in uniaxial tension. Relative density, strut cross-sectional profile, solid’s fracture strain, and fracture energy are varied systematically. The effect of the specific fracture energy on the peak behavior has been shown to be captured by the ratio of the fracture energy to the stored elastic energy. We have also explored the net section strength variation in the presence of a central crack at two different fracture energies. Comparison is made between two structurally identical quasi-brittle and ductile strain hardening foams to identify the differences in the damage mechanisms.     

Modelling unilateral damage effect in strongly anisotropic materials by the introduction of the loading mode in damage mechanics

Damage weakens the mechanical characteristics of materials. But this weakening can disappear if the cracks close again: this is called the unilateral effect of damage. We propose a model of this phenomenon using damage mechanics in the case of a diffuse network of identical microcracks. The microcracks state is defined with two internal state variables. These two variables are control parameters of the geometry of the microcracks. They also define the loading mode in damage mechanics as in fracture mechanics. In order to limit the anisotropy induced by the microcracks, hence by the loading, we suppose that the geometry of damage spreads into preferential directions. Therefore, this model is essentially applied to materials with a strong anisotropy where the defects follow the constituents of the material. An application is given for composite laminates.     

Microcrack interaction brittle rock subjected to uniaxial tensile loads

A micromechanics-based model is proposed to describe unstable damage evolution in microcrack-weakened brittle rock material. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress–strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to be consistent with the experimental results.     

SBFEM与非局部宏微观损伤模型相结合的准脆性材料开裂模拟

混凝土是一种被广泛应用于土木和水利工程中的准脆性材料, 在各种内外部因素的作用下, 开裂是混凝土结构最为普遍的破坏形式, 准确模拟结构的开裂过程, 对于结构的安全评估至关重要. 将比例边界有限元与非局部宏微观损伤模型相结合提出一种准脆性材料开裂模拟新方法. 以比例边界有限元子域的比例中心作为物质点, 通过两比例中心(物质点对)之间的物质键的正伸长率来定义微细观损伤, 将某点影响域内物质键的微细观损伤加权平均得到该点的宏观拓扑损伤. 再引入能量退化函数, 将宏观拓扑损伤嵌入到比例边界有限元的基本框架中. 充分利用比例边界有限元网格允许存在悬挂节点的优势, 采用四叉树网格离散技术进行快速、高质量的多级网格划分与过渡. 通过一个I型开裂与一个混合型开裂的两个典型算例, 验证了该方法可捕获结构裂纹扩展路径与荷载变形曲线. 与现有的方法相比, 本文的损伤模型可得到更准确的局部开裂损伤带, 结果更为合理, 且具有更高的计算精度和计算效率. 当损伤过程区网格尺寸小于影响域半径的1/5时, 计算结果不存在网格敏感性问题.      

准脆性材料抗压强度能量平衡尺寸效应模型

针对现有尺寸效应模型难以体现准脆性材料完整的抗压强度尺寸效应变化规律及其内在机理, 本文通过分析准脆性材料单轴压缩破坏过程中能量输入、储存、整体和局部能量耗散, 建立体现整体和局部损伤的力学模型及描述上述能量演化过程的双线性名义和真实应力应变曲线, 在此基础上确定了名义应力最大时输入能量、储存弹性能、整体和局部能量耗散的表达式, 最后基于能量平衡原理建立抗压强度尺寸效应模型. 抗压强度能量平衡尺寸效应模型能完整体现名义抗压强度尺寸效应, 即随试样尺寸增大, 名义抗压强度在试样尺寸小于等于局部损伤区尺寸时为真实强度, 然后逐渐减小, 最终当试样尺寸趋于无穷大时趋于弹性极限强度; 抗压强度能量平衡尺寸效应模型也能同时体现高径比和试样直径对名义强度的影响, 其包含的参数具有明确的物理意义, 可以反映真实强度、弹性极限强度、名义损伤模量非线性、局部损伤区大小和方向对准脆性材料名义抗压强度尺寸效应的影响; 通过把抗压强度能量平衡尺寸效应模型和现有尺寸效应模型应用于预测各种材料尺寸效应试验和数值模拟数据, 结果表明: 抗压强度能量平衡尺寸效应模型能很好描述试验和数值模拟尺寸效应的非线性变化规律及内在机理, 和现有尺寸效应模型相比, 其总体平均误差最小, 且小于5%.      

Micromechanical concept for the analysis of damage evolution in thermo-viscoelastic and quasi-brittle materials

The aim of this paper is to develop a thermodynamically consistent micromechanical concept for the damage analysis of viscoelastic and quasi-brittle materials. As kinematical damage variables a set of scalar-, vector-, and tensor-valued functions is chosen to describe isotropic and anisotropic damage. Since the process of material degradation is governed by physical mechanisms on levels with different length scale, the macro- and mesolevel, where on the mesolevel microdefects evolve due to microforces, we formulate in this paper the dynamical balance laws for macro- and microforces and the first and second law of thermodynamics for macro- and mesolevel.Assuming a general form of the constitutive equations for thermo-viscoelastic and quasi-brittle materials, it is shown that according to the restrictions imposed by the Clausius–Duhem inequality macro- and microforces consist of two parts, a non-dissipative and a dissipative part, where on the mesolevel the latter can be regarded as driving forces on moving microdefects. It is shown that the non-dissipative forces can be derived from a free energy potential and the dissipative forces from a dissipation pseudo-potential, if its existence can be assured.The micromechanical damage theory presented in this paper can be considered as a framework which enables the formulation of various weakly nonlocal and gradient, respectively, damage models. This is outlined in detail for isotropic and anisotropic damage.     

A micromechanics-based strain gradient damage model for fracture prediction of brittle materials – Part II: Damage modeling and numerical simulations

In this paper, we established a strain-gradient damage model based on microcrack analysis for brittle materials. In order to construct a damage-evolution law including the strain-gradient effect, we proposed a resistance curve for microcrack growth before damage localization. By introducing this resistance curve into the strain-gradient constitutive law established in the first part of this work (Li, 2011), we obtained an energy potential that is capable to describe the evolution of damage during the loading. This damage model was furthermore implemented into a finite element code. By using this numerical tool, we carried out detailed numerical simulations on different specimens in order to assess the fracture process in brittle materials. The numerical results were compared with previous experimental results. From these studies, we can conclude that the strain gradient plays an important role in predicting fractures due to singular or non-singular stress concentrations and in assessing the size effect observed in experimental studies. Moreover, the self-regularization characteristic of the present damage model makes the numerical simulations insensitive to finite-element meshing. We believe that it can be utilized in fracture predictions for brittle or quasi-brittle materials in engineering applications.     

Statistical damage theory of 2D lattices: Energetics and physical foundations of damage parameter

The paper presents an in-depth analysis of two-dimensional disordered lattices of statistical damage mechanics for the study of quasi-brittle materials. The strain energy variation in correspondence to damage formation is thoroughly examined and all the different contributions to the net energy changes are identified and analyzed separately. We demonstrate that the introduction of a new defect in the microstructure produces a perturbation of the microscopic random fields according to a principle of maximum energy dissipation. A redistribution parameter η is introduced to measure the load redistribution capability of the microstructure. This parameter can be estimated from simulation data of detailed models. This energetic framework sets the stage for the investigation of the statistical foundations of the damage parameter as well as the damage localization. Logical statistical arguments are developed to derive two analytical models (a maximum field model and a mean field one) for the estimate of the damage parameter via a bottom-up approach that relates the applied load to the microstructural disorder. Simulation data provided input to the statistical models as well as the means of validation. Simulated tensile tests of honeycomb lattices with mechanical disorder demonstrate that long-range interactions amongst sets of microcracks with different orientations play a fundamental role already in damage nucleation as well as in the homogeneous–heterogeneous transition. A functional “hierarchy of sets” of grain boundaries, based on their orientation in relation to the applied load, seems to emerge from this study. Results put in evidence the ability of discrete models of capturing seamlessly the damage anisotropy. The ideas exposed inhere should be useful to develop a full rational model for disordered lattices and, later, to extend the approach to discrete models with solid elements. The findings suggest that statistical damage mechanics might aid in the quest of reliable and physically sound constitutive relations of damage, even in synergy with micromechanics.     

考虑岩石微元强度呈对数正态分布的冻融岩石损伤特性

冻融环境下准脆性岩石损伤力学特性对寒区工程建设有重要影响。由于准脆性岩石存在特征长度,其尺度效应不再符合Weibull尺度效应统计理论。本文以准脆性岩石为研究对象,考虑准脆性岩石微单元强度服从对数正态分布的条件下,依据D-P破坏准则,运用损伤力学理论建立了围压作用下冻融准脆性岩石的损伤本构模型,对冻融准脆性岩石损伤力学特性进行了探讨。研究表明,（1）采用对数正态分布建立的冻融准脆性岩石损伤力学本构模型与试验结果一致,该模型反映出冻融环境下准脆性岩石总损伤沿应变、围压和冻融循环路径相互作用,能够揭示出准脆性岩石损伤演化的本质特征。（2）冻融与围压共同作用下,准脆性岩石应力-应变变化关系可划分为线弹性阶段、非线性强化阶段、应力减弱阶段以及应变衰退阶段。（3）当准脆性岩石所受围压一定时,冻融次数是使准脆性岩石冻融受荷总损伤增大的关键因素。冻融循环次数对准脆性岩石的损伤劣化起促进作用;当准脆性岩石所受冻融循环次数不变时,冻融受荷总损伤随围压增大而减小,即围压对准脆性岩石力学性能劣化起抑制作用。相关研究成果对寒区岩体工程安全稳定性评价有一定的理论参考价值。