3-D numerical simulation of fracture processes in heterogeneous brittle materials

By using the lattice model combined with finite element methods and statistical techniques, a numerical approach is developed to establish mechanical models of three-dimensional heterogeneous brittle materials. A special numerical code is introduced, in which a lattice model and statistical approaches are used to simulate the initial heterogeneity of material properties. The size of displacement-load step is adaptively determined so that only few elements would fail in each load step. When the tensile principal strain in an element exceeds the ultimate strain of this element, the element is considered broken and its Young's modulus is set to be very low. Some important behaviors of heterogeneous brittle materials are indicated using this code. Load-displacement curves and figures of three-dimensional fracture patterns are also numerically obtained, which are similar to those observed in laboratory tests.     

三维非均匀脆性材料破坏过程的数值模拟

采用有限元方法模拟了三维均匀固体材料在宏观等效力学性质和破坏过程。首先采用格形（lattice）方法把试件离散成三维均匀网格,在每个单元格中将材料按照均匀处理,根据给定的统计规则来确定不同单元格中的材料常数以反映材料的非均匀性。然后对非均匀脆性材料选用简单的本构关系与断裂准则,采用自适应选取载荷步长对试件进行加载,通过非平衡迭代技术对刚度矩阵进行不断修正,实现了非均匀脆性材料的弹性行为及破坏过程的数值模拟。在此基础上,通过数值计算研究了材料的非均匀分布对宏观等效力学性质和破坏过程的影响,给出了破坏全过程的非线性载荷－位移曲线以及不同载荷阶段的三维损伤破坏的演化图。     

The application of the Eshelby equivalent inclusion method for unifying modulus and transformation toughening

When a crack is lodged in an inclusion, both difference between the modulus of the inclusion and matrix material and stress-free transformation strain of the inclusion will cause the near-tip stress intensity factor to be greater (amplification effect) or less (shielding or toughening effect) than that prevailing in a homogeneous material. In this paper, the inclusion may represent a second phase particle in composites and a transformation or microcracked process zone in brittle materials, which may undergo a stress-free transformation strain induced by phase transformation, microcracking, thermal expansion mismatch and so forth. A close form of solution is derived for predicting the toughening (or amplification) effect. The derivation is based on Eshelby equivalent inclusion approach that provides rigorous theoretical basis to unify the modulus and transformation contributions to the near-tip field. As validated by numerical examples, the developed formula has excellent accuracy for different application cases.     

热环境中旋转运动功能梯度圆板的强非线性固有振动

研究热环境中旋转运动功能梯度圆板的非线性固有振动问题．针对金属－陶瓷功能梯度圆板,考虑几何非线性、材料物理属性参数随温度变化以及材料组分沿厚度方向按幂律分布的情况,应用哈密顿原理推得热环境中旋转运动功能梯度圆板的非线性振动微分方程．考虑周边夹支边界条件,利用伽辽金法得到了横向非线性固有振动方程,并确定了静载荷引起的静挠度．用改进的多尺度法求解强非线性方程,得出非线性固有频率表达式．通过算例,分析了旋转运动功能梯度圆板固有频率随转速、温度等参量的变化情况．结果表明,非线性固有频率随金属含量的增加而降低;随转速和圆板厚度的增大而升高;随功能梯度圆板表面温度的升高而降低．     

2D编织陶瓷基复合材料应力-应变行为的试验研究和模拟

本文对2D编织陶瓷基复合材料拉伸应力-应变行为进行了试验研究和理论模拟。将2D编织结构简化为：正交铺层结构和纤维束波动结构。基于基体随机开裂、纤维随机断裂的统计分布理论,得到正交铺层结构的应力-应变关系;基于体积平均方法,将纤维束波动部分分割为若干子单元;由于纤维束的波动使各子单元材料方向与加载方向不一致,因此考虑了各子单元的线性行为和非线性行为对材料响应的影响,同时引入强度分析模型,得到纤维束波动部分的应力-应变关系。结合正交铺层部分和纤维束波动部分的应力-应变关系,得到2D编织结构的应力-应变行为,理论与试验吻合较好。     

