内聚力模型裂纹问题分析的解析奇异单元

内聚力模型已经被广泛应用于需要考虑断裂过程区的裂纹问题当中,然而常用的数值方法应用于分析内聚力模型裂纹问题时还存在着一些不足,比如不能准确的给出断裂过程区的长度、需要网格加密等。为了克服这些缺点,论文构造了一个新型的解析奇异单元,并将之应用于基于内聚力模型的裂纹分析当中。首先将虚拟裂纹表面处的内聚力用拉格拉日插值的方法近似表示为多项式的形式,而多项式表示的内聚力所对应的特解可以被解析地给出。然后利用一个简单的迭代分析,基于内聚力模型的裂纹问题就可以被模拟出来了。最后,给出二个数值算例来证明本文方法的有效性。     

数学网格和物理网格分离的有限单元法(I):基本理论

常规有限单元法在复杂边界问题的网格剖分、可移动边界和非连续变形问题的数值模拟等方面存在困难.本文将常规的有限单元分离为几何上相互独立的数学单元和物理单元,基于数学单元构造近似函数,引入位移模式关联法则以确定物理单元的位移模式,提出了在现有有限单元法框架内、基于数学网格和物理网格分离的强化有限单元法(FEM++).与常规有限单元法(SFEM)比较表明,强化有限单元法不仅很好地克服了常规有限单元法网格剖分上的困难,而且提供了一条更简便、更自然的分析移动边界问题和非连续变形问题的新途径.最后,通过数值算例验证了强化有限单元法的适用性和有效性.     

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

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

混凝土断裂过程区的虚拟裂纹粘聚力奇异性

混凝土断裂过程区视为具有粘聚阻力作用的虚拟裂纹，其非线性断裂和尺寸效应特性是与该虚拟裂纹粘聚力分布规律密切相关的。通过得到的粘聚应力分布函数解析结果，对该粘聚力分布特征的分析得知，在基于断裂过程区之外用线弹性场的力学模型上，该粘聚力随距离虚拟裂纹尖点的靠近，仍具有平方根奇异性。从而本文提出一个能够反映裂纹发展状态的粘聚应力奇异性强度参数，它是无粘聚力的线弹性裂纹应力强度因子和表征裂纹张开位移分布多项式参数的函数；因此，该参数可以作为混凝土非线性断裂的一个参量。文中就已有断裂试验测试结果进行了算例分析和相应的讨论。     

Cohesive zone model of intermediate crack-induced debonding of FRP-plated concrete beam

External bonding of FRP plates or sheets has emerged as a popular method for strengthening reinforced concrete structures. Debonding along the FPR–concrete interface can lead to premature failure of the structures. In this study, debonding induced by a flexural crack in a FRP-plated concrete beam is analyzed through a nonlinear fracture mechanics method. The concrete beam and FRP plate are modeled as linearly elastic simple beams connected together through a thin layer of FRP–concrete interface. A bi-linear cohesive (bond-slip) law, which has been verified by experiments, is used to model the FRP–concrete interface as a cohesive zone. Thus a cohesive zone model for intermediate crack-induced debonding is established with a unique feature of unifying the debonding initiation and growth into one model. Closed-form solutions of interfacial stress, FRP stress and ultimate load of the plated beam are obtained and then verified with the numerical solutions based on finite element analysis. Parametric studies are carried out to demonstrate the significant effect of FRP thickness on the interface debonding. The bond-slip shape is examined specifically. In spite of its profound effect on softening zone size, the bond-slip shape has been found to have little effect on the ultimate load of the plated beam. By making use of such a unique feature, a simplified explicit expression is obtained to determine the ultimate load of the plated concrete beam with a flexural crack conveniently. The cohesive zone model in this study also provides an efficient and effective way to analyze more general FRP–concrete interface debonding.     

Analysis on the cohesive stress at half infinite crack tip

The nonlinear fracture behavior of quasi-brittle materials is closely related with the cohesive force distribution of fracture process zone at crack tip. Based on fracture character of quasi-brittle materials, a mechanical analysis model of half infinite crack with cohesive stress is presented. A pair of integral equations is established according to the superposition principle of crack opening displacement in solids, and the fictitious adhesive stress is unknown function . The properties of integral equations are analyzed, and the series function expression of cohesive stress is certified. By means of the data of actual crack opening displacement, two approaches to gain the cohesive stress distribution are proposed through resolving algebra equation. They are the integral transformation method for continuous displacement of actual crack opening, and the least square method for the discrete data of crack opening displacement. The calculation examples of two approaches and associated discussions are give     

A new quasi-continuum constitutive model for crack growth in an isotropic solid

