基于单元破裂的岩石裂纹扩展模拟方法

传统离散元方法在处理破裂问题时, 采用界面上的准则进行判断, 裂纹只能沿着单元边界扩展. 当物理问题存在宏观或微观裂隙时, 在界面上应用准则具有其合理性; 而裂纹沿着单元边界扩展, 使得裂纹路径受网格影响较大, 扩展方向受到限制. 针对上述情况, 可以基于单元破裂的方式, 构建连续- 非连续单元法, 并应用于岩石裂纹扩展问题的模拟. 该方法在连续计算时, 将单元离散为具有物理意义的弹簧系统, 在局部坐标系下由弹簧特征长度、面积求解单元变形和应力, 通过更新局部坐标系和弹簧特征量, 可进一步计算块体大位移、大转动, 连续问题计算结果与有限元一致, 同时提高了计算效率. 在此基础上, 引入最大拉应力与莫尔—库伦的复合准则, 判断单元破裂状态和破裂方向, 并采用局部块体切割的方式, 在单元内形成初始裂纹. 裂纹两侧相应增加新的计算节点, 同时引入内聚力模型描述裂纹两侧的法向、切向作用与张开度及滑移变形之间的关系. 按此方式, 裂纹尖端处的扩展路径可穿过单元内部和单元边界, 在扩展方向的选取上更为准确. 最后, 通过三点弯曲梁、单切口平板拉伸、双切口试样等典型数值试验, 模拟裂纹在拉伸、压剪等各种应力状态下的扩展问题, 并对岩石单轴压缩试验的破坏过程进行模拟, 分析裂纹形成与应力—应变曲线各阶段之间的对应关系. 结果表明: 连续—非连续单元法通过单元内部破裂的方式, 可以显示模拟裂纹萌生、扩展、贯通直至形成宏观裂缝的过程.      

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.     

Experimental and numerical analysis of fretting crack formation based on 3D X-FEM frictional contact fatigue crack model

Nowadays, numerical simulation of 3D fatigue crack growth is easily handled using the eXtended Finite Element Method coupled with level set techniques. The finite element mesh does not need to conform to the crack geometry. Most difficulties associated to complex mesh generation around the crack and the re-meshing steps during the possible propagation are hence avoided. A 3D two-scale frictional contact fatigue crack model developed within the X-FEM framework is presented in this article. It allows the use of a refined discretization of the crack interface independent from the underlying finite element mesh and adapted to the frictional contact crack scale. A stabilized three-field weak formulation is also proposed to avoid possible oscillations in the local solution linked to the LBB condition when tangential slip is occurring. Two basic three-dimensional numerical examples are presented. They aim at illustrating the capacities and the high level of accuracy of the proposed X-FEM model. Stress intensity factors are computed along the crack front. Finally an experimental 3D ball/plate fretting fatigue test with running conditions inducing crack nucleation and propagation is modeled. 3D crack shapes defined from actual experimental ones and fretting loading cycle are considered. This latter numerical simulation demonstrates the model ability to deal with challenging actual complex problems and the possibility to achieve tribological fatigue prediction at a design stage based on the fatigue crack modeling.     

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

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

面内裂纹扩展问题的虚拟裂纹闭合法和离散内聚力模型

论文介绍了处理面内裂纹扩展问题的两种断裂单元,分别为处理线性问题的VCCT哑节点断裂单元和处理非线性问题的DCZM断裂单元.并与有限元软件ABAQUS的断裂单元相对比,以验证其计算准确性和对网格、步长等不敏感性,以及计算效率高等特点.     

模拟裂纹扩展的一种有限元局部动态子划分方法

提出了一种有限元子划分结合子结构的方法来模拟裂纹扩展问题。提出的方法中,将单元分为三类:被裂纹贯穿的单元,包含裂尖的单元和常规单元。对前两类单元进行子划分,每个单元的归类随裂纹的扩展而动态变化。覆盖一条裂纹的前两类单元子划分后构成一个子结构,子结构也是动态的,跟随裂纹的扩展而逐步扩大。本文的方法可以使裂纹沿任意路径扩展而不受初始网格的限制,裂纹扩展后无需对结构整体的网格重划分,结构整体分析的总自由度也不变。用该方法计算无限大平面中心裂纹的应力强度因子,模拟三点弯梁跨中裂纹的扩展,验证了计算精度,并进一步用该方法模拟了非均质材料中裂纹的扩展,考核了对复杂裂纹扩展问题的适用性。     

混凝土三点弯曲梁裂缝断裂全过程数值模拟研究

考虑裂缝黏聚力的作用,基于Paris位移公式推导出混凝土三点弯曲梁裂缝扩展过程中断裂过程区上的裂缝张开位移的解析表达式.采用起裂韧度作为裂缝起裂及扩展的判断标准,提出了荷载作用下混凝土裂缝起裂、扩展及失稳破坏全过程的数值模拟方法,并分别与国内外断裂试验实测值及有限元计算值进行了比较.结果表明,本文提出的数值模拟方法形式简单且精度较好.     

