Eigenvalue analysis of size effect for cohesive crack model

The paper analyses the effect of structure size on the nominal strength of the structure that is implied by the cohesive (or fictitious) crack model proposed for concrete by Hillerborg et al. A new method to calculate the maximum load of geometrically similar structures of different sizes without calculating the entire load-deflection curves is presented. The problem is reduced to a matrix eigenvalue problem, in which the structure size for which the maximum load occurs at the given (relative) length of the cohesive crack is obtained as the smallest eigenvalue. Subsequently, the maximum load, nominal strength and load-point displacement are calculated from the matrix equilibrium equation. The nonlinearity of the softening stress-displacement law is handled by iteration. For a linear softening law, the eigenvalue problem is linear and independent of the matrix equilibrium equation, and the peak load can then be obtained without solving the equilibrium equation. The effect of the shape of the softening law is studied, and it is found that the size effect curve is not very sensitive to it. The generalized size effect law proposed earlier by Baant, which describes a transition between the horizontal and inclined asymptotes of strength theory and linear elastic fracture mechanics, is found to fit the numerical results very well. Finally some implications for the determination of fracture energy from the size effect tests are discussed. The results are of interest for quasibrittle materials such as concrete, rocks, sea ice and modern tough ceramics.     

Effect of acid corrosion on crack propagation of concrete beams

The effect of acid corrosion on crack propagation of concrete beams was theoretically studied by the method of crack extension resistance curve. Based on this method, a calculation approach was proposed to determine fracture stress intensity factors in crack propagation of concrete beams. Loop iteration analysis was carried out to calculate maximum bearing capacity load, unstable crack toughness, resistance toughness curve, cohesive toughness curve and load–crack mouth opening displacement. Both bilinear and nonlinear softening traction–separation curves were adopted for each of these calculation parameters. The analysis results of each showed the effect of acid corrosion degrees. The influence of acid corrosion on fracture properties was discussed through the calculated results of cohesive toughness curves. These five kinds of simulated results were basically consistent, before the load attained the maximum value. However, with further crack propagation, cohesive toughness of nonlinear softening model was significantly larger than that of bilinear softening model, and the descending branch of P–CMOD curve by nonlinear law is higher than that by bilinear law. To validate the approach, tests of specimens under six different corrosion periods were experimentally studied, using three-point bending notched concrete beams soaked in sulphuric acid solution. The Double-K fracture parameters were investigated based on the test results, and load–crack mouth opening displacement curves for different acid conditions were obtained using synchronous sampling of a load sensor and clip-gauge. Numerical results by bilinear softening model showed a good correlation with the experimental ones.     

Determination of the kink point in the bilinear softening model for concrete

The characterization of the softening curve from experimental results is essential for predicting the fracture behavior of quasi-brittle materials like concrete. Among various shapes (e.g. linear, exponential) to describe the softening behavior of concrete, the bilinear softening relationship has been extensively used and is the model of choice in this work. Currently, there is no consensus about the location of the kink point in the bilinear softening curve. In this study, the location of the kink point is proposed to be the stress at the critical crack tip opening displacement. Experimentally, the fracture parameters required to describe the bilinear softening curve can be determined with the “two-parameter fracture model” and the total work of fracture method based on a single concrete fracture test. The proposed location of the kink point compares well with the range of kink point locations reported in the literature, and is verified by plotting stress profiles along the expected fracture line obtained from numerical simulations with the cohesive zone model. Finally, prediction of experimental load versus crack mouth opening displacement curves validate the proposed location of the kink point for different concrete mixtures and also for geometrically similar specimens with the same concrete mixture. The experiments were performed on three-point bending specimens with concrete mixtures containing virgin coarse aggregate, recycled concrete coarse aggregate (RCA), and a 50-50 blend of RCA and virgin coarse aggregate. The verification and validation studies support the hypothesis of the kink point occurring at the critical crack tip opening displacement.     

Eigenvalue method for computing size effect of cohesive cracks with residual stress,with application to kink-bands in composites

