Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials: A 3D study

In a recent publication (Yang et al., 2009. Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials. Int. J. Solids Struct. 46 (17) 3222–3234), we developed a finite element method capable of modelling complex two-dimensional (2D) crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. The present study extends the method to model three-dimensional (3D) problems. First, 3D cohesive elements are inserted into the initial mesh of solid elements to model potential crack surfaces by a specially designed, flexible and efficient algorithm and corresponding computer program. The softening constitutive laws of the cohesive elements are modelled by spatially-varying 3D Weibull random fields. Monte Carlo simulations are then carried out to obtain statistical information of structural load-carrying capacity. A concrete cube under uniaxial tension was analysed as an example. It was found that as the 2D heterogeneous model, the 3D model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing heterogeneity in terms of the variance in the tensile strength random fields resulted in lower mean and higher standard deviation of peak loads. Due to constraint effects and larger areas of unsmooth, non-planar fracture surfaces, 3D modelling resulted in higher mean and lower standard deviation of peak loads than 2D modelling.     

大开口复合材料层合板强度破坏研究

复合材料层合板的各向异性及非均质,使得复合材料层合板内部的破坏形式非常复杂.在复合材料结构的设计中,为满足制造及使用功能上的需求,在复合材料层合板承力结构件上不可避免地需要设计各种开口.然而,含大开口复合材料层合板的强度破坏问题变得更为复杂,使得现有的强度理论面临新的挑战.针对碳纤维增强复合材料大开口层合板受单向拉伸载荷作用下的强度破坏问题进行了数值分析和实验研究.首先,根据Hashin准则和刚度退化模型,对含不同圆形开口尺寸的[0]_(10)单向铺层、[0/90]_5和[±45]_5正交铺层的层合板,进行了单向拉伸载荷作用下渐进失效的数值模拟分析,获得了对应结构的极限载荷和破坏模式.在此基础上,采用数字图像相关方法,进行复合材料大开口层合板强度破坏的实验研究.研究结果表明,大开口复合材料层合板在单向拉伸加载下主要呈现脆性破坏形式,破坏起始位置处于应力集中区.此外,破坏强度和失效模式与复合材料铺层方式和开口尺寸大小密切相关.其中[±45]_5铺层的开口层合板承载能力最弱,分层破坏最严重.开口尺寸越大,结构的极限载荷值越低.同实验测试结果相比,数值模拟对复合材料层合板的损伤失效分析略显不足,往往很难全面分析复合材料层合板破坏失效过程中的各种因素的影响.     

一种随机裂纹结构的裂纹扩展路径分析方法

裂纹结构中存在大量不确定性因素,如裂纹长度、材料性质、外部载荷等,裂纹扩展路径的不确定性分析对研究随机裂纹结构损伤和断裂的力学特性并预测其性能及可靠性具有重要意义。本文提出了一种适应于混合载荷模式下随机裂纹结构的裂纹扩展路径分析方法。该方法考虑了裂纹长度、材料性质和外部载荷等的随机性,并通过蒙特卡洛方法对随机参数空间进行采样。采用比例边界有限元方法计算结构应力强度因子,进而模拟单次裂纹扩展路径。在此基础上,通过概率分析方法获得随机裂纹结构中裂纹扩展路径的统计特性。最后给出了两个数值算例验证了本文方法的有效性。     

Discrete fracture modeling of asphalt concrete

A heterogeneous fracture approach is presented for modeling asphalt concrete that is composed of solid inclusions and a viscous matrix, and is subjected to mode-I loading in the fracture test configuration. A heterogeneous fracture model, based on the discrete element method (DEM), is developed to investigate various fracture toughening mechanisms of asphalt materials using a high-resolution image processing technique. An energy-based bilinear cohesive zone model is used to model the crack initiation and propagation of materials, and is implemented as a user-defined model within the discrete element method. Experimental fracture tests are performed to investigate various fracture behavior of asphalt concrete and obtain material input parameters for numerical models. Also, bulk material properties are necessary for each material phase for heterogeneous numerical models; these properties are determined by uniaxial complex modulus tests and indirect tensile strength tests. The main objective of this study is to integrate the experimental tests and numerical models in order to better understand the fracture mechanisms of asphaltic heterogeneous materials. Experimental results and numerical simulations are compared at different test conditions with excellent agreement. The heterogeneous DEM fracture modeling approach has the potential capability to understand various crack mechanisms of quasi-brittle materials.     

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

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

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

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