粘聚裂纹扩展的强化有限元h型网格自适应模拟

非连续变形分析和非规则节点处理是基于单元细划的粘聚裂纹扩展网格自适应模拟的关键。首先,利用强化有限单元法中数学单元和物理单元分离的特点,通过引入过渡单元,将适用于非连续变形描述的数学模式覆盖法和方便处理非规则节点的物理模式重构法结合,提出了强化有限单元法的统一关联法则,并导出了相应的单元列式。其次,基于数学裂纹尖端影响域和裂尖单元尺寸,提出了基于强化有限单元法的粘聚裂纹扩展过程模拟的h型网格自适应策略。最后,通过两个算例验证了本文方法的合理性和有效性。     

数学网格和物理网格分离的有限单元法(II):粘聚裂纹扩展问题中的应用

强化有限单元法将物理网格与数学网格分离开来,可以方便地描述非连续变形;粘聚区域模型是模拟断裂过程区作用最简单有效的方法,且可以避免裂纹尖端的应力奇异性.本文以平面问题为例,将强化有限单元法与粘聚区域模型相结合,利用富集数学节点描述任意粘聚裂纹扩展过程中的非连续变形问题,提出了裂纹扩展过程中数学节点富集和数学单元定义的方法.本文还导出了与平面4～8节点平面等参单元对应的8～16节点粘聚裂纹单元列武.最后,通过三点弯梁的裂纹扩展过程模拟验证了本文提出的粘聚裂纹扩展模拟方法的有效性.     

Numerical simulation of quasi-brittle fracture using damaging cohesive surfaces

The cohesive surface methodology is used in a numerical study of fracture of concrete. The traction vs. separation response is governed by an isotropic damage law in which damage evolves according to a prescribed one-dimensional linear or exponential softening law. Cohesive surfaces are immersed in the continuum to allow for a maximum freedom of crack path selection. The single edge notched four point shear beam and the double edge notched tensile bar are used to study: (i) the influence of the tangential cohesive response on the development of the fracture path and (ii) the mesh alignment sensitivity. It is shown that in the present formulation, the tangential cohesive response has a minor influence on both crack path and global characteristics. Mesh alignment does have a significant influence on the outcome of the numerical analysis.     

Process zone in the Single Cantilever Beam under transverse loading - Part II: Experimental

This paper describes an experimental arrangement to evaluate stress/strain fields in the process zone of asymmetric adhesively bonded joints. A transparent polycarbonate flexible beam was bonded to an aluminium alloy rigid block with an epoxy adhesive in a Single Cantilever Beam (SCB) configuration. The flexible adherend was loaded in the direction parallel to the initial crack front at constant rate. To monitor strains induced by bending and shear along the beam, electric strain gauges were attached to the upper surface of the flexible adherend. Thus strain distribution was measured above the bonded surface, which could be used to monitor crack propagation and investigate stress redistribution in the process zone. A Timoshenko beam lying on a Pasternak elastic foundation model was used for the analysis of experimental findings. Subsequently, the Digital Image Correlation technique was used to measure the flexible substrate in-plane displacement field in the vicinity of the crack front and to assess the specimen kinematics. We found that strain gauge instrumentation of the fracture mechanics specimen was a very sensitive technique for experimental analysis of crack propagation under complex loading, offering fine investigation of stress distribution in the cohesive zone.     

Nonlinear interface shear fracture of end notched flexure specimens

The nonlinear analytical solutions of an end notched flexure adhesive joint or fracture test specimen with identical or dissimilar adherends are investigated. In the current study, a cohesive zone model (with arbitrary nonlinear cohesive laws) based analytical solution is obtained for the interface shear fracture of an end notched flexure (ENF) specimen with sufficiently long bond length. It is found that the scatter and inconsistency in calculating Mode II toughness may be significantly reduced by this model. The present work indicates that the Mode II toughness G_IIc under pure shear cracking condition is indeed very weakly dependent on the initial crack length. And this conclusion is well supported by the experimental results found in the literature. The parametric studies show that the interface shear strength is the most dominant parameter on the critical load. It is also interesting to note that with very short initial crack length and identical interface shear strength, higher Mode II toughness indeed cannot increase the critical load. Unlike the high insensitivity of critical load to the detailed shape of the cohesive law for Mode I peel fracture, the shape of the cohesive law becomes relatively important for the critical load of joints under pure Mode II fracture conditions, especially for joints with short initial crack length. The current study may help researchers deepen the understanding of interface shear fracture and clarify some previous concepts on this fracture mode.     

Simulation of local instabilities during crack propagation in the ductile–brittle transition region

