FRP-混凝土三点受弯梁损伤粘结模型有限元分析

该文采用双线形损伤粘结模型研究带切口FRP-混凝土三点受弯梁(3PBB)I型加载下的界面断裂性能。通过有限元参数分析,详细讨论了界面粘结强度、界面粘结能、混凝土抗拉强度、混凝土断裂能对3PBB受力性能的影响。数值模拟表明,FRP-混凝土界面有两种破坏形式,包括FRP-混凝土界面的损伤脱粘和界面混凝土的损伤脱粘破坏,与实验所观察到的现象一致。两种破坏形式尽管在宏观上均表现为界面脱粘,但破坏机制却不同。FRP-混凝土界面的损伤粘结模型与混凝土的拉伸塑性损伤模型相结合,不但再现了3PBB的宏观力学性能,数值分析得到的荷载-位移曲线接近实验结果,而且还能详细展示FRP-混凝土界面的损伤、断裂破坏过程以及损伤在FRP-混凝土界面和界面混凝土之间的转移,能够预测构件的承载力,有助于界面优化设计,这是单纯以能量判据预测裂纹发展的经典断裂力学方法所无法做到的。     

An embedded cohesive crack model for finite element analysis of mixed mode fracture of concrete*

An embedded cohesive crack model is proposed for the analysis of the mixed mode fracture of concrete in the framework of the Finite Element Method. Different models, based on the strong discontinuity approach, have been proposed in the last decade to simulate the fracture of concrete and other quasi‐brittle materials. This paper presents a simple embedded crack model based on the cohesive crack approach. The predominant local mode I crack growth of the cohesive materials is utilized and the cohesive softening curve (stress vs. crack opening) is implemented by means of a central force traction vector. The model only requires the elastic constants and the mode I softening curve. The need for a tracking algorithm is avoided using a consistent procedure for the selection of the separated nodes. Numerical simulations of well‐known experiments are presented to show the ability of the proposed model to simulate the mixed mode fracture of concrete.     

Adaptive mesh refinement and coarsening for cohesive zone modeling of dynamic fracture

Adaptive mesh refinement and coarsening schemes are proposed for efficient computational simulation of dynamic cohesive fracture. The adaptive mesh refinement consists of a sequence of edge‐split operators, whereas the adaptive mesh coarsening is based on a sequence of vertex‐removal (or edge‐collapse) operators. Nodal perturbation and edge‐swap operators are also employed around the crack tip region to improve crack geometry representation, and cohesive surface elements are adaptively inserted whenever and wherever they are needed by means of an extrinsic cohesive zone model approach. Such adaptive mesh modification events are maintained in conjunction with a topological data structure (TopS). The so‐called PPR potential‐based cohesive model (J. Mech. Phys. Solids 2009; 57 :891–908) is utilized for the constitutive relationship of the cohesive zone model. The examples investigated include mode I fracture, mixed‐mode fracture and crack branching problems. The computational results using mesh adaptivity (refinement and coarsening) are consistent with the results using uniform mesh refinement. The present approach significantly reduces computational cost while exhibiting a multiscale effect that captures both global macro‐crack and local micro‐cracks. Copyright © 2012 John Wiley & Sons, Ltd.     

A damage-mechanics-based approach for modelling decohesion in adhesively bonded assemblies

A cohesive interface model formulated within the framework of damage mechanics is presented for the simulation of decohesion in adhesively bonded assemblies. Characteristic features of the model are: the introduction of a single energy-based damage variable for describing the damage state of the interface; use of a decohesion propagation condition relying upon the linear elastic fracture mechanics (LEFM) energy balance; a treatment for the mixed-mode situation based on the definition of an equivalent energy release rate whose expression is consistently derived from the formulation. The comparisons between numerical and experimental responses obtained for typical test problems illustrate the capabilities of the proposed approach.     

Potential-based constitutive models for cohesive interfaces: Theory,implementation and examples