Experimental and numerical investigation of fracture in a cast aluminium alloy

This paper describes an experimental and numerical investigation on the fracture behaviour of a cast AlSi9MgMn aluminium alloy. In the experiments, a modified Arcan test set-up was used to study mixed-mode fracture. During testing, the tension load and the displacement of the actuator of the test machine were recorded, simultaneously as a high-resolution digital camera was used to record a speckle-patterned surface of the specimen. The recorded images were post-processed using an in-house digital image correlation (DIC) software to obtain information of the displacement and strain fields in the specimen during the test. In addition, some newly implemented features in the DIC software allowed us to detect and follow the crack propagation in the material. The numerical calculations were carried out with a user-defined material model implemented in an explicit finite element code. In the model, the material behaviour is described by the classical J₂ flow theory, while fracture was modelled by the Cockcroft–Latham criterion, assuming the fracture parameter to follow a modified weakest-link Weibull distribution. With the proposed probabilistic fracture modelling approach, the fracture parameter can be introduced as a random variable in the finite element simulations. Crack propagation was modelled by element erosion, and a non-local damage formulation was used to reduce mesh-size sensitivity. To reveal the effect of mesh density and meshing technique on the force–displacement curves and the crack propagation, several different meshes were used in the numerical simulations of the modified Arcan tests. The numerical results were finally compared to the experimental data and the agreement between the measured and predicted response was evaluated.     

非常规态型近场动力学热黏塑性模型及其应用

在非常规态型近场动力学(non-ordinary state-based peridynamics, NOSB-PD) 理论框架下构建了考虑应变率效应、塑性硬化、热软化效应和材料断裂特征的非局部三维热黏塑性固体本构模型以及相应的非局部空间积分型数值算法, 并应用于金属类材料和构件在冲击载荷作用等工况下的高应变率热黏塑性变形与破坏分析. 通过对经典含初始裂纹Kalthoff-Winkler板冲击试验进行三维近场动力学模拟, 可得到裂纹的起裂角度、扩展路径、扩展速度以及裂纹扩展过程中靶板等效应力和温度分布, 所得结果与已有试验结果和其他数值方法结果吻合较好. 在此基础上, 应用该模型分析了不同冲击速度作用下金属靶板的变形与裂纹扩展过程, 结果表明: 该模型能较好地模拟不同冲击速度(应变率)情况下靶板的变形与破坏全过程. 随着冲击速度变化, 初始裂纹的起裂时间、扩展方向和扩展速度呈一定规律变化. 冲击速度越低, 起裂时间越晚(直至冲击速度低于某值时初始裂纹不扩展), 裂纹扩展速度峰值越低, 冲击过程中靶板温度峰值越低, 完全扩展所需时间越长.      

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.     

基于反向传播神经网络的疲劳裂纹扩展分析

材料发生疲劳断裂时往往会引起重大安全事故,而基于传统数值模拟方法求解疲劳裂纹扩展问题时模 型复杂、计算量大．本文基于包含多隐层的反向传播神经网络分析金属材料疲劳裂纹扩展行为,计算了裂纹扩 展过程中的 von Mises应力场和位移场,并与数值解和实验解进行对比,误差分析结果表明其求解精度高．并 基于该神经网络有效预测了裂纹扩展中裂纹长度及裂纹扩展速率的变化过程,预测精度高．该神经网络分析方 法可为材料剩余寿命和疲劳强度预测提供研究基础．     

Inverse Estimation of Cohesive Fracture Properties of Asphalt Mixtures Using an Optimization Approach

