水力压裂过程的扩展有限元数值模拟方法

建立了基于扩展有限元法的水力压裂数值模拟方法,使水力裂缝独立于网格存在,无需预设裂缝扩展方位。在扩展有限元计算框架下,将裂缝面处理为求解域内边界,将缝内水压力转化为相关单元等效节点力;运用考虑缝内水压力作用的相互积分法来数值求解缝尖应力强度因子;采用最大能量释放率准则确定裂缝是否继续扩展及扩展方位;最终编制了计算机程序。利用该方法数值模拟了单条水力裂缝在恒定水压力作用下作非平面扩展,所得结果分别与室内试验和解析模型相对比。结果表明,数值结果与室内试验和解析解吻合较好,缝尖应力强度因子最大相对误差不高于0.45%,验证了该方法的可行性和准确性。     

利用无网格方法分析钢筋混凝土梁开裂问题

裂缝处理一直是混凝土数值分析中的一个难点,而无网格方法由于不需要单元网格,非常适用于分析断裂问题,因此可以考虑将其引入混凝土断裂分析领域。首先将混凝土中的裂缝分为微观裂缝和宏观裂缝,对于微观裂缝仍然使用传统的弥散裂缝模型;而宏观裂缝则利用无网格方法可以非常方便地调整节点分布的特点,通过增加裂缝节点和裂面边界的方法加以模拟。并给出了宏观裂缝产生,扩展的具体模拟方法。通过算例表明,利用提出的方法,可以较准确地分析混凝土宏观裂缝的产生、扩展以及裂缝宽度变化,得到一些传统有限元方法难以得到的计算结果。     

多边形比例边界有限元非线性高效分析方法

比例边界有限元法作为一种高精度的半解析数值求解方法,特别适合于求解无限域与应力奇异性等问题,多边形比例边界单元在模拟裂纹扩展过程、处理局部网格重剖分等方面相较于有限单元法具有明显优势。目前,比例边界有限元法更多关注的是线弹性问题的求解,而非线性比例边界单元的研究则处于起步阶段。该文将高效的隔离非线性有限元法用于比例边界单元的非线性分析,提出了一种高效的隔离非线性比例边界有限元法。该方法认为每个边界线单元覆盖的区域为相互独立的扇形子单元,其形函数以及应变-位移矩阵可通过半解析的弹性解获得;每个扇形区的非线性应变场通过设置非线性应变插值点来表达,引入非线性本构关系即可实现多边形比例边界单元高效非线性分析。多边形比例边界单元的刚度通过集成每个扇形子单元的刚度获取,扇形子单元的刚度可采用高斯积分方案进行求解,其精度保持不变。由于引入了较多的非线性应变插值点,舒尔补矩阵维数较大,该文采用Woodbury近似法对隔离非线性比例边界单元的控制方程进行求解。该方法对大规模非线性问题的计算具有较高的计算效率,数值算例验证了算法的正确性以及高效性,将该方法进行推广,对实际工程分析具有重要意义。     

奇异薄层单元及其在双材料界面断裂分析中的应用

构造了一种平面应力奇异薄层单元并证明了其具有-1/2阶奇异性。用此单元研究了双材料界面层的刚度对界面裂纹尖端场的影响。研究发现:对于I 型界面断裂, 减小界面法向刚度对K_Ⅰ、K_Ⅱ的影响远大于减小界面切向刚度, 且法向和切向刚度的减小对K_Ⅱ的影响均大于对K_Ⅰ的影响, 降低法向刚度会显著改变裂尖正应力和剪应力的分布, 而降低切向刚度只明显改变剪应力的分布, 对正应力的分布影响不大;对于界面II 型断裂, 则减小切向刚度对K_Ⅰ、K_Ⅱ的影响远大于减小法向刚度, 且切向、法向刚度的减小对K_Ⅰ影响均大于对K_Ⅱ影响, 降低切向刚度会显著改变裂尖正应力和剪应力的分布, 而降低法向刚度只明显改变裂尖正应力的分布, 对剪应力的分布影响不大。随着界面刚度增大, 应力强度因子和裂尖应力分布均趋近无厚度理想界面情况。     

设周边缝对小湾高拱坝工作性态的影响