Cohesive interfaces appear in many materials or structures, e.g. composites or adhesive bonds. Originally introduced to model crack tips in fracture mechanics, cohesive zone models are used to describe the constitutive behavior of cohesive interfaces since the early 1990s. In the present contribution, the concept of generalized standard materials (GSM) is transferred from the modeling of bulk behavior to cohesive zones. The potential-based framework obtained hereby is referred to as standard dissipative cohesive zones (SD-CZ). Within this framework, an irreversible interface model is developed for the one-dimensional case and subsequently extended to three-dimensional transverse isotropy. While the potential structure of the constitutive law may be required for certain applications, it also brings benefits with regard to the numerical implementation. To the best knowledge of the authors, this is the first approach to interface modeling in a GSM-like framework, where dissipative effects like damage and softening are considered as well as normal-tangential coupling for mixed-mode decohesion. Comparisons to experimental data and existing cohesive zone models outline the capabilities of the approach.     

Hybrid equilibrium element with interelement interface for the analysis of delamination and crack propagation problems

This article proposes a formulation for the analysis of delamination and fracture propagation problems at the interelement interface, with perfect adhesion at the pre-failure condition and with linear softening at the post-failure regime. The proposed formulation is based on the hybrid equilibrium element (HEE) model, with stress fields which strongly verify the homogeneous equilibrium equations and interelement equilibrium equations. The HEE can easily model high-order stress fields and can implicitly model the initially rigid behavior of an extrinsic interface at the element sides. The interface model is defined as a function of the same degrees of freedom of the HEE (generalized stresses) and the pre- and post-failure behavior of the interface can be modeled without any additional degree of freedom. The proposed formulation is developed for a nine-node triangular HEE with quadratic stress fields. This article also proposes an extrinsic cohesive model, rigorously developed in the damage mechanics framework.     

Cracks and Plasticity

The paper gives a review of the state‐of‐the‐art of modeling deformation, damage, and fracture of ductile metals on various length scales, namely numerical simulation of micro‐void nucleation and growth, porosity‐based constitutive equations, the cohesive zone model and the macroscopic fracture‐mechanics approach. Links between the various models are outlined and their respective favorable and preferential applications are discussed. The potential of a hierarchical modeling approach from the micro‐ to the macroscale is emphasized.     

An adaptive spacetime discontinuous Galerkin method for cohesive models of elastodynamic fracture

This paper describes an adaptive numerical framework for cohesive fracture models based on a spacetime discontinuous Galerkin (SDG) method for elastodynamics with elementwise momentum balance. Discontinuous basis functions and jump conditions written with respect to target traction values simplify the implementation of cohesive traction–separation laws in the SDG framework; no special cohesive elements or other algorithmic devices are required. We use unstructured spacetime grids in a h‐adaptive implementation to adjust simultaneously the spatial and temporal resolutions. Two independent error indicators drive the adaptive refinement. One is a dissipation‐based indicator that controls the accuracy of the solution in the bulk material; the second ensures the accuracy of the discrete rendering of the cohesive law. Applications of the SDG cohesive model to elastodynamic fracture demonstrate the effectiveness of the proposed method and reveal a new solution feature: an unexpected quasi‐singular structure in the velocity response. Numerical examples demonstrate the use of adaptive analysis methods in resolving this structure, as well as its importance in reliable predictions of fracture kinetics. Copyright © 2009 John Wiley & Sons, Ltd.     

An extended finite element method with analytical enrichment for cohesive crack modeling

A recent approach to fracture modeling has combined the extended finite element method (XFEM) with cohesive zone models. Most studies have used simplified enrichment functions to represent the strong discontinuity but have lacked an analytical basis to represent the displacement gradients in the vicinity of the cohesive crack. In this study enrichment functions based upon an existing analytical investigation of the cohesive crack problem are proposed. These functions have the potential of representing displacement gradients in the vicinity of the cohesive crack and allow the crack to incrementally advance across each element. Key aspects of the corresponding numerical formulation and enrichment functions are discussed. A parameter study for a simple mode I model problem is presented to evaluate if quasi‐static crack propagation can be accurately followed with the proposed formulation. The effects of mesh refinement and mesh orientation are considered. Propagation of the cohesive zone tip and crack tip, time variation of the cohesive zone length, and crack profiles are examined. The analysis results indicate that the analytically based enrichment functions can accurately track the cohesive crack propagation of a mode I crack independent of mesh orientation. A mixed mode example further demonstrates the potential of the formulation. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical techniques for the analysis of crack propagation in cohesive materials

