基于一类非局部宏-微观损伤模型的裂纹模拟

结合近场动力学和统一相场理论的基本思想,最近提出了一类非局部宏-微观损伤模型,为固体裂纹扩展模拟提供了新途径.本文在此基础上改进了微观损伤准则,并给出损伤的■语言以刻画固体破坏过程中位移场的不连续程度.在改进模型中,首先根据两物质点(即物质点对)之间的变形量,基于相对临界伸长量的历史最大超越程度,给出表征物质键性能退化的微细观损伤.进而,对影响域内的物质键损伤进行空间局部加权平均,获得宏观拓扑损伤.通过引入能量退化函数,建立基于能量的损伤与宏观拓扑损伤之间的关系,由此将其嵌入连续损伤力学基本框架,形成了问题求解的基本方程.该模型是一类非局部化模型,可采用有限单元法进行离散求解,避免了经典局部损伤力学所面临的网格敏感性问题.文中,进一步将其应用于具有强非线性回弹特性的裂纹扩展模拟问题.实例分析表明,本文方法不仅可以把握裂纹扩展模式,而且能够定量刻画裂纹扩展过程中的载荷-变形关系.最后指出了需要进一步研究的问题.     

基于一类非局部宏-微观损伤模型的混凝土典型试件力学行为模拟

混凝土是一类典型的准脆性材料, 其受力过程中的非线性分析与裂纹模拟依然是具有挑战性的问题. 经典的断裂力学与损伤力学分别从间断与连续的视角对裂纹拓扑进行了描述, 是早期人们研究固体破坏问题的有力工具. 21世纪以来, 相场理论和近场动力学在预测裂纹的萌生、扩展与非线性分析方面取得了重要的进展. 最近, 结合统一相场理论与近场动力学的基本思想, 发展了一类非局部宏-微观损伤模型. 该模型引入物质点偶的概念来刻画由于变形引起的微细观损伤, 对微细观损伤在作用域中进行加权平均得到定量描述物质不连续程度的拓扑损伤. 通过具有物理机制的能量退化函数, 将拓扑损伤嵌入到连续介质-损伤力学的框架中, 这使得该模型在进行非线性分析的同时可以自然地进行裂纹模拟, 而毋须预设初始裂纹与裂纹扩展路径. 本文考虑细观物理参数的空间变异性, 采用非局部宏-微观损伤模型进行混凝土试件受力全过程的精细化模拟. 通过一维建模标定模型细观参数, 并探讨了细观参数与混凝土材料细观物理-几何特性之间的内在关联, 在此基础上采用二维模型进行精细化分析. 进而, 考察了材料参数空间变异性对混凝土单轴受拉试件和带缺口三点弯曲试件力学行为的重要影响. 本文的研究工作为非局部宏-微观损伤模型细观参数的试验标定与复杂应力状态下混凝土等准脆性材料的非线性力学行为研究提供了有意义的参考.      

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

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

基于结构化变形驱动的非局部宏-微观损伤模型的真Ⅱ型裂纹模拟

Ⅱ型载荷作用下裂纹变形模式也为Ⅱ型的破坏问题称为真Ⅱ型破坏.准确定量地把握真Ⅱ型破坏的全过程是具有挑战性的问题.本文采用结构化变形驱动的非局部宏-微观损伤模型对真Ⅱ型破坏问题进行了模拟.根据结构化变形理论将点偶的非局部应变分解为弹性应变与结构化应变两部分,进而利用Cauchy-Born准则与结构化应变计算点偶的结构化正伸长量.在本文中,结构化应变取为非局部应变的偏量部分.当点偶的结构化正伸长量超过临界伸长量时,微细观损伤开始在点偶层次发展.将微细观损伤在作用域中进行加权求和得到拓扑损伤,并通过能量退化函数将其嵌入到连续介质-损伤力学框架中进行数值求解.进一步地,本文采用Gauss-Lobatto积分格式计算点偶的非局部应变,将积分点数目降低到4个,显著降低了前处理和非线性分析的计算成本.通过对Ⅱ型加载下裂尖应变场的分析揭示了采用偏应变作为结构化应变的原因.对两个典型真Ⅱ型破坏问题的模拟结果表明,本文方法不仅可以把握Ⅱ型加载下的真Ⅱ型裂纹扩展模式,同时可以定量刻画加载过程中的载荷-变形曲线,且不具有网格敏感性.最后指出了需要进一步研究的问题.     