Tensile cracking in asphalt pavements due to vehicular and thermal loads has become an experimental and numerical research focus in the asphalt materials community. Previous studies have used the discrete element method (DEM) to study asphalt concrete fracture. These studies used trial-and-error to obtain local fracture properties such that the DEM models approximate the experimental load-crack mouth opening displacement response. In the current study, we identify the cohesive fracture properties of asphalt mixtures via a nonlinear optimization method. The method encompasses a comparative investigation of displacement fields obtained using both digital image correlation (DIC) and heterogeneous DEM fracture simulations. The proposed method is applied to two standard fracture test geometries: the single-edge notched beam test, SE(B), under three-point bending, and the disk-shaped compact tension test, DC(T). For each test, the Subset Splitting DIC algorithm is used to determine the displacement field in a predefined region near the notch tip. Then, a given number of DEM simulations are performed on the same specimen. The DEM is used to simulate the fracture of asphalt concrete with a linear softening cohesive contact model, where fracture-related properties (e.g., maximum tensile force and maximum crack opening) are varied within a predefined range. The difference between DIC and DEM displacement fields for each set of fracture parameters is then computed and converted to a continuous function via multivariate Lagrange interpolation. Finally, we use a Newton-like optimization technique to minimize Lagrange multinomials, yielding a set of fracture parameters that minimizes the difference between the DEM and DIC displacement fields. The optimized set of fracture parameters from this nonlinear optimization procedure led to DEM results which are consistent with the experimental results for both SE(B) and DC(T) geometries.     

螺杆钻具传动轴接头表面裂纹断裂性能研究

以螺杆钻具的作业载荷为依据,对传动轴上端接头螺纹进行断裂性能研究。首先,通过强度分析,确定接头螺纹易发生断裂的危险位置;然后,考虑材料的弹塑性以及螺纹升角,建立螺纹根部含横向半椭圆表面裂纹的外螺纹接头和相啮合水帽的有限元模型,定量评价弹塑性和螺纹升角对接头断裂性能的影响。结果表明,材料塑性变形和螺纹升角对裂纹前缘J积分的影响较大,忽略这两个因素得到的结果偏于不安全。最后,利用所建立的弹塑性有限元模型,对于上扣预紧力单独作用以及上扣预紧力和工作扭矩共同作用的两种载荷条件,得到螺纹根部裂纹前缘的临界形状比及其对应的量纲为一的J积分随裂纹扩展的变化,为进一步研究螺杆钻具传动轴接头螺纹的疲劳寿命提供依据。     

局部高密度钢纤维混凝土疲劳裂纹扩展实验研究

在混凝土弯曲构件底部用高密度钢纤维局部增强称为局部高密度钢纤维混凝土（PHPFRC）,与同样纤维掺量的传统钢纤维混凝土（SFRC）比较,它可以用近似的价格获得高得多的力学性能,本文对PHPFRC缺口试件在循环载下的裂纹扩展规律进行了研究,发现裂纹并非是由一条主裂纹循序渐进地发展,而是随着循环次数的增加,呈多条短裂纹陆续出现,弥散式发展的趋势,说明在高密度纤维的增强下,混凝土的疲劳断裂由脆性向专心性转变,根据实验及分析结果得出了两种纤维掺量下PHPFRC度件的Paris参数,并且用一个分维值D的变化描述了裂纹的演经过程,发现裂纹失稳扩展前的临界分维值D随着纤维掺量的增加而增加。     

横观各向同性油气藏水力压裂裂纹扩展规律研究

针对横观各向同性与各向同性油气藏水力压裂裂纹扩展的差异性,基于扩展有限元法建立水力压裂力学模型,通过ABAQUS子程序开发了各向同性和横观各向同性岩体的起裂判据。在各向同性岩体数值模拟结果与解析解以及现场压裂典型曲线对比吻合的基础上,得到了包含层理构造的横观各向同性岩体水力压裂过程中裂纹扩展规律。层理类岩体水力压裂的裂缝扩展方向由地应力状态、层理方向以及岩体与层理界面抗拉强度共同决定;水力压裂过程中,注水压力在裂纹尖端产生应力集中,层理面法向分量先达到界面抗拉强度时,裂纹沿层理方向开裂,反之裂纹沿垂直最小地应力的方向扩展;裂纹扩展速度随层理抗拉强度的增加而降低;由于地层的滤失,随压裂液的注入,裂纹长宽尺度增长速率降低。     

高强混凝土热-力响应的断裂相场分析方法研究

高强混凝土(High-Strength Concrete, HSC)的力学性能受温度影响较大,尤其在高温环境下性能衰退剧烈,因此在结构分析中宜采用随温度变化的材料性能参数．断裂相场方法基于Griffith变分断裂准则,无需复杂的裂纹拓展追踪技术便可处理多裂纹的问题,可方便地模拟裂纹的萌生、扩展、分岔及汇合过程．本文采用适宜进行结构断裂分析的断裂相场方法,基于能量泛函变分原理,将温度对混凝土材料弹性模量及断裂能的影响引入断裂相场分析方法中,用于高强混凝土高温环境下的强度和破坏分析．以高温作用下高强混凝土梁三点弯曲试验为算例,进行方法验证,通过与实测结果对比,证实了算法的有效性．     