The crack propagation for a cohesive zone model within an elastic-plastic material under small-scale yielding conditions is simulated numerically. The resulting crack growth resistance curves show local instabilities, so-called pop-ins even for homogeneous cohesive properties if the cohesive strength lies sufficiently close to the maximum stress of the corresponding blunting solution. The formation of secondary cracks and unloading zones in front of the actual crack tip is identified as the underlying mechanism. It is found that the shape of the cohesive law has a considerable influence on the crack arrest behavior. Furthermore, requirements to the mesh resolution are derived.     

基于有限断裂法和比例边界有限元法的裂纹扩展模拟

基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。     

The simulation of dynamic crack propagation using the cohesive segments method

The cohesive segments method is a finite element framework that allows for the simulation of the nucleation, growth and coalescence of multiple cracks in solids. In this framework, cracks are introduced as jumps in the displacement field by employing the partition of unity property of finite element shape functions. The magnitude of these jumps are governed by cohesive constitutive relations. In this paper, the cohesive segments method is extended for the simulation of fast crack propagation in brittle solids. The performance of the method is demonstrated in several examples involving crack growth in linear elastic solids under plane stress conditions: tensile loading of a block; shear loading of a block and crack growth along and near a bi-material interface.     

A crack tip tracking algorithm for cohesive interface element analysis of fatigue delamination propagation in composite materials

A novel approach is proposed for the use of cohesive elements in the analysis of delamination propagation in composite materials under high-cycle fatigue loading. The method is applicable to delamination propagation within the Paris-law regime and is suitable for the analysis of three-dimensional structures typical of aerospace applications. The major advantages of the proposed formulation are its complete independence of the cohesive zone length – which is a geometry-dependent parameter – and its relative insensitivity to mesh refinement. This is only possible via the introduction of three nonlocal algorithms, which provide (i) automated three-dimensional tracking of delamination fronts, (ii) an estimation of direction of crack propagation and (iii) accurate and mesh-insensitive integration of strain energy release rate. All calculations are updated at every increment of an explicit time-integration finite element solution, which models the envelopes of forces and displacements with an assumption of underlying constant cyclic loading. The method was implemented as a user-defined subroutine in the commercial finite element software LS-Dyna and supports the analysis of complex three-dimensional models. Results are presented for benchmark cases such as specimens with central cut plies and centrally-loaded circular plates. Accurate predictions of delamination growth rates are observed for different mesh topologies in agreement with the Paris-laws of the material.     

压应力对压剪裂纹扩展的影响研究

本文针对压剪裂纹的启裂及扩展问题,研究了脆性材料中裂纹面压应力变化对其扩展的影响规律.借助双轴加载试验机可自由调节两个轴位移和力的优势,设计了一种单边对称双裂纹压剪试样.试验中,施加裂纹面压应力至不同的预设值后,使剪应力单调增大直至裂纹启裂及扩展,得到不同预设压应力下压剪裂纹启裂及扩展规律. 随着预设压应力增大,启裂角增大,裂纹扩展路径与初始裂纹的偏角也增大.但随着压应力增大,启裂角增大至一定值后趋于稳定,实验发现,依据裂纹是否闭合,压应力对压剪裂纹扩展的作用大致可分为两个阶段:闭合前,压应力对裂纹启裂载荷及启裂角、扩展路径均有影响,预设压应力增大,裂纹启裂载荷增大、启裂角增大,扩展路径愈来愈偏离初始裂纹方向;闭合后,压应力虽然增大,启裂角和临界压剪应力比始终恒定,压应力对临界剪力和扩展路径存在一定影响.研究表明,裂纹启裂角与启裂时的压剪应力比存在一定的对应关系.启裂时的压剪应力比增大,启裂角增大,启裂时的压剪应力比恒定,启裂角不变.      

复合砂浆钢绞线网加固RC梁受弯性能计算研究

分析了复合砂浆钢绞线网加固钢筋混凝土梁的受弯性能, 不同材料加固梁的受弯性能可采用同一公式进行计算, 研究表明加固梁受弯性能受本体梁的初始受力状态和是否采用三面U型加固影响较大. 采用正截面极限平衡法, 推导出复合砂浆钢绞线网加固钢筋混凝土梁极限抗弯承载力计算公式; 在混凝土结构设计规范公式基础上, 提出加固梁裂缝宽度和裂缝间距计算公式, 加固梁公式计算值与试验值吻合较好.     

利用无网格方法分析斜拉破坏钢筋混凝土梁

从传统的混凝土分离裂缝模型和弥散裂缝模型出发，使用不同的模型模拟不同发展阶段的混凝土裂缝，并利用无网格方法可以灵活布置节点和边界的优点，使分离裂缝模型的使用得到很大的简化。通过在宏观裂缝表面布置基于试验结果的混凝土与混凝土界面单元，正确地模拟了裂缝表面抗剪能力的变化。算例表明，使用本文方法可以准确模拟斜拉破坏混凝土梁的破坏过程，且精度和数值稳定性要高于传统的有限元弥散裂缝方法。