岩石紧凑拉伸断裂过程及其尺寸效应的三维数值模拟

采用最近开发的三维岩石破裂过程分析软件RFPA3D模拟单边裂纹紧凑拉伸断裂过程。试验中五个不同尺寸的岩样具有相同的力学性质参数分布,模拟结果得到了裂纹扩展中的应力场、位移场和声发射的空间分布以及单边裂纹扩展贯通的过程。单边裂纹拉伸断裂的路径是一个复杂的空间三维曲面,三维裂纹比二维裂纹更为复杂。分析了岩石试样的峰值强度和试样尺寸之间的关系。随着岩样尺寸的增加,峰值强度逐渐减小,并且延性破坏特征更加明显,模拟结果满足岩石的尺寸效应规律。最后模拟了三组不同均匀性的试样拉伸破坏过程,结果表明细观上的非均匀性对岩石尺寸效应有很大影响,随着非均匀性的增加,岩石宏观强度随之提高,即使在均匀材料中一样存在尺寸效应。     

岩石单轴压缩作用下变形局部化的梯度塑性解

采用梯度塑性理论研究单轴压缩作用下岩石变形局部化,得到了单轴压缩作用下岩石变形局部化带宽度的一维、二维解析解,为实验测定内部材料长度参数提供了理论依据.     

含水率和围压对安家岭泥岩峰后力学特性影响的试验研究

为分析含水率和围压对泥岩峰后力学特性的影响,从山西安家岭矿取泥岩,制成不同含水率的5组试样,在YAW2000电液伺服试验机上开展了三轴试验,获得了不同含水率泥岩试样三轴全程应力应变曲线,使用激光共聚焦显微镜观测了泥岩增水过程中微观结构变化。利用试验结果,分析了围压和含水率对泥岩峰值强度、残余强度、弹性模量、破坏应变和脆性模量的影响规律和泥岩增水过程中微观结构的变化规律。引入脆性模量系数和强度退化指数来描述围压对泥岩峰后强度退化过程和残余强度的影响,与FLAC中的SS模型结合,建立了考虑围压影响的泥岩应变软化力学模型,模拟了围压对泥岩应变软化行为的影响。研究结果表明:(1)随着围压增加,泥岩的峰值强度、偏应力峰值、破坏应变和残余强度都增长,峰后强度降低速率趋缓,强度退化指数和脆性模量系数可以较好地描述围压对泥岩残余强度和峰后强度退化过程的影响。(2)泥岩增水过程中,岩样内微裂隙及尺寸增长,泥岩的力学性质劣化。随着含水率增加,泥岩的弹性模量、峰值强度和残余强度降低,破坏应变增长。含水率与泥岩的弹性模量、峰值强度和破坏应变之间近似服从线性关系。(3)本文基于脆性模量系数的岩石应变软化模型能较好地描述三轴压缩泥岩的全程变形行为。     

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

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

低应变率下硅铝质泡沫陶瓷复合材料的力学性能研究

以泡沫陶瓷复合材料在防护工程中的应用为背景,利用MTS(Material Test System,材料试验机)对该型材料进行了准静态压缩实验。得到了应变率在10-5~10-3s-1范围内的应力应变曲线,并对实验结果进行了理论分析和数值模拟。研究表明,泡沫陶瓷复合材料的力学性能在准静态一维应力压缩条件下显示出明显的应变率效应,同时其应力应变曲线可用一种经验的脆性材料本构模型进行较好地拟合。而在一维应变压缩条件下,材料的应力应变曲线则显示出明显的三段式特征:弹性段、平台段和密实段,同时材料的吸能幅值随着应变率的增大而增加。     

单轴压缩下页岩类脆性材料细观损伤演化的离散元分析