用非线性有限元对小湾高拱坝(坝高292．0m)设周边缝时的工作性态进行研究。模型中考虑了缝面应力与变位的非线性重调整,考虑了坝体、坝基材料非线性和周边缝接触非线性的耦合,模拟了缝的开裂、滑移等非线性。对设缝深度、范围、以及干缝、湿缝等对高拱坝拱冠、坝踵和坝肩等关键部位的位移场和应力场的影响进行了分析研究,并对设缝前后建基面、坝踵附近的开裂区和屈服区进行分析,研究结果表明,设周边缝可以有效地降低高拱坝坝踵的粱向拉应力,对坝踵附近应力影响较大,坝体其它部位应力和位移变化不大。     

混凝土非线性断裂韧度G_(Ic)及其尺寸效应

就混凝土而言，在裂缝失稳扩展前，由于裂缝尖端的微裂区而产生非线性变形并伴有裂缝的亚临界扩展。这种非线性变形与亚临界扩展，使混凝土材料在断裂前吸收更多的能量。本文从试验得出的荷载－位移的非线性曲线出发，求得了材料抵抗裂缝扩展的非线性断裂韧度，并探讨了试件缝高比及高度对的影响。     

基于粘结单元的三维随机细观混凝土离散断裂模拟

开发MATLAB和ABAQUS Python脚本程序,建立含随机三维多面体骨料的混凝土细观有限元模型。通过自编高效C++程序插设三维零厚度的粘结界面单元,以模拟复杂三维离散多裂缝的起裂与扩展。对典型混凝土试件单轴拉伸断裂的模拟,分析粘性界面单元的主要材料参数(抗拉强度与断裂能)对应力-位移曲线、断裂过程、裂缝面特征的影响。结果表明:宏观应力-位移曲线主要受砂浆、界面粘性裂缝单元的抗拉强度、断裂能绝对数值的影响;裂缝面的位置与形态决定于砂浆、界面粘结裂缝单元的抗拉强度相对比值以及断裂能相对比值;混凝土的力学响应反映其裂缝发展特征,二者既决定于断裂材料参数,也受到骨料大小、形状等细观结构因素的影响;建立的有限元模型能有效描述混凝土复杂三维断裂过程。     

线粘弹性断裂分析的增量加料有限元法

提出了一种用于解决线粘弹性断裂问题的增量加料有限元法。为了反映裂纹尖端的应力奇异性,在裂尖附近的应力奇异区采用若干四边形加料单元和过渡单元,非奇异区采用常规四边形单元,三种单元分区混合使用形成求解域网格划分。加料单元通过引入裂尖渐近位移场,构造出可以较好反映裂尖奇异性的单元位移模式,过渡单元在加料单元基础上引入调整函数构造单元位移模式,用于连接加料单元和常规单元,以消除加料单元和常规单元间位移不协调。基于Boltzmann叠加原理,推导了粘弹性材料的增量型本构关系,进而获得了增量加料有限元列式,并基于节点位移外推法计算粘弹性介质中裂纹应变能释放率。数值算例验证了该文方法的正确性和有效性。     

基于X-SBFEM的裂纹体非网格重剖分耦合模型研究

扩展有限元法利用了非网格重剖分技术,但需要基于裂尖解析解构造复杂的插值基函数,计算精度受网格疏密和插值基函数等因素影响。比例边界有限元法则在求解无限域和裂尖奇异性问题优势明显,两者衔接于有限元法理论内,可建立一种结合二者优势的断裂耦合数值模型。该文从虚功原理出发,利用位移协调与力平衡机制,提出了一种断裂计算的新方法X-SBFEM,达到了扩展有限元模拟裂纹主体、比例边界有限元模拟裂尖的目的。在数值算例中,通过边裂纹和混合型裂纹的应力强度因子计算,并与理论解对比,验证了该方法的准确性和有效性。     

基于水平集算法的扩展有限元方法研究

扩展有限元是一种以单位分解思想为基础,在常规有限元位移中加入跳跃函数和渐近位移场函数,以处理不连续问题的数值方法。将水平集算法应用到裂纹界面的描述及加强单元类型的判别,并与扩展有限元相结合,用于分析材料断裂问题。相比传统有限元,有限元网格与裂纹面位置相互独立,不需满足裂纹为单元边、裂尖为单元节点和在裂纹附近进行高密度的...     

A novel approach based on ALE and delamination fracture mechanics for multilayered composite beams