The previously developed eigenvalue method for computing the size effect of cohesive crack model is extended to the cohesive crack model with a finite residual stress. In this model, the structure size for which a specified relative length of kink-band corresponds to the maximum load is obtained as an eigenvalue of a homogeneous Fredholm integral equation. This new method is direct and much more efficient than the classical finite element approach in which the entire load-deflection history must be computed to obtain the maximum load. A secondary purpose of the paper is to apply the new method to the effect of structure size on the compressive strength of unidirectional fiber–polymer composites failing by propagation of kink-band with fiber microbuckling. The kink-band is simulated by a cohesive crack model with a linear compressive softening law and a finite residual stress. The simulation shows that the specimens tested have a negative–positive geometry, i.e., the energy release rate of the kink-band for a unit load first decreases but at a certain length of propagation begins to increase. Finally the effect of shape of the softening law of cohesive crack on the size effect curve is studied by using the new eigenvalue method. It is shown that, for a negative–positive geometry, the size effect on the peak load depends on the entire softening curve if the specimens is not too small.     

基于扩展有限元法的钢筋混凝土梁复合断裂过程模拟研究

针对钢筋混凝土梁裂纹扩展问题,基于扩展有限元法,建立了预置裂纹的简支混凝土梁三维模型,用粘聚裂纹模型描述裂纹面间的力学行为,采用线性的软化曲线表示裂纹尖端断裂过程区的应变软化行为,分别对素混凝土梁和钢筋混凝土梁的复合断裂过程进行模拟,分析了纵向钢筋对裂纹扩展路径、荷载-挠度和荷载 -CMOD (裂缝开口处张开位移)曲线的影响,并与文献中的试验结果进行对比,计算结果与试验结果吻合良好,展示了扩展有限元法在结构断裂破坏分析方面的独特优势。     

GENERALIZED SIZE EFFECT EQUATION FOR QUASIBRITTLE MATERIALS

Abstract— Structures made of quasibrittle materials show a marked decrease in strength as their size increases; this is the well known size effect on strength. This contribution proposes a general formula to predict with great accuracy the size effect of pre-cracked specimens. Such a formula provides the basis to determine experimentally the fracture parameters of a cohesive crack model from the measurement of the peak load in specimens with an initial crack. The formula is applicable whenever the fracture of the material can be described by a cohesive crack model with initial linear softening, as is generally accepted for concrete. Experimental evidence is presented showing that the formula is reasonably accurate also for rock.     

Influence of the fracture softening behaviour of wood on load-COD curve and R-curve

A cohesive crack model is used to analyse failure of wood in mode I along the grain. Several configurations of the gradual fracture softening behaviour of an interface, meshed with joint-elements located on the potential crack path, are investigated. Different constitutive laws, obtained from a single normalized polynomial function, are tested in order to estimate the influence of parameters such as, the tensile strength, the fracture energy or the ultimate opening of the interface, on the macroscopic response of a fracture specimen. Numerical results are compared with experimental data obtained on DCB specimen. We argue that the fracture energy related to the constitutive law must correspond to the plateau value of the R-curve. Moreover, this study reveals that the peak load of a load-COD (Crack Opening Displacement) curve is strongly affected by the slope of the softening behaviour. Finally, we present a review of the influence of each parameter describing the softening function on: (1) the load-COD curve and (2) the corresponding R-curve.     

Size effect and inverse analysis in concrete fracture

This paper summarizes the basic experimental and numerical results supporting an easy procedure to determine up to two fracture parameters based on numerically computed size effect curves. Furthermore, it supplies closed-form expressions to determine the initial linear segment approximation of the (stress vs. crack opening) softening curve of cohesive crack models for concrete, based only on the peak loads determined in splitting-tension (Brazilian) tests and in three-point-bending test on notched specimens. Knowledge of the initial segment, although not enough to describe all the fracture process of concrete structures, is enough to predict the fracture behavior of unnotched concrete structures prior and around the peak load. The same is true for notched structures provided their size is less than a limiting size, approximately defined in the paper. This revised version was published online in July 2006 with corrections to the Cover Date.     

Asymptotic analysis of a cohesive crack: 2. Influence of the softening curve