准脆性材料损伤破坏的细宏观联结分析

准脆性材料在外力作用下不仅产生应力与变形,而且还会内部缺陷或微裂纹随着变化和出现宏观裂缝扩展.根据材料的微细观测,通过固体力学原理建立模型分析其微裂纹与宏观裂缝之间的关联.材料中宏观裂缝端部通常存在损伤区,且在该区内存在许多相对微小缺陷或微裂纹,从而使得材料局部变形增加而力量减弱.若把该损伤带视为裂缝一部分,那么在该虚拟裂缝两岸上将存在分布黏聚力;虚拟裂缝两岸的相对位移将是微裂纹区变形及扩展的综合体现.为此,对包含有细观裂纹的代表性体积单元进行力学计算,然后将其嵌入到宏观裂缝端部损伤带的变形计算中,以探寻固体失效的细观与宏观尺度力学分析的联结,并与实测数据相对比.     

考虑应变率效应的混凝土随机损伤本构模型研究

混凝土材料组分复杂且具有随机分布的特点, 其受力力学行为不可避免地存在非线性和随机性. 同时, 在动力荷载作用下, 混凝土材料具有不可忽视的率敏感性. 为了综合反映混凝土受力力学行为中的非线性、随机性与率敏感性, 本文从对材料的纳-微观裂纹扩展分析入手, 引入速率过程理论描述纳观裂纹的扩展速率, 并研究了对应的能量耗散过程. 在此基础上通过裂纹层级模型将纳观分析推演到微观尺度, 建立了微观能量耗散的基本表达式. 进而与微-细观随机断裂模型相结合, 形成了混凝土纳-微-细观随机损伤本构模型. 同时, 基于速率相关势垒的分析, 揭示了动力强度的提高源自加载速率和原子键断裂速率的竞争机制. 据此, 假定微裂纹间相互作用与应变率比值的相关关系以建立微弹簧能量耗散速率与应变率的联系, 实现了从静力本构模型向动力本构模型的扩展. 数值算例表明, 建议模型能够同时反映混凝土材料力学行为中的非线性、随机性和率敏感性. 最后通过与相关试验结果的对比, 验证了建议模型的正确性.      

复合材料刚度性能表征的协同损伤力学模型

针对复合材料层合板的弥散型损伤,提出一个刚度性能表征的协同损伤力学模型. 该模型兼顾了微观物理损伤响应和宏观材料刚度性能表征. 从微观角度,建立细观RVE 模型求解裂纹表面张开位移和滑开位移,以此定义损伤张量,并在宏观上通过对材料应变和损伤表面位移进行均匀化处理,建立单向板或层合板的损伤刚度矩阵和损伤张量之间的联系. 以基体裂纹为例,详细分析并建立了横向裂纹和纵向裂纹的损伤本构. 计算了[±θ/90₄]_S 铺层层合板中基体横向裂纹对刚度性能的影响,结果表明该方法能够准确地预测复合材料层合板由损伤导致的刚度性能衰减.      

基于能量耗散分析的混凝土随机疲劳损伤模型

基于微-细观随机断裂模型(MMSF)发展了一类混凝土随机疲劳损伤模型。该模型将MMSF中的微弹簧视为一能量耗散单元并考察了其跨尺度的能量耗散过程。在纳观尺度,引入速率过程理论描述裂纹的扩展速度,并基于裂纹层级模型和自相似假定完成从纳观尺度到微观尺度的过渡。由此建立了微弹簧在疲劳荷载作用下的多尺度耗能描述。此外,为了考虑疲劳加载中混凝土多条裂纹的相互影响,引入了包含损伤扩展效应和损伤愈合效应的疲劳损伤因子来修正耗能表达。通过与相关试验的对比,证明了该模型能反映疲劳荷载作用下混凝土的主要力学行为,如疲劳损伤三阶段特点、疲劳寿命的离散性和疲劳寿命随加载频率的变化趋势等。     