A new quasi-continuum constitutive model is established based on the randomized cohesive bonds model proposed by Gao and Klein (1998). This model bridges the microscopic discrete constitution characters and the macroscopic mechanical properties of material. In the presented constitutive model, both the bond stretch energy potential and the rotation energy potential are considered, which makes the presented constitutive model applicable to different Poisson-ratio and Young's modulus materials. By establishing a phenomenological bond stiffness function according to the complete stress–strain relationship of uniaxial tension test, the fracture criterion is directly incorporated into the constitutive model. The method requires no external fracture criterion when simulating fracture initiation and propagation, which brings convenience in the numerical simulation. At last, the presented constitutive model is applied to an example of crack growth in an isotropic solid.     

A thermo-mechanical cohesive zone formulation for ductile fracture

The paper addresses the possibility to project both mechanical and thermal phenomena pertinent to the fracture process zone into a cohesive zone. A wider interpretation of the notion cohesive zone is thereby suggested to comprise not only stress degradation due to micro-cracking but also heat generation and energy transport. According to our experience, this widening of the cohesive zone concept allows for a more efficient finite element simulation of ductile fracture. The key feature of the formulation concerns the thermo-mechanical cohesive zone model, evolving within the thermo-hyperelastoplastic continuum, allowing for the concurrent modelling of both heat generation, due to the fracture process, and heat transfer across the fracture process zone. This is accomplished via thermodynamic arguments to obtain the coupled governing equation of motion, energy equation, and constitutive equations. The deformation map is thereby defined in terms of independent continuous and discontinuous portions of the displacement field. In addition, as an extension of the displacement kinematics, to represent the temperature field associated with the discontinuous heat flux across the fracture interface, a matching discontinuous temperature field involving the interface (or band) temperature is proposed. In the first numerical example, concerning dynamic quasi-brittle crack propagation in a thermo-hyperelastoplastic material, we capture the initial increase in temperature close to the crack surface due to the energy dissipating fracture process. In the second example, a novel application of ductile fracture simulation to the process of high velocity (adiabatic) cutting is considered, where some general trends are observed when varying the cutting velocity.     

Size effect in quasibrittle failure: Analytical model and numerical simulations using cohesive zone model

The recent rewriting of the Ba?ant’s size effect law (Morel, 2008) which has suggested the existence of an additional asymptotic regime for intermediate structure sizes is now compared to numerical simulations of fracture of geometrically similar notched structures of different sizes extending over 2.4 decades. The quasibrittle fracture behavior is simulated through cohesive zone model (bilinear softening) using a constant set of cohesive parameters whatever the specimen size D is. The R-curves resulting from the load–displacement responses are estimated and appear as size-independent. On this basis, the different asymptotic regimes expected for the size effect on fracture properties at peak load such as the relative crack length, the resistance to crack growth and the nominal strength are shown in fair agreement with the size effect observed on the results obtained from numerical simulations.     

Analytical representation of the non-square-root singular stress field at a finite angle sharp notch

The stress field near the tip of a finite angle sharp notch is singular. However, unlike a crack, the order of the singularity at the notch tip is less than one-half. Under tensile loading, such a singularity is characterized by a generalized stress intensity factor which is analogous to the mode I stress intensity factor used in fracture mechanics, but which has order less than one-half. By using a cohesive zone model for a notional crack emanating from the notch tip, we relate the critical value of the generalized stress intensity factor to the fracture toughness. The results show that this relation depends not only on the notch angle, but also on the maximum stress of the cohesive zone model. As expected the dependence on that maximum stress vanishes as the notch angle approaches zero. The results of this analysis compare very well with a numerical (finite element) analysis in the literature. For mixed-mode loading the limits of applicability of using a mode I failure criterion are explored.     

Introducing new cracked finite elements and a method for SIF calculation of cracks

Two different types of 8-node cracked quadrilateral finite element are presented for fracture applications. The first element contains a central crack and the other one includes an edge crack. The introduced elements are applicable in 2D problems. The crack is not physically modeled within the element, but instead, its effects on the stiffness matrix are taken into account by utilizing linear fracture mechanics laws. Furthermore, a simple and practical procedure is proposed for calculation of stress intensity factor (SIF) by employing proposed cracked elements. Several numerical examples are presented to evaluate the capabilities of the proposed elements and procedure.     

Pure mode II fracture characterization of composite bonded joints

A new data reduction scheme is proposed for measuring the critical fracture energy of adhesive joints under pure mode II loading using the End Notched Flexure test. The method is based on the crack equivalent concept and does not require crack length monitoring during propagation, which is very difficult to perform accurately in these tests. The proposed methodology also accounts for the energy dissipated at the Fracture Process Zone which is not negligible when ductile adhesives are used. Experimental tests and numerical analyses using a trapezoidal cohesive mixed-mode damage model demonstrated the good performance of the new method, namely when compared to classical data reduction schemes. An inverse method was used to determine the cohesive properties, fitting the numerical and experimental load–displacement curves. Excellent agreement between the numerical and experimental R-curves was achieved demonstrating the effectiveness of the proposed method.     