了解广泛存在的类似页岩的脆性材料各向异性对工程安全具有重要意义。本研究将页岩视为粘结颗粒材料,基于离散单元方法研究了横观各向同性脆性页岩的损伤演化。再现了不同层理角的页岩试样的破坏模式,并对比了实验和数值模拟的抗压强度和弹性模量。引入微裂纹的概念,通过定义裂纹密度函数,系统地研究了单轴压缩条件下,页岩层理角对细观结构的影响。此外基于平均配位数建立了配位数变化与细观损伤的联系,并根据配位数的变化与裂纹数量将加载过程分为三个阶段,分析了不同阶段配位数与裂纹数量的对应变化关系。研究表明,页岩的裂纹密度随着层理角的增加而增加,而试样的平均配位数在加载过程中先上升后剧烈下降,颗粒集合体在单轴压缩条件下的应力应变及裂纹数量曲线与平均配位数曲线有良好的一致性。该研究揭示了横观各向同性脆性岩石的破坏过程和内在机理,将为页岩类脆性材料的工程应用提供理论指导。     

Strength and fracture toughness of heterogeneous blocks with joint lognormal modulus and failure strain

We obtain analytical approximations to the probability distribution of the fracture strengths of notched one-dimensional rods and two-dimensional plates in which the stiffness (Young’s modulus) and strength (failure strain) of the material vary as jointly lognormal random fields. The fracture strength of the specimen is measured by the elongation, load, and toughness at two critical stages: when fracture initiates at the notch tip and, in the 2D case, when fracture propagates through the entire specimen. This is an extension of a previous study on the elastic and fracture properties of systems with random Young’s modulus and deterministic material strength (Dimas et al., 2015a). For 1D rods our approach is analytical and builds upon the ANOVA decomposition technique of (Dimas et al., 2015b). In 2D we use a semi-analytical model to derive the fracture initiation strengths and regressions fitted to simulation data for the effect of crack arrest during fracture propagation. Results are validated through Monte Carlo simulation. Randomness of the material strength affects in various ways the mean and median values of the initial strengths, their log-variances, and log-correlations. Under low spatial correlation, material strength variability can significantly increase the effect of crack arrest, causing ultimate failure to be a more predictable and less brittle failure mode than fracture initiation. These insights could be used to guide design of more fracture resistant composites, and add to the design features that enhance material performance.     

Tailored heterogeneity increases overall stability regime of composites having a negative-stiffness inclusion

Recent theoretical and experimental results have shown the possibility of enormous increases in composite material overall elastic stiffness, damping, thermal expansion, piezoelectricity, etc., when the composite contains a tuned non-positive-definite (i.e., negative stiffness) constituent. For such composite materials to have practical utility, they must be stable. Recent research has shown they can be, for a limited range of constituent negative stiffness. This research has treated linear elastic composite materials with homogeneous phases, via the energy method and full dynamic stability analyses.In the present work, we first show how to analyze the composites previously treated by the comprehensive but simpler static stability approach, obtaining closed-form results. We then employ this approach to show that permitting heterogeneity of the positive-definite phase can substantially increase the range of constituent negative stiffness while maintaining overall composite stability. We first treat the positive-definite phase heterogeneity as piecewise homogeneous, and then treat it as continuously-varying. In the continuously-varying heterogeneity case, we seek the radially optimal distribution of the elastic moduli in the coatings, under constant coating average moduli constraint, to permit the most negative possible inclusion stiffness while maintaining overall composite stability. This is accomplished for three coating cases: constant bulk modulus but arbitrarily radially-varying shear modulus; constant shear modulus but arbitrarily radially-varying bulk modulus; and both moduli arbitrarily radially varying. We find the optimal coatings to be: a heterogeneous one with shear modulus being a specific continuously decreasing function of radius for the first case; a homogeneous one for the second case; and a heterogeneous one with both moduli being either Dirac-delta or Heaviside-step decreasing functions of radius for the last case (if the coating moduli are unrestricted in magnitude or have upper limits, respectively). The results show a substantial increase in the permissible inclusion negative stiffness range is provided by coating heterogeneity, while maintaining overall composite stability. Such an increased range of constituent negative stiffness provides an enlarged tuning parameter range for the development of novel, high-performance composite materials.     