脆性岩石蠕变裂纹成核宏细观力学机理研究

脆性岩石内部细观裂纹扩展、贯通及成核影响下的脆性蠕变行为, 对深部地下工程围岩微震及岩爆事件评价有着重要意义. 然而, 目前能够解释裂纹成核损伤突变影响下, 脆性岩石完整蠕变宏细观力学机理模型的研究很少. 本文基于脆性岩石亚临界裂纹扩展模型、裂纹-应变-声发射事件相关的损伤模型及裂纹成核损伤时间演化路径函数, 提出了一种脆性岩石裂纹成核损伤突变影响下的蠕变宏细观力学模型. 裂纹成核损伤时间演化路径函数通过岩石内部裂纹成核损伤突变大小$\Delta D_{CN}$及相邻裂纹成核损伤时间差$\Delta t$进行定义, 该函数可以结合岩石声发射监测试验数据定义的岩石损伤数据确定. 通过与试验结果对比分析验证模型的合理性. 并讨论了裂纹成核损伤大小、相邻裂纹成核损伤时间间隔、及裂纹成核数量对脆性岩蠕变裂纹长度、裂纹速率、轴向应变及应变率的影响. 该结果对于更加合理、经济、高效的深部地下工程施工及设计提供了一定的理论支持.      

脆性固体变形响应与裂纹扩展的同步试验观测及其定量分析

通过对砂浆试样的材料变形响应和表面裂纹图象的同步观测实验研究、试样表面细观裂纹萌生和裂纹扩展图象信息的采集和量化方法的研究，探讨脆性固体材料的损伤状态与其宏观力学响应之间的关联。用本文提出的实验及其分析处理方法，可以得到砂浆、混凝土、岩石一类脆性固体材料在承载和变形过程中，表面裂纹扩展的定量结果。试样的表面裂损度代表了砂浆试样在承载和变形过程中表面裂纹扩展的情况，也在一定程度上反映了试样内部的损伤情况。     

Why continuum damage is nonlocal: Justification by quasiperiodic microcrack array

Strain-softening damage due to distributed cracking is modeled by an elastic continuum with a quasiperiodic array of cracks of regular spacing but varying sizes. As a model for the initial stage, the cracks are penny-shaped and small compared to their spacing, and as a model for the terminal stage the uncracked ligaments between the cracks are circular and small compared to their spacing. The strain due to cracks and the compliance per crack are calculated. The cracked material is homogenized in such a manner that the macroscopic continuum strains satisfy exactly the condition of compatibility with the actual strains due to cracks, and the macroscopic continuum stress satisfies exactly the condition of work equivalence with the actual stresses in the cracked material. The results show that, contrary to the existing theories, the damage variable used in continuum damage mechanics should be nonlocal, while the elastic part of the response should be local. In particular, the nonlocal continuum damage should be considered as a function of the spatial average of the cracking strain rather than its local value. The size of the averaging region is determined by the crack spacing.     

Crack shielding and material deterioration in damaged materials: an antiplane shear fracture problem

Summary A simple damage evolution model is proposed for a quasibrittle material in the framework of continuum damage mechanics. The model is used to obtain a closed form solution for a mode-III stationary crack under small scale damage conditions. It is found that the crack tip stress intensity factor is reduced, i.e., the crack is shielded by the damage. However, this shielding effect is completely offset by the material deterioration caused by the damage. It also holds for steady state crack growth. When the most effective shielding is reached for the stationary crack, the zone dominated by the stress intensity factor shrinks to the crack tip. The results for the antiplane shear problem should shed some light on the in- plane fracture problem. Received 4 August 1997; accepted for publication 7 October 1997     

Adaptive coupling between damage mechanics and peridynamics: A route for objective simulation of material degradation up to complete failure