On the determination of the cohesive zone properties of an adhesive layer from the analysis of the wedge-peel test

An extensive numerical study of the mechanics of the “wedge-peel test” is performed in order to analyze the mode I steady state debonding of a sandwich structure made of two thin plastically deforming metallic plates bonded with an adhesive. The constitutive response of the metallic plates is modeled by J₂ flow theory, and the behavior of the adhesive layer is represented with a cohesive zone model characterized by a maximum separation stress and the fracture energy. A steady-state finite element code accounting for finite rotation has been developed for the analysis of this problem. Calculations performed with the steady-state formulation are shown to be much faster than simulations involving both crack initiation and propagation within a standard, non-steady-state code. The goal of this study is to relate the measurable parameters of the test to the corresponding fracture process zone characteristics for a representative range of adherent properties and test conditions. An improved beam bending model for the energy release rate is assessed by comparison with the numerical results. Two procedures are proposed for identifying the cohesive zone parameters from experimental measurements.     

基于图像四叉树的改进型比例边界有限元法研究

为提高数值计算的精度,断裂力学问题的数值模拟需要在裂纹扩展的局部区域采用较密的网格,而远离裂纹扩展的区域可采用较疏的网格,且对于裂纹扩展问题的数值模拟,大多数数值方法又存在局部网格重剖分的问题.论文提出了一种基于图像四叉树的改进型比例边界有限元法用于模拟裂纹扩展问题,该方法可根据结构域几何外边界的图像全自动进行四叉树网格剖分,无需任何人工干预,网格剖分效率极高,由于比例边界有限元法本身的优势,四叉树网格的悬挂节点可以直接地视为新的节点,无需任何特殊处理.通过引入虚节点的思想,将裂纹与四叉树单元边界交叉点作为虚节点,虚节点的自由度作为附加自由度处理,并采用水平集函数表征材料内部的裂纹面,含不连续裂纹面的子域可通过节点水平集函数识别,使得裂纹扩展时无需进行网格重剖分,界面的几何特征通过比例边界有限元子域的附加自由度表征.最后,通过若干算例验证了该方法的性能,建议的改进型比例边界有限元法在求解复合型应力强度因子和模拟材料内部裂纹扩展路径时均具有较高的精度.     

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

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

An isolated cohesive crack in tension

The shape of the back-calculated stress-separation curves obtained from the in-situ fracture of first-year sea ice in the Arctic and Antarctic is convex, and radically different from those for concrete. In these tests, the process zone size changes with crack growth, but the nature of this change differs with the test conditions, load-control or displacement control. These results prompted a closer examination of the cohesive crack model using the simplest cracked configuration, a finite cohesive crack in an infinite elastic medium loaded in tension by a uniform stress at infinity. Different types of strain softening are examined: rectangular softening, linear softening, prescribed cohesive stresses, and prescribed cohesive crack-opening displacements. For each of these cases, crack nucleation is examined; close attention is paid to test control conditions, be they load-control or fixed-grip. The test control conditions alter the fracture, crack nucleation, crack growth, and process zone size behavior significantly. Accurate approximate solutions to linear softening are presented and examined.     

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.     

准脆性材料中强不连续问题的数值方法若干应用及其研究进展

综述了模拟准脆性材料开裂过程的数值计算方法的研究进展和工程应用,比较了表征强不连续问题的显式非连续模型和隐式非连续模型的优缺点.结合混凝土粘结裂纹, 重点讨论了嵌入非连续模型,扩展有限元方法和富集有限元技术等非连续方法的构造特征和本质区别.从各种富集方法的理论完备性考察,以假定发展应变为基础的嵌入非连续方法虽然可以解决混凝土开裂过程中的应力锁死,满足内部边界的静力平衡条件以及反映开裂后的位移不连续问题,但嵌入非连续所采用的富集函数在开裂单元中并不能满足协调条件,使非连续两侧的应变不独立. 其局限性是由于富集自由度在单元的水平上引入,而以单位分解为基础的扩展有限元和富集有限元的富集函数以节点自由度的方式引入,除具有嵌入非连续的优点, 还可以有效消除嵌入非连续引起裂纹两侧应变的相互影响.文中同时指出了网格重构技术,弥散裂纹模型的局限性以及扩展有限元和富集有限元技术在构造方式上的细微差别.对于节点自由度方式引入的富集函数, 其操作困难性在文中也作了说明.