This paper presents a numerical method well suited to solve the integral equation governing the asymptotic behavior of a cohesive crack, and uses it to analyze the influence of the softening curve on the cracking response of large specimens. The analysis is performed with two main objectives in mind: (1) providing criteria to determine when a simplified linear elastic fracture mechanics (LEFM) approach can be applied, and (2) providing possible procedures of extracting information on the softening behavior from experimental data. The main conclusion is that the effective crack extension prior to peak is nearly determined by the length of the softening curve (the critical crack opening) and so is the deviation from LEFM. Furthermore, a simplified curve approach is proposed as an approximate alternative to solving the governing integral equation.     

Loading‐rate dependence of mode I crack growth in concrete

Mode I crack propagation process of concrete under relatively low loading rates which cover four orders of magnitude (0.2 μm/s to 2.0 mm/s) is investigated with three‐point bending (TPB) beams. All measured material properties exhibit rate sensitivity and follow a log‐linear relationship with the loading rate. A rate‐sensitive softening curve is established. The complete load‐crack mouth opening displacement (P‐CMOD) curve, crack propagation length, and fracture process zone (FPZ) length are simulated based on crack growth criterion with the fitted material parameters under those loading rates. Results show that the simulated P‐CMOD curves agree well with those of experimental measurements. It is clear that the peak load increases with the loading rate and so is the critical crack mouth opening displacement. Moreover, under the same load level, the length of the FPZ and the cohesive stress at the initial crack tip also increase with the increasing loading rate.     

Nonlinear fracture of cohesive materials

The cohesive crack is a useful model for describing a wide range of physical situations from polymers and ceramics to fiber and particle composite materials. When the cohesive zone length is of the order of the specimen size, the influence method—based on finite elements—may be used to solve the fracture problem. Here a brief outline of an enhanced algorithm for this method is given. For very large specimen sizes, an asymptotic analysis developed by the authors allows an accurate treatment of the cohesive zone and provides a powerful framework for theoretical developments. Some recent results for the zeroth order and first order asymptotic approaches are discussed, particularly the effective crack concept and the maximum load size effect. These methods are used to analyze the effect of the size and of the shape of the softening curve on the value at the peak load of several variables for three point bent notched beams. The results show, among other things, that for intermediate and very large sizes the size effect curves depend strongly on the shape of the softening curve, and that only the simultaneous use of asymptotic and influence methods may give an adequate estimate of the size effect in the intermediate range.     

Cohesive Fracture Model Based on Necking

A linear hardening model together with a linear elastic background material is first used to discuss some aspects of the mathematical and physical limitations and constraints on cohesive laws. Using an integral equation approach together with the cohesive crack assumption, it is found that in order to remove the stress singularity at the tip of the cohesive zone, the cohesive law must have a nonzero traction at the initial zero opening displacement. A cohesive zone model for ductile metals is then derived based on necking in thin cracked sheets. With this model, the cohesive behavior including peak cohesive traction, cohesive energy density and shape of the cohesive traction–separation curve is discussed. The peak cohesive traction is found to vary from 1.15 times the yield stress for perfectly plastic materials to about 2.5 times the yield stress for modest hardening materials (power hardening exponent of 0.2). The cohesive energy density depends on the critical relative plate thickness reduction at the root of the neck at crack initiation, which needs to be determined by experiments. Finally, an elastic background medium with a center crack is employed to re-examine the shape effect of cohesive traction–separation curve, and the relation between the linear elastic fracture mechanics (LEFM) and cohesive zone models by considering the cohesive zone development and crack growth in the infinite elastic medium. It is shown that the shape of the cohesive curve does affect the cohesive zone size and the apparent energy release rate of LEFM for the crack growth in the elastic background material. The apparent energy release rate of LEFM approaches the cohesive energy density when the crack extends significantly longer than the characteristic length of the cohesive zone.     

Asymptotic stress intensity factor density profiles for smeared-tip method for cohesive fracture