混凝土拉伸软化曲线折线近似的逆解方法

研究基于Hillerborg的虚拟裂纹模型，利用有限元分析方法，求得折线近似的拉伸软化曲线的逆解方法。对弹性模量，初始开裂应力的决定方法进行了研究。以双直线模型的计算结果为算例进行了逆推分析，算例符合得很好。也较好地从实验得到的荷载位移曲线再现了拉伸软化曲线。这对于研究混凝土的断裂能，尺寸效应等问题很具意义。     

随机超载下疲劳裂纹扩展寿命的计算

针对在基本循环载荷上加入随机超载序列的疲劳裂纹扩展问题，应用裂纹闭合的概念考虑超载的迟滞效应，将延迟时间描述成纯离散型马尔可夫过程。对相应的柯尔莫哥洛夫－费勒微积分方程进行了分步求解，结合疲劳断裂分别出现在基本循环峰载作用时和超载作用时两种情况，计算出不同可靠度下裂纹的扩展寿命，并研究了超载大小和发生强度对扩展寿命的影响。     

基于Cosserat近场动力学的纤维混凝土破坏模拟

本文基于Cosserat近场动力学发展了一种纤维混凝土的近场动力学模型,对纤维混凝土的破坏现象进行研究。该模型考察了物质点独立的转动自由度和物质点间的力偶作用,而且有表征微结构尺寸的内禀长度,相比于传统的近场动力学模型,更适合描述纤维混凝土这类胶结颗粒材料的力学行为。本文采用完全离散的方式对纤维进行建模,引入了纤维拔出实验中拔出位移和切应力的关系,并且采用组构张量描述纤维的局部分布。数值算例部分模拟了单纤维拔出实验、带预制裂纹的平板拉伸实验和三点弯曲梁实验。数值结果和已有的数值模型以及实验进行了对比,验证了所提出模型的正确性。此外,本文还调查了微结构对纤维混凝土破坏的影响,数值结果显示Cosserat剪切模量和内禀长度会影响裂纹的局部分布,但是不会改变裂纹的主方向。     

单轴荷载下饱水岩石静态和动态抗压强度的细观力学分析

由于单轴荷载下饱水岩石的动态力学特性与静态力学特性存在很大差异,从宏观上进行力学分析 存在局限性。根据岩石受压全应力应变曲线的细观机制,分析了静态及动态单轴荷载条件下孔隙水影响饱水 岩石裂纹扩展的情况。在静态单轴压缩条件下,初始裂隙受压使自由水产生孔隙水压力,自由水对翼裂纹有 向外挤压的应力,促进裂纹扩展。在动态单轴压缩条件下,自由水会产生粘结力,抑制裂纹扩展。根据翼裂纹 受压扩展原理,推导出饱水单轴条件下动态抗压强度、静态抗压强度的计算公式,在相同断裂韧度下,饱水岩 石静态抗压强度风干岩石静态抗压强度饱水岩石动态抗压强度。对自然风干和饱水砂岩进行单轴静态、 动态压缩实验,结果与理论模型的结果相符。     

复合载荷下脆性材料紧凑拉伸试件疲劳裂纹的预制

本文针对脆性材料的紧凑拉伸（CT）试件,提出了在复合疲劳载荷作用下预制疲劳裂纹的新方法,用此方法预制出7060-T6铝合金的CT试件疲劳裂纹并测试其断裂韧性。由于采用了框架结构和螺旋结构在CT试件的两侧面施加压缩载荷,在拉伸疲劳载荷作用下产生与之方向相反的压缩疲劳载荷,形成复合疲劳载荷,引起应力强度因子范围ΔK的下降,导致裂纹扩展速率的降低,从而控制了疲劳预制裂纹长度。此压缩疲劳载荷值用贴电阻片方法和组成全桥原理测定出。在一定拉伸疲劳载荷和压缩载荷下,可成功地预制出所需的疲劳裂纹长度。