A novel approach able to predict debonding or fracture phenomena in multilayered composite beams is proposed. The structural model is based on the first-order shear deformable laminated beam theory and moving mesh strategy developed in the framework of Arbitrary Lagrangian–Eulerian (ALE) formulation. The former is utilized to evaluate fracture parameters by using a multilayer approach, in which a low number of interface elements are introduced along the thickness, whereas the latter is utilized to reproduce crack tip motion due to the crack extension produced by moving boundaries. The model is able to avoid computational complexities introduced by an explicit crack representation in bi-dimensional structures, in which typically high computational efforts are expected for handling moving boundaries. To this aim, a moving mesh strategy is proposed for the first time in the context of beam modeling based on a multilayered configuration. Such an approach, essentially based on ALE formulation, is able to reproduce interfacial crack paths by using a low number of computational elements. The numerical method is proposed in the framework of the finite element formulation for a quasi-static or dynamic evolution of the crack tip front. In order to investigate the accuracy and to validate the proposed methodology, comparisons with experimental data and existing formulations available from the literature are developed. Moreover, a parametric study in the framework of dynamic fracture is developed to investigate the capability of the proposed model to reproduce more complex loading cases.     

Non-orthogonal stress modes for interfacial fracture based on local plastic dissipation

Interfacial cracks have several features which are different from those of cracks in homogeneous materials. Among those, the loading mode dependency of interfacial toughness has been a main obstacle to the widespread utilization of interfacial fracture mechanics. In this study, plasticity-induced toughening of an interface crack between an elastic–plastic material and an elastic material is studied. A useful relationship between the plastic dissipation and the plastic zone size is derived via an effective crack length model. Non-orthogonal stress modes for interface cracks are proposed on the basis of the plastic dissipation mechanism and a mixed-mode criterion for interfacial crack growth is also proposed using these stress modes. The non-orthogonal stress modes are able to represent the asymmetric behavior, mode-dependent toughening and ε-dependency of interfacial crack growth.     

Non-orthogonal stress modes for interfacial fracture based on local plastic dissipation

Interfacial cracks have several features which are different from those of cracks in homogeneous materials. Among those, the loading mode dependency of interfacial toughness has been a main obstacle to the widespread utilization of interfacial fracture mechanics. In this study, plasticity-induced toughening of an interface crack between an elastic-plastic material and an elastic material is studied. A useful relationship between the plastic dissipation and the plastic zone size is derived via an effective crack length model. Non-orthogonal stress modes for interface cracks are proposed on the basis of the plastic dissipation mechanism and a mixed-mode criterion for interfacial crack growth is also proposed using these stress modes. The non-orthogonal stress modes are able to represent the asymmetric behavior, mode-dependent toughening and ε-dependency of interfacial crack growth.     

A study on the influence of cohesive zone interface element strength parameters on mixed mode behaviour

This paper presents a detailed study of the influence of maximum interfacial stress on interface element analyses for composites delamination. The development of the non-linear cohesive zone ahead of a crack tip is analysed with respect to length, stress distribution and mode ratio. The energy absorbed by interface elements is compared with the crack tip strain energy release rate from fracture mechanics analyses. These studies are performed initially on standard fracture toughness specimens, where mode-ratio is fixed by the applied displacement constraints. Results show close agreement with linear elastic fracture mechanics solutions. A simple ply drop specimen is then modelled, where the mode ratio is not constrained by the boundary conditions, and results are compared with the Virtual Crack Closure Technique. In this case maximum interfacial stress has a far greater influence on the numerical results, due to its significant influence on cohesive zone length, mode ratio and energy absorbed.     

纤维增强复合材料界面脱层和基体裂纹的模拟分析

基于Ghosh提出的Voronoi单元有限元方法,构造能同时反映纤维增强复合材料界面脱层和基体裂纹扩展的单元(X-VCFEM单元);应用界面力学理论和断裂力学理论,建立界面脱层、界面裂纹扩展方向和基体裂纹扩展的判断准则;结合网格重划分技术,模拟分析了只含有一个夹杂时界面脱层和基体裂纹扩展的过程,并通过与传统有限元计算结果的比较,验证X-VCFEM单元的可靠性和有效性;同时,模拟分析含任意随机分布夹杂的纤维增强复合材料界面脱层和基体裂纹的产生和扩展过程。结果表明：应用该方法模拟复杂多相复合材料裂纹问题具有计算速度快和精度高的优越性。     

Experimental and numerical investigations on fracture process zone of rock–concrete interface