The paper presents a computational approach and numerical data which facilitate the use of the smeared-tip method for cohesive fracture in large enough structures. In the recently developed K-version of the smeared tip method, the large-size asymptotic profile of the stress intensity factor density along a cohesive crack is considered as a material characteristic, which is uniquely related to the softening stress-displacement law of the cohesive crack. After reviewing the K-version, an accurate and efficient numerical algorithm for the computation of this asymptotic profile is presented. The algorithm is based on solving a singular Abel's integral equation. The profiles corresponding to various typical softening stress-displacement laws of the cohesive crack model are computed, tabulated and plotted. The profiles for a certain range of other typical softening laws can be approximately obtained by interpolation from the tables. Knowing the profile, one can obtain with the smeared-tip method an analytical expression for the large-size solution to fracture problems, including the first two asymptotic terms of the size effect law. Consequently, numerical solutions of the integral equations of the cohesive crack model as well as finite element simulations of the cohesive crack are made superfluous. However, when the fracture process zone is attached to a notch or to the body surface and the cohesive zone ends with a stress jump, the solution is expected to be accurate only for large-enough structures.     

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.     

Damage distribution and size effect in numerical concrete from lattice analyses

Size effect on structural strength of concrete prisms subjected to three-point bending has been studied using the lattice model, which has been extended and now contains a realistic aggregate structure of concrete. The aggregate structure was obtained from CT-scans of real concrete prisms and overlaying the obtained image with a 3-dimensional hcp-lattice. The numerical analyses show that a size effect on structural strength exists for all studied aggregate densities and aggregate shapes. The size effect can be approximated with a Weibull model, where the main parameter, the Weibull modulus, depends on the concrete composition. The crack size distributions have been calculated and show a similar distribution as hypothesized before for fracture in ceramics. The results from the crack size distribution are helping to provide insight into the nature of the fracture process, which seems to differ from that hitherto assumed in cohesive crack models. After a weakening of the material through a multitude of microcracks, at peak load a single large crack propagates while loading continues in the softening regime. The presumed ‘cloud of microcracks’ advancing ahead of the macro-crack tip has not been found. Instead an alternative macroscopic model strategy, referred to as the 4-stage fracture model, is proposed.     

Fracture Analysis of Concrete Structural Components Accounting for Tension Softening Effect

This paper presents methodologies for fracture analysis of concrete structural components with and without considering tension softening effect. Stress intensity factor (SIF) is computed by using analytical approach and finite element analysis. In the analytical approach, SIF accounting for tension softening effect has been obtained as the difference of SIF obtained using linear elastic fracture mechanics (LEFM) principles and SIF due to closing pressure. Superposition principle has been used by accounting for non-linearity in incremental form. SIF due to crack closing force applied on the effective crack face inside the process zone has been computed using Green's function approach. In finite element analysis, the domain integral method has been used for computation of SIF. The domain integral method is used to calculate the strain energy release rate and SIF when a crack grows. Numerical studies have been conducted on notched 3-point bending concrete specimen with and without considering the cohesive stresses. It is observed from the studies that SIF obtained from the finite element analysis with and without considering the cohesive stresses is in good agreement with the corresponding analytical value. The effect of cohesive stress on SIF decreases with increase of crack length. Further, studies have been conducted on geometrically similar structures and observed that (i) the effect of cohesive stress on SIF is significant with increase of load for a particular crack length and (iii) SIF values decreases with increase of tensile strength for a particular crack length and load.     

Equivalence between stress intensity factor and energy approach based fracture parameters of concrete

The paper presents numerical study and relationship between the double-K fracture parameters and the double-G fracture parameters using two standard tests. The data required for calculation is obtained using cohesive crack model. It is observed that both the corresponding parameters of the double-K fracture model and the double-G fracture model at the onset of crack initiation and unstable fracture are equivalent. This observation agrees well with experimental results available in literature. It is also found that the fracture parameters of double-K fracture criterion and double-G fracture criterion are influenced by initial notch-length/depth ratio, specimen shape, size and softening function.     

Evaluation of Fracture Parameters by Double-G,Double-K Models and Crack Extension Resistance for High Strength and Ultra High Strength Concrete Beams