Application of a discrete polycrystal model to the analysis of cyclic straining in copper

A discrete polycrystal model, designed to simulate a metal aggregate macro-element, is applied to the study of cyclic straining in copper. The numerical method of solution (an adaptation of the “finite element method”) incorporates a convergent discrete Green's function within the constrained minimum principle which governs the (crystallographic) plastic shear increments at each load step. Isothermal elastic moduli of copper crystals and Taylor's hardening rule with constant hardening modulus are used in the calculations. Numerical results are obtained for macroscopic elastic properties, cyclic stress-strain curves (which indicate the contribution of aggregate heterogeneity to macroscopic hardening), macroscopic plastic work, and residual (latent) strain energy through four loading cycles between fixed macrostrain limits. Other estimates for elastic properties also are included, and all results are compared, both qualitatively and quantitatively, with published experiments. The predictions of the model are in general satisfactory.     

Discontinuous bifurcation analysis in fracture energy-based gradient plasticity for concrete

Conditions for discontinuous bifurcation in limit states of selective non-local thermodynamically consistent gradient theory for quasi-brittle materials like concrete are evaluated by means of both geometrical and analytical procedures. This constitutive formulation includes two internal lengths, one related to the strain gradient field that considers the degradation of the continuum in the vicinity of the considered material point. The other characteristic length takes into account the material degradation in the form of energy release in the cracks during failure process evolution.The variation from ductile to brittle failure in quasi-brittle materials is accomplished by means of the pressure dependent formulation of both characteristic lengths as described by Vrech and Etse (2009).In this paper the formulation of the localization ellipse for constitutive theories based on gradient plasticity and fracture energy plasticity is proposed as well as the explicit solutions for brittle failure conditions in the form of discontinuous bifurcation. The geometrical, analytical and numerical analysis of discontinuous bifurcation condition in this paper are comparatively evaluated in different stress states and loading conditions.The included results illustrate the capabilities of the thermodynamically consistent selective non-local gradient constitutive theory to reproduce the transition from ductile to brittle and localized failure modes in the low confinement regime of concrete and quasi-brittle materials.     

The Three-Stage Model Based on Strain Strength Distribution for the Tensile Failure Process of Rock and Concrete Materials

A three-stage model is introduced to describe the tensile failure process of rock and concrete materials.Failure of the material is defined to contain three stages in the model,which include elastic deformation stage,body damage stage and localization damage stage.The failure mode change from uniform body damage to localization damage is expressed.The heterogeneity of material is described with strain strength distribution.The fracture factor and intact factor,defined as the distribution function of strain strength,are used to express the fracture state in the failure process.And the distributive parameters can be determined through the experimental stress-strain curve.     

二维编织陶瓷基复合材料应力-应变行为

对二维编织陶瓷基复合材料拉伸应力-应变行为进行了试验研究和理论模拟. 将二维 编织结构简化为：正交铺层结构和纤维束波动结构. 基于基体随机开裂、纤维随机断裂分布 理论,得到正交铺层结构的应力-应变关系;基于体积平均方法,将纤维束波动部分进行分割, 引入强度分析模型,得到纤维束波动部分的应力-应变关系. 结合正交铺层部分和纤维束波动 部分的应力-应变关系,得到二维编织结构的应力-应变行为,理论与试验吻合较好.     

影响冲击载荷下脆性材料碎片尺度的因素

概括回顾近年我们针对脆性材料受冲击拉伸载荷时的断裂及破碎现象的理论分析及数值模拟成果。重点分析碎片尺度对材料参数及加载速率的依赖关系。通过量纲分析,动力学模拟以及数值实验,建立了一个普适性的无量纲关系用以评估碎片尺度。将本研究结果与其他理论及现有实验结果进行的比较表明本模型更为完善合理。