The objective (mesh-independent) simulation of evolving discontinuities, such as cracks, remains a challenge. Current techniques are highly complex or involve intractable computational costs, making simulations up to complete failure difficult. We propose a framework as a new route toward solving this problem that adaptively couples local-continuum damage mechanics with peridynamics to objectively simulate all the steps that lead to material failure: damage nucleation, crack formation and propagation. Local-continuum damage mechanics successfully describes the degradation related to dispersed microdefects before the formation of a macrocrack. However, when damage localizes, it suffers spurious mesh dependency, making the simulation of macrocracks challenging. On the other hand, the peridynamic theory is promising for the simulation of fractures, as it naturally allows discontinuities in the displacement field. Here, we present a hybrid local-continuum damage/peridynamic model. Local-continuum damage mechanics is used to describe “volume” damage before localization. Once localization is detected at a point, the remaining part of the energy is dissipated through an adaptive peridynamic model capable of the transition to a “surface” degradation, typically a crack. We believe that this framework, which actually mimics the real physical process of crack formation, is the first bridge between continuum damage theories and peridynamics. Two-dimensional numerical examples are used to illustrate that an objective simulation of material failure can be achieved by this method.     

各向异性材料疲劳损伤模拟

根据连续损伤理论，考虑到了损伤能释放率等引起损伤的重要因素，提出了一个各向异性疲劳损伤模型。此模型结构简单，适用于各种各向异性材料。利用该模型对玻璃聚酯复合材料层板在单向受力情况下进行了寿命预测，预测值与实验值符合较好     

Fractal effect and anisotropic constitutive model for concrete

An elasto-anisotropic damage constitutive model for concrete is developed in this work. Disregarding the coupling between the isotropic and the anisotropic damage, the isotropic damage variables are defined as functions of the microcrack fractal dimension, and the anisotropic parts are expressed by the lengths of cracks in concrete which various in different directions. The Helmholtz free energy is decomposed into the elastic deforming, damage and irreversible deforming components, with the last component used to replace the plastic deformation. Therefore the damage constitutive formulas for concrete are derived based on continuum damage mechanics. Evolution laws for both isotropic and anisotropic damage variables are derived, in which the anisotropic parts are obtained by modifying an empirical model. The critical fracture stress and the fracture toughness are investigated for materials with a single fractal crack based on the fractal geometry and the Griffith fracture criterion. Numerical computation is conducted for concrete under the uniaxial and the biaxial compression. The results indicate that the material stiffness degradation can be well addressed when the anisotropic damage is incorporated; the irreversible deformation is greatly related to the behavior of the descending branch beyond the peak load. The validation of the presented model is proofed by comparing results with the experimental data. This model provides an approach to link the macro properties of a material with its micro-structure change.     

A fracture damage model for Jointed rock masses and its application to the stability analysis of dam abutments

Based on continuum damage mechanics, for jointed rock masses, a fracture damage model is presented in this paper. First, the damage tensors are defined through the elastic-flexibility of intact rock and the equivalent elastic-damage flexibility for rock mass. Then, by the self-consistent principle of solid mechanics, the equivalent elastic-damage flexibility tensors involving the interaction between multicracks are deduced. The damage evolution law is proposed involving the mechanism of crack propagation process: frictional sliding, crack kinking, growing of branched tension cracks, interlinking of the microcracks near branched crack tips leading to the breakthrough of macro-cracks and finally the failure of rock mass. Thus the evolution of damage variables reasonably unified with the process of crack propagation is given. Finally, a plastic-brittle damage constitutive relation including brittle coupled strain rate, developed and applied to the stability analysis of complicated rock foundation of a dam in China, is described in this paper.     

一个新的普遍形式的局部损伤断裂准则

利用损伤函数概念，建立了一个普遍形式的局部断裂准则。该准则考虑了局掊应力，应变和损伤历史对断裂的影响，根据损伤力学理选取了一个新的连续损伤函数，从而导出一个新的连续损伤断裂准则。新的临界断裂参数ＷＤＣ，具有明显的物理意义，且易通过试验测得，是一个不依赖于应力状态的材料九。文中还从细观力学理论和有关的试验资料出发，选取了相应的损伤函数，再现了前人的细观力学准则和经验准则。