This paper presents the advanced analytical methodologies such as Double- G and Double - K models for fracture analysis of concrete specimens made up of high strength concrete (HSC, HSC1) and ultra high strength concrete. Brief details about characterization and experimentation of HSC, HSC1 and UHSC have been provided. Double-G model is based on energy concept and couples the Griffith's brittle fracture theory with the bridging softening property of concrete. The double-K fracture model is based on stress intensity factor approach. Various fracture parameters such as cohesive fracture toughness (K_Ic^c), unstable fracture toughness (K_Ic^un) and initiation fracture toughness (K_Icⁱⁿⁱ) have been evaluated based on linear elastic fracture mechanics and nonlinear fracture mechanics principles. Double-G and double-K method uses the secant compliance at the peak point of measured P-CMOD curves for determining the effective crack length. Bi-linear tension softening model has been employed to account for cohesive stresses ahead of the crack tip. From the studies, it is observed that the fracture parameters obtained by using double - G and double - K models are in good agreement with each other. Crack extension resistance has been estimated by using the fracture parameters obtained through double - K model. It is observed that the values of the crack extension resistance at the critical unstable point are almost equal to the values of the unstable fracture toughness K_Ic^un of the materials. The computed fracture parameters will be useful for crack growth study, remaining life and residual strength evaluation of concrete structural components.     

Cohesive Crack with Rate-Dependent Opening and Viscoelasticity: I. Mathematical Model and Scaling

The time dependence of fracture has two sources: (1) the viscoelasticity of material behavior in the bulk of the structure, and (2) the rate process of the breakage of bonds in the fracture process zone which causes the softening law for the crack opening to be rate-dependent. The objective of this study is to clarify the differences between these two influences and their role in the size effect on the nominal strength of stucture. Previously developed theories of time-dependent cohesive crack growth in a viscoelastic material with or without aging are extended to a general compliance formulation of the cohesive crack model applicable to structures such as concrete structures, in which the fracture process zone (cohesive zone) is large, i.e., cannot be neglected in comparison to the structure dimensions. To deal with a large process zone interacting with the structure boundaries, a boundary integral formulation of the cohesive crack model in terms of the compliance functions for loads applied anywhere on the crack surfaces is introduced. Since an unopened cohesive crack (crack of zero width) transmits stresses and is equivalent to no crack at all, it is assumed that at the outset there exists such a crack, extending along the entire future crack path (which must be known). Thus it is unnecessary to deal mathematically with a moving crack tip, which keeps the formulation simple because the compliance functions for the surface points of such an imagined preexisting unopened crack do not change as the actual front of the opened part of the cohesive crack advances. First the compliance formulation of the cohesive crack model is generalized for aging viscoelastic material behavior, using the elastic-viscoelastic analog (correspondence principle). The formulation is then enriched by a rate-dependent softening law based on the activation energy theory for the rate process of bond ruptures on the atomic level, which was recently proposed and validated for concrete but is also applicable to polymers, rocks and ceramics, and can be applied to ice if the nonlinear creep of ice is approximated by linear viscoelasticity. Some implications for the characteristic length, scaling and size effect are also discussed. The problems of numerical algorithm, size effect, roles of the different sources of time dependence and rate effect, and experimental verification are left for a subsequent companion paper. This revised version was published online in July 2006 with corrections to the Cover Date.     

Modelling cohesive crack growth using a two-step finite element-scaled boundary finite element coupled method

A two-step method, coupling the finite element method (FEM) and the scaled boundary finite element method (SBFEM), is developed in this paper for modelling cohesive crack growth in quasi-brittle normal-sized structures such as concrete beams. In the first step, the crack trajectory is fully automatically predicted by a recently-developed simple remeshing procedure using the SBFEM based on the linear elastic fracture mechanics theory. In the second step, interfacial finite elements with tension-softening constitutive laws are inserted into the crack path to model gradual energy dissipation in the fracture process zone, while the elastic bulk material is modelled by the SBFEM. The resultant nonlinear equation system is solved by a local arc-length controlled solver. Two concrete beams subjected to mode-I and mixed-mode fracture respectively are modelled to validate the proposed method. The numerical results demonstrate that this two-step SBFEM-FEM coupled method can predict both satisfactory crack trajectories and accurate load-displacement relations with a small number of degrees of freedom, even for crack growth problems with strong snap-back phenomenon. The effects of the tensile strength, the mode-I and mode-II fracture energies on the predicted load-displacement relations are also discussed.