The aim of the paper is the development, assessment and use of suitable numerical procedures for the analysis of the crack evolution in cohesive materials. In particular, homogeneous as well as heterogeneous materials, obtained by embedding short stiff fibres in a cohesive matrix, are considered. Two‐dimensional Mode I fracture problems are investigated. The cohesive constitutive law is adopted to model the process zone occurring at the crack tip. An elasto‐plastic constitutive relationship, able to take into account the processes of fibre debonding and pull‐out, is introduced to model the mechanical response of the short fibres. Two numerical procedures, based on the stress and on the energy approach, are developed to investigate the crack propagation in cohesive as well as fibre‐reinforced materials, characterized by a periodic crack distribution. The results obtained using the stress and energy approaches are compared in order to evaluate the effectiveness of the procedures. Investigations on the size effect for microcracked periodic cohesive materials, and on the beneficial effects of the fibres in improving the composite material response, are developed. Copyright © 2003 John Wiley & Sons, Ltd.     

A material model to reproduce mixed‐mode fracture in concrete

This paper presents a material model to reproduce crack propagation in cement‐based material specimens under mixed‐mode loading. Its numerical formulation is based on the cohesive crack model, proposed by Hillerborg, and extended for the mixed‐mode case. This model is inspired by former works by Gálvez et al but implemented for its use in a finite element code at a material level, that is to say, at an integration point level. Among its main features, the model is able to predict the crack orientation and can reproduce the fracture behaviour under mixed‐mode fracture loading. In addition, several experimental results found in the literature are properly reproduced by the model.     

A New Procedure for Mode I Fracture Characterization of Cement‐Based Materials

Fracture characterization under mode I loading of a cement‐based material using the single‐edge‐notched beam loaded in tree‐point‐bending was performed. A new method based on beam theory and crack equivalent concept is proposed to evaluate the Resistance‐curve, which is essential to determine fracture toughness with accuracy. The method considers the existence of a stress relief region in the vicinity of the crack, dispensing crack length monitoring during experiments. A numerical validation was performed by finite element analysis considering a bilinear cohesive damage model. Experimental tests were performed in order to validate the numerical procedure. Digital image correlation technique was used to measure the specimen displacement with accuracy and without interference. Excellent agreement between numerical and experimental load–displacement curves was obtained, which validates the procedure.     

Theory and numerics for finite deformation fracture modelling using strong discontinuities

A general finite element approach for the modelling of fracture is presented for the geometrically non‐linear case. The kinematical representation is based on a strong discontinuity formulation in line with the concept of partition of unity for finite elements. Thus, the deformation map is defined in terms of one continuous and one discontinuous portion, considered as mutually independent, giving rise to a weak formulation of the equilibrium consisting of two coupled equations. In addition, two different fracture criteria are considered. Firstly, a principle stress criterion in terms of the material Mandel stress in conjunction with a material cohesive zone law, relating the cohesive Mandel traction to a material displacement ‘jump’ associated with the direct discontinuity. Secondly, a criterion of Griffith type is formulated in terms of the material‐crack‐driving force (MCDF) with the crack propagation direction determined by the direction of the force, corresponding to the direction of maximum energy release. Apart from the material modelling, the numerical treatment and aspects of computational implementation of the proposed approach is also thoroughly discussed and the paper is concluded with a few numerical examples illustrating the capabilities of the proposed approach and the connection between the two fracture criteria. Copyright © 2005 John Wiley & Sons, Ltd.     

Interface fracture in adhesively bonded shell structures

Two methods for the prediction of crack propagation through the interface of adhesively bonded shells are discussed. One is based on a fracture mechanics approach; the other is based on a cohesive zone approach. Attention is focussed on predicting the shape of the crack front and the critical stress required to propagate the crack under quasi-static conditions. The fracture mechanical model is theoretically sound and it is accurate and numerically stable. The cohesive zone model has some advantages over the fracture mechanics based model. It is easier to generalise the cohesive zone model to take into account effects such as plastic deformation in the adhering shells, and to take into account effects of large local curvatures of the interface crack front. The comparison shows a convergence of the results based on the cohesive zone model towards the results based on a fracture mechanics approach in the limit where the size of the cohesive zone becomes smaller than other relevant geometrical lengths for the problem. However, convergence issues and numerical stability must be addressed.     

