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.     

Disk-shaped compact tension test for asphalt concrete fracture

A disk-shaped compact tension (DC(T)) test has been developed as a practical method for obtaining the fracture energy of asphalt concrete. The main purpose of the development of this specimen geometry is the ability to test cylindrical cores obtained from in-place asphalt concrete pavements or gyratory-compacted specimens fabricated during the mixture design process. A suitable specimen geometry was developed using the ASTM E399 standard for compact tension testing of metals as a starting point. After finalizing the specimen geometry, a typical asphalt concrete surface mixture was tested at various temperatures and loading rates to evaluate the proposed DC(T) configuration. The variability of the fracture energy obtained from the DC(T) geometry was found to be comparable with the variability associated with other fracture tests for asphalt concrete. The ability of the test to detect changes in the fracture energy with the various testing conditions (temperature and loading rate) was the benchmark for determining the potential of using the DC(T) geometry. The test has the capability to capture the transition of asphalt concrete from a brittle material at low temperatures to a more ductile material at higher temperatures. Because testing was conducted on ungrooved specimens, special care was taken to quantify deviations of the crack path from the pure mode I crack path. An analysis of variance of test data revealed that the prototype DC(T) can detect statistical differences in fracture energy resulting for tests conducted across a useful range of test temperatures and loading rates. This specific analysis also indicated that fracture energy is not correlated to crack deviation angle. This paper also provides an overview of ongoing work integrating experimental results and observations with numerical analysis by means of a cohesive zone model tailored for asphalt concrete fracture behavior.     

Experiment and simulation of slurry flow in irregular channels to understand proppant transport in complex fractures

Slurry flow and proppant placement in irregular fractures are crucial to evaluate hydraulic fracturing stimulation but need to be better understood. This study aims to investigate how irregular fracture affects proppant transport and distribution using laboratory experiments and micro-scale numerical models. The unresolved method of the computational fluid dynamics (CFD) and the discrete element method (DEM) considers Saffman lift force, Magnus force, and virtual mass force to accurately capture the frequent interaction between proppant and slickwater. Experimental results validated the reliability of the optimized CFD-DEM model and calibrated primary parameters. The effects of crack height and width, bending angle, and distance between the crack and inlet on particle distribution were studied. The results indicated that the improved numerical method could rationally simulate proppant transport in fractures at a scale factor. The small crack height causes downward and upward flows, which wash proppant to the fracture rear and form isolated proppant dunes. A wider region in the fracture is beneficial to build up a large dune, and the high dune can hinder particle transport into the fracture rear. When the crack is close to the inlet, the primary fracture without proppants will close to hinder gas production. The smaller the bending angle, the smaller the proppant dune. A regression model can precisely predict the dune coverage ratio. The results fundamentally understand how complex fractures and natural cracks affect slurry flow and proppant distribution.     

结构破坏的尺度律

文中综述了结构破坏的尺度律和尺寸效应的研究进展,尤其将重点放在准脆性材料的分析上,因为它们的尺寸效应是重要和复杂的．在回顾了尺寸效应研究的悠远发展史以后。着重讨论了三种主要类型的尺寸效应,即由于强度随机性引起的统计尺寸效应、能量释放的尺寸效应和由于微裂纹或断裂的分形特性可能引起的尺寸效应．得出了这些理论应用的明确结论．之后讨论了如何运用已知的尺寸效应律来测量材料的断裂特性,并采用内聚裂纹模型（cohesivecrackmodel）、非局域化有限元模型和离散元模型等对尺寸效应进行模化．文中还进而分析了尺寸效应在压缩失效和车相关材料行为下的有关问题,并讨论了在断裂扩展区描述含微裂纹材料所需的损伤本构关系．最后也讨论了尺寸效应对准脆性材料的多种应用,这些材料包括,如混凝土、海冰、纤维复合材料、岩石和陶瓷等．本文包含了参考文献377篇     

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.     

基于离散元法的涂层硬质合金刀具建模及其加工损伤预测

首先建立了TiCN涂层硬质合金刀具基体材料(WC)的离散元模型,根据单轴压缩、三点弯曲以及断裂韧性等数值试验方法校准了基体材料离散元模型的微观参数,然后采用划痕法校准了基体与涂层的界面结合强度。根据Merchant切削模型,建立了涂层刀具切削过程中的刀-屑接触模型,通过对切屑施加周期性边界条件来模拟实际的切削加工过程;模拟了涂层刀具加工过程中的裂纹扩展和破坏情况,并预测了切削加工用量对涂层裂纹扩展及破坏的影响。     

Modeling of ductile fracture: Application of the mechanism-based concepts