A crack propagation criterion for a rock–concrete interface is employed to investigate the evolution of the fracture process zone (FPZ) in rock–concrete composite beams under three‐point bending (TPB). According to the criterion, cracking initiates along the interface when the difference between the mode I stress intensity factor at the crack tip caused by external loading and the one caused by the cohesive stress acting on the fictitious crack surfaces reaches the initial fracture toughness of a rock–concrete interface. From the experimental results of the composite beams with various initial crack lengths but equal depths under TPB, the interface fracture parameters are determined. In addition, the FPZ evolution in a TPB specimen is investigated by using a digital image correlation technique. Thus, the fracture processes of the rock–concrete composite beams can be simulated by introducing the initial fracture criterion to determine the crack propagation. By comparing the load versus crack mouth opening displacement curves and FPZ evolution, the numerical and experimental results show a reasonable agreement, which verifies the numerical method developed in this study for analysing the crack propagation along the rock–concrete interface. Finally, based on the numerical results, the effect of ligament length on the FPZ evolution and the variations of the fracture model during crack propagation are discussed for the rock–concrete interface fracture under TPB. The results indicate that ligament length significantly affects the FPZ evolution at the rock–concrete interface under TPB and the stress intensity factor ratio of modes II to I is influenced by the specimen size during the propagation of the interfacial crack.     

Numerical analysis of the stress–strain distributions in the particle reinforced metal matrix composite SiC/6064Al

The axisymmetric cell model consisting of interface, matrix and reinforced particle is used to simulate the tensile test of particle reinforced metal matrix composite for predicting the micro stress/strain field and macro tensile stress/strain curve. In simulation of the tensile test, the cohesive element model is selected to model interfacial crack growth. It mainly analyzed the effects of interfacial properties, reinforcement volume fractions and aspect ratios on the stress–strain states of particle reinforced metal matrix composite. The results show that the peak micro stress and plastic strain occur at the interface in which it is a certain angle from the tensile stress direction; with the interfacial fracture toughness and reinforcement volume fraction increasing, the flow stress increases firstly and then decreases. The tensile stress–strain properties of SiC/6064Al are good when the interfacial fracture toughness is equal to 60 J/m and the reinforcement fraction volume is equal to 20%. Smaller reinforcement aspect ratio leads to smaller micro stress in composites.     

Delamination analysis of composites by new orthotropic bimaterial extended finite element method

The extended finite element method (XFEM) is further improved for fracture analysis of composite laminates containing interlaminar delaminations. New set of bimaterial orthotropic enrichment functions are developed and utilized in XFEM analysis of linear‐elastic fracture mechanics of layered composites. Interlaminar crack‐tip enrichment functions are derived from analytical asymptotic displacement fields around a traction‐free interfacial crack. Also, heaviside and weak discontinuity enrichment functions are utilized in modeling discontinuous fields across interface cracks and bimaterial weak discontinuities, respectively. In this procedure, elements containing a crack‐tip or strong/weak discontinuities are not required to conform to those geometries. In addition, the same mesh can be used to analyze different interlaminar cracks or delamination propagation. The domain interaction integral approach is also adopted in order to numerically evaluate the mixed‐mode stress intensity factors. A number of benchmark tests are simulated to assess the performance of the proposed approach and the results are compared with available reference results. Copyright © 2011 John Wiley & Sons, Ltd.     

Stress intensity factors of a surface crack near an interface end

Due to the singular behavior of the stress field near the interface edge of bonded dissimilar materials, fracture generally initiates near the interface edge, or just from the interface edge point. In this paper, an edge crack near the interface, which can be considered as being induced by the edge singularity and satisfying two conditions, is analyzed theoretically, based on the singular stress field near the interface edge and the superposition principle. It is found that the stress intensity factor can be expressed by the stress intensity coefficient of the edge singular stress field, the crack length, the distance between the interface and the crack, as well as the material combination. Boundary element method analysis is also carried out. It is found that the theoretical result coincides well with the numerical result when the crack length is small. Therefore, the theoretical representation obtained by this study can be used to simply evaluate the stress intensity factor of an edge singularity induced crack for this case. However, when the crack length becomes larger than a certain value, a significant difference appears, especially for the case with large edge singularity.     

FRP片材与混凝土粘结性能的精细有限元分析

FRP-混凝土界面粘结性能是外贴FRP片材加固混凝土结构技术的关键问题。基于FRP与混凝土界面面内剪切试验,采用精细单元有限元模型对其界面粘结性能进行了研究。在该模型中,混凝土和FRP片材都使用非常小的单元加以模拟,通过调整混凝土材料的本构模型来考虑单元尺寸的影响。FRP单元和混凝土单元直接连接,通过混凝土单元的断裂破坏来模拟FRP和混凝土界面的宏观剥离破坏过程。通过与大量面内剪切试验结果对比,验证了该精细有限元模型的正确性,并基于精细有限元分析结果,对界面剥离破坏机理进行了讨论。