Cohesive-zone modelling of the deformation and fracture of spot-welded joints

The deformation and failure of spot‐welded joints have been successfully modelled using a cohesive‐zone model for fracture. This has been accomplished by implementing a user‐defined, three‐dimensional, cohesive‐zone element within a commercial finite‐element package. The model requires two material parameters for each mode of deformation. Results show that the material parameters from this type of approach are transferable for identical spot welds in different geometries where a single parameter (such as maximum stress) is not. The approach has been demonstrated using a model system consisting of spot‐welded joints made from 5754 aluminium sheets. The techniques for determining the cohesive fracture parameters for both nugget fracture and nugget pullout are described in this paper. It has been demonstrated that once the appropriate cohesive parameters for a weld are determined, quantitative predictions can be developed for the strengths, deformations and failure mechanisms of different geometries with nominally identical welds.     

Modelling complex cracks with finite elements: a kinematically enriched constitutive model

A continuum constitutive framework with embedded cohesive interface model is presented to describe the failure of quasi-brittle materials. Both cohesive behaviour for cracking inside the fracture process zone and elastic bulk behaviour are treated at integration points making implementation straightforward. In this sense, the proposed approach is simpler than existing ones that focus on element enrichments, such as the extended finite element method, while share similarities with smeared crack models, and offers the capability to correctly model quasi-brittle failure in post-peak regime at constitutive level. In this work, the formulation is introduced, numerical algorithms described and static and dynamic fracture simulations with complex crack patterns are conducted to demonstrate the capability and advantage of the proposed approach.     

A hybrid experimental/numerical technique to extract cohesive fracture properties for mode-I fracture of quasi-brittle materials

We propose a hybrid technique to extract cohesive fracture properties of a quasi-brittle (not exhibiting bulk plasticity) material using an inverse numerical analysis and experimentation based on the optical technique of digital image correlation (DIC). Two options for the inverse analysis were used—a shape optimization approach, and a parameter optimization for a potential-based cohesive constitutive model, the so-called PPR (Park-Paulino-Roesler) model. The unconstrained, derivative free Nelder-Mead algorithm was used for optimization in the inverse analysis. The two proposed schemes were verified for realistic cases of varying initial guesses, and different synthetic and noisy displacement field data. As proof of concept, both schemes were applied to a Polymethyl-methacrylate (PMMA) quasi-static crack growth experiment where the near tip displacement field was obtained experimentally by DIC and was used as input to the optimization schemes. The technique was successful in predicting the applied load-displacement response of a four point bend edge cracked fracture specimen.     

A fracture framework for Euler–Bernoulli beams based on a full discontinuous Galerkin formulation/extrinsic cohesive law combination

A new full Discontinuous Galerkin discretization of Euler–Bernoulli beam is presented. The main interest of this framework is its ability to simulate fracture problems by inserting a cohesive zone model in the formulation. With a classical Continuous Galerkin method, the use of the cohesive zone model is difficult because inserting a cohesive element between bulk elements is not straightforward. On one hand if the cohesive element is inserted at the beginning of the simulation, there is a modification of the structure stiffness and on the other hand inserting the cohesive element during the simulation requires modification of the mesh during computation. These drawbacks are avoided with the presented formulation as the structure is discretized in a stable and consistent way with full discontinuous elements and inserting cohesive elements during the simulation becomes straightforward. A new cohesive law based on the resultant stresses (bending moment and membrane) of the thin structure discretization is also presented. This model allows propagating fracture while avoiding through‐the‐thickness integration of the cohesive law. Tests are performed to show that the proposed model releases, during the fracture process, an energy quantity equal to the fracture energy for any combination of tension‐bending loadings. Copyright © 2010 John Wiley & Sons, Ltd.     

Extended finite element method for cohesive crack growth

The extended finite element method allows one to model displacement discontinuities which do not conform to interelement surfaces. This method is applied to modeling growth of arbitrary cohesive cracks. The growth of the cohesive zone is governed by requiring the stress intensity factors at the tip of the cohesive zone to vanish. This energetic approach avoids the evaluation of stresses at the mathematical tip of the crack. The effectiveness of the proposed approach is demonstrated by simulations of cohesive crack growth in concrete.