This paper summarizes our recent studies on modeling ductile fracture in structural materials using the mechanism-based concepts. We describe two numerical approaches to model the material failure process by void growth and coalescence. In the first approach, voids are considered explicitly and modeled using refined finite elements. In order to predict crack initiation and propagation, a void coalescence criterion is established by conducting a series of systematic finite element analyses of the void-containing, representative material volume (RMV) subjected to different macroscopic stress states and expressed as a function of the stress triaxiality ratio and the Lode angle. The discrete void approach provides a straightforward way for studying the effects of microstructure on fracture toughness. In the second approach, the void-containing material is considered as a homogenized continuum governed by porous plasticity models. This makes it possible to simulate large amount of crack extension because only one element is needed for a representative material volume. As an example, a numerical approach is proposed to predict ductile crack growth in thin panels of a 2024-T3 aluminum alloy, where a modified Gologanu–Leblond–Devaux model [Gologanu, M., Leblond, J.B., Devaux, J., 1993. Approximate models for ductile metals containing nonspherical voids – Case of axisymmetric prolate ellipsoidal cavities. J. Mech. Phys. Solids 41, 1723–1754; Gologanu, M., Leblond, J.B., Devaux, J., 1994. Approximate models for ductile metals containing nonspherical voids – Case of axisymmetric oblate ellipsoidal cavities. J. Eng. Mater. Tech. 116, 290–297; Gologanu, M., Leblond, J.B., Perrin, G., Devaux, J., 1995. Recent extensions of Gurson’s model for porous ductile metals. In: Suquet, P. (Ed.) Continuum Micromechanics. Springer-Verlag, pp. 61–130] is used to describe the evolution of void shape and void volume fraction and the associated material softening, and the material failure criterion is calibrated using experimental data. The calibrated computational model successfully predicts crack extension in various fracture specimens, including the compact tension specimen, middle crack tension specimens, multi-site damage specimens and the pressurized cylindrical shell specimen.     

DEM与FEM动态耦合算法研究

离散单元法作为一种有效的数值分析方法,能够模拟脆性材料的裂纹扩展及碎片飞散等破坏特性,但是无法从根本上克服计算效率低下的诟病;传统有限单元法具有计算高效稳定的优点,却难以描述脆性材料冲击破坏过程中连续体向非连续体的转化。本文首先提出一种基于罚函数法的改进型离散单元和有限单元耦合方法,以提高耦合分析精度。在此基础上提出了动态耦合算法:即在初始阶段,模型全部为有限单元,当局部即将发生破坏时,仅使即将发生破坏的有限单元及相邻单元自动转化为离散单元,在离散单元区域研究破坏问题。这种算法充分利用有限单元法计算高效的优点,同时最大限度克服了离散单元法计算效率的不足。最后,通过两个简单算例验证了改进型耦合算法和动态耦合算法的有效性。     

Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially-varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.     

A numerical study on the deformation and fracture modes of steel projectiles during Taylor bar impact tests

A heterogeneous material model based on macro-mechanical observations is proposed for simulation of fracture in steel projectiles during impact. A previous experimental study on the deformation and fracture of steel projectiles during Taylor bar impact tests resulted in a variety of failure modes. The accompanying material investigation showed that the materials used in the impact tests were heterogeneous on scales ranging from microstructure as investigated with SEM to variation in fracture strains from tensile tests. A normal distribution is employed to achieve a heterogeneous numerical model with respect to the fracture properties. The proposed material model is calibrated based on the tensile tests, and then used to independently simulate the Taylor bar impact tests. A preliminary investigation showed that the simulations are sensitive to assumptions regarding the anvil behaviour and friction properties. A flexible anvil and a yield-limited friction law are shown to be necessary to correctly reproduce the experimental behaviour. The proposed model is then shown to be capable of correctly reproducing all fracture modes but one, and also predict critical impact velocities for projectile fracture with reasonable accuracy. Fragmentation at velocities above the critical velocity is not well reproduced due to excessive element erosion. Measures to make the element erosion process more physical are proposed and discussed with their respective drawbacks. The use of a simple fracture criterion in combination with an element erosion technique accentuates the effect of distributing the fracture parameter.     

颗粒材料多尺度分析的连接尺度方法

本文提出了耦合细尺度上基于离散颗粒集合体模型的离散单元法（DEM）和粗尺度上基于Cosserat连续体模型的有限元法（FEM）的连接尺度方法（BSM）以研究颗粒材料的力学行为。采用Cosserat连续体模型和FEM模拟的粗尺度域覆盖全域,而采用离散颗粒集合体模型的DEM模拟的细尺度域仅限于需特别关注材料微结构演变和非连续变形行为的局部区域。对这两个区域间的界面提出了适当的界面条件及其实施方案。通过采用适当的连接尺度投影算子,空间离散的粗、细尺度耦合系统多尺度运动方程具有解耦和允许分别求解、因而也允许分别采用不同时间步长对粗、细尺度计算的特点,可极大地提高BSM的计算效率。文中二维地基数值算例结果说明了所陈述方法的可应用性,以及相对基于Cosserat连续体模型的FEM和基于离散颗粒集合体模型的DEM的优越性。      

A Bridging Scale Method for Multi-scale Analysis of Granular Materials

本文提出了耦合细尺度上基于离散颗粒集合体模型的离散单元法（DEM）和粗尺度上基于Cosserat连续体模型的有限元法（FEM）的连接尺度方法（BSM）以研究颗粒材料的力学行为。采用Cosserat连续体模型和FEM模拟的粗尺度域覆盖全域,而采用离散颗粒集合体模型的DEM模拟的细尺度域仅限于需特别关注材料微结构演变和非连续变形行为的局部区域。对这两个区域间的界面提出了适当的界面条件及其实施方案。通过采用适当的连接尺度投影算子,空间离散的粗、细尺度耦合系统多尺度运动方程具有解耦和允许分别求解、因而也允许分别采用不同时间步长对粗、细尺度计算的特点,可极大地提高BSM的计算效率。文中二维地基数值算例结果说明了所陈述方法的可应用性,以及相对基于Cosserat连续体模型的FEM和基于离散颗粒集合体模型的DEM的优越性。     

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

水力劈裂问题的态型近场动力学建模

将态型近场动力学理论引入水力劈裂问题的模拟。构建了能反映岩土类材料准脆性断裂特征的态型近场动力学本构模型,并在物质点间相互作用力模型中加入等效水压力项,以实现在新生裂纹面上跟踪施加水压力。同时,考虑裂纹面间的接触,引入物质点间的短程排斥力作用,并设计了相应的接触算法。通过自编程序将模型和算法应用于含初始裂纹、不含初始裂纹以及含坝基软弱结构面的混凝土重力坝在高水头作用下的水力劈裂过程模拟,并与扩展有限元等模拟结果对比,验证了本文模型和算法的可行性和准确性。     

混凝土断裂过程的力学模型与数值模拟

混凝土是一种典型的非均匀准脆性材料,其断裂过程是非常复杂的.评述了混凝土断裂以及断裂过程研究的进展状况,主要总结了一些研究混凝土断裂的宏观力学方法和模型,介绍了从混凝土的细观层次出发模拟其断裂过程的数值模型.最后,简单介绍了作者用细观数值方法模拟混凝土宏观断裂过程的研究成果.      

基于离散元的杆件断裂全过程数值模拟

杆件的断裂会涉及到大变形、非线性以及不连续等问题,通常的数值计算方法模拟这种复杂力学行为具有局限性。本文基于颗粒离散元法DEM,将接触粘结处的分布式弹簧用梁纤维进行等效,提出了一种适于结构弹塑性分析的DEM纤维梁模型,然后在此基础上构建了构件断裂模拟算法以及纤维破环准则。将该模型应用于悬臂梁结构,模拟了悬臂梁从弹性到弹塑性阶段,再到断裂破坏的全过程,数值模拟得到的结构响应和截面开裂破坏形态均较合理。最后将该方法应用于单层网壳倒塌破坏模拟,并与网壳振动台倒塌试验进行对比,结果表明,数值模拟得到的杆件断裂过程及结构倒塌模式与试验现象一致,验证了该模型的正确性和适用性。     

粗粒土离散元非球趋真颗粒模型

颗粒形状对粗粒土的物理力学特性有着显著的影响。离散元法广泛应用于粗粒土宏观物理力学特性的细观机理研究。为了考虑颗粒形状的影响,亟待发展计算高效的离散元非球趋真颗粒模型。本文基于X射线CT扫描技术并结合数字图像处理技术,对光滑和棱角性两类典型粗粒土（鹅卵石与碎石）进行三维重构,并提出了两类趋真颗粒模型,分别采用扩展超椭球模型和球多面体模型进行趋真逼近;开展了两类颗粒试样的3D打印和单轴压缩试验,分析了配位数和局部孔隙率分布等细观特性;基于离散元开源程序SudoDEM开展了两类试样的离散元模拟,并将模拟细观分析结果与物理试验进行了对比。结果表明,提出的两类趋真颗粒模型能够较好地对粗粒土颗粒进行离散元建模。     

脆性材料破坏过程分析的数值试验方法

文中运用作者新近开发的材料破坏过程分析MFPA2D系统,以岩石、混凝土等非均匀脆性材料的破坏过程分析为例,说明了数值模拟方法给脆性材料破坏理论发展所带来的契机,并简述了MFPA2D在煤层移动、地下工程稳定性、地震孕育机制,以及复合材料破坏问题研究中的应用前景．