Comparison of continuous and discontinuous constitutive models to simulate concrete behaviour under mixed‐mode failure conditions

The paper presents detailed FE simulation results of concrete elements under mixed‐mode failure conditions according to the so‐called shear‐tension test by Nooru‐Mohamed, characterized by curved cracks. A continuous and discontinuous numerical two‐dimensional approach was used. In order to describe the concrete's behaviour within continuum mechanics, two different constitutive models were used. First, an elasto‐plastic model with isotropic hardening and softening was assumed. In a compression regime, a Drucker–Prager criterion with a non‐associated flow rule was used. In turn, in a tensile regime, a Rankine criterion with an associated flow rule was adopted. Second, an isotropic damage constitutive model was applied with a single scalar damage parameter and different definitions of the equivalent strain. Both constitutive laws were enriched by a characteristic length of micro‐structure to capture properly strain localization. As an alternative approach, the extended finite element method was used. Our results were compared with the experimental ones and with results of other FE simulations reported in the literature. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of crack opening from a two‐dimensional continuum‐based finite element computation

Damage models are capable of representing crack initiation and mimicking crack propagation within a continuum framework. Thus, in principle, they do not describe crack openings. In durability analyses of concrete structures however, transfer properties are a key issue controlled by crack propagation and crack opening. We extend here a one‐dimensional approach for estimating a crack opening from a continuum‐based finite element calculation to two‐dimensional cases. The technique operates in the case of mode I cracking described in a continuum setting by a nonlocal isotropic damage model. We used the global tracking method to compute the idealized crack location as a post‐treatment procedure. The original one‐dimensional problem devised in Dufour et al. [4] is recovered as profiles of deformation orthogonal to the idealized crack direction are computed. An estimate of the crack opening and an error indicator are computed by comparing finite element deformation profiles and theoretical profiles corresponding to a displacement discontinuity. Two estimates have been considered: In the strong approach, the maxima of the profiles are assumed to be equal; in the weak approach, the integrals of each profile are set equal. Two‐dimensional numerical calculations show that the weak estimates perform better than do the strong ones. Error indicators, defined as the distance between the numerical and theoretical profiles, are less than a few percentages. In the case of a three‐point bending, test results are in good agreement with experimental data, with an error lower than 10% for widely opened crack (> 40µm). Copyright © 2011 John Wiley & Sons, Ltd.     

A macroscopic damage‐plastic constitutive law for modeling quasi‐brittle fracture and ductile behavior of concrete

A new phenomenological macroscopic constitutive model for the numerical simulation of quasi‐brittle fracture and ductile concrete behavior, under general triaxial stress conditions, is presented. The model is particularly addressed to simulate a wide range of confinement stress states, as also, to capture the strong influence of the mean stress value in the concrete failure mechanisms. The model is based on a two‐surface damage‐plastic formulation. The mechanical behavior in different domains of the stress space is separately described by means of a quasi‐brittle or ductile material response:

相似文献   

3‐D mesoscale simulation of crack‐permeability coupling in the Brazilian splitting test

In this paper, mesoscale hydromechanical simulations are performed to study (1) fracture features and (2) crack‐gas permeability coupling evolution in the context of the tensile splitting test. The mesostructure is based on a 2‐phase 3‐D representation of heterogeneous materials, such as concrete, where stiff aggregates are embedded into a mortar matrix. To take into account these heterogeneities without any mesh adaptation, a weak discontinuity is introduced into the strain field. In addition, a strong discontinuity is also added to take into account microcracking. This mechanical model is cast into the framework of the enhanced finite element method. Concerning the coupling with gas permeability, a double‐porosity method is used to simulate the flow through the cracks and the porosity. The apparent gas permeability is afterwards evaluated by a homogenization method. On the basis of finite element simulations, influence of aggregate size on ultimate crack opening, macroscopic ultimate tensile stress, total dissipated energy, and gas permeability evolution is numerically investigated. Furthermore, gas permeability evolution is also compared with experimental results from the literature. In addition, in the spirit of a sequential multiscale approach, macroscale gas permeability equations are identified from the hydromechanical results coming from the mesoscale computations. These equations lead to a relation between macroscale gas permeability evolution and crack opening. Besides, we show how the aggregate size influences the percolation threshold and that after this threshold, a cubic relation between macroscale gas permeability and crack opening is obtained.     

Simulation of localization failure with strain‐gradient‐enhanced damage mechanics

Strain gradient implies an important characteristic in localized damage deformation, which can be observed in the softening state of brittle materials, and strain gradients constitute the basic behaviours of localization failure area of the materials. The most important point in strain gradient is its damaging function including an internal length scale, which can be used to express the scale effects of mechanical responses of brittle rock mass. By extending the strain gradient theory and introducing an intrinsic material length scale into the constitutive law, the authors develop an isotropic damage model as well as a micro‐crack‐based anisotropic damage model for rock‐like materials in this paper. The proposed models were used to simulate the damage localization under uniaxial tension and plain strain compression, respectively. The simulated results well illustrated the potential of these models in dealing with the well‐known mesh‐sensitivity problem in FEM. In the computation, elements with C₁ continuity have been implemented to incorporate the proposed models for failure localization. When regular rectangle elements are encountered, the coupling between finite difference method (FDM) and conventional finite element method (FEM) is used to avoid large modification to the existing FEM code, and to obtain relatively higher efficiency and reasonably good accuracy. Application of the anisotropic model to the 3D‐non‐linear FEM analysis of Ertan arch dam has been conducted and the results of its numerical simulation coincide well with those from the failure behaviours obtained by Ertan geophysical model test. In this paper, new applications of gradient theories and models for a feasible approach to simulate localized damage in brittle materials are presented. Copyright © 2002 John Wiley & Sons, Ltd.     

Three‐dimensional simulations of tensile cracks in geomaterials by coupling meshless and finite element method

Failure in geotechnical engineering is often related to tension‐induced cracking in geomaterials. In this paper, a coupled meshless method and FEM is developed to analyze the problem of three‐dimensional cracking. The radial point interpolation method (RPIM) is used to model cracks in the smeared crack framework with an isotropic damage model. The identification of the meshless region is based on the stress state computed by FEM, and the adaptive coupling of RPIM and FEM is achieved by a direct algorithm. Mesh‐bias dependency, which poses difficulties in FEM‐based cracking simulations, is circumvented by a crack tracking algorithm. The performance of our scheme is demonstrated by two numerical examples, that is, the four‐point bending test on concrete beam and the surface cracks caused by tunnel excavation. Copyright © 2014 John Wiley & Sons, Ltd.     

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.     

A practical method to estimate crack openings in concrete structures

Numerical modelling of concrete cracking requires robust models able to describe opening and propagation of cracks. Structural concrete codes provide practical relations to describe crack openings. However, these empirical methods were developed for specific structures and cannot be used for general applications. Here, a continuous modelling approach based on damage mechanics is used to compute crack openings in a tie‐beam concrete structure. We propose a post‐processing method to extract crack openings from a continuum damage finite element computation. This method can be applied to all continuum damage/plasticity models. The tie‐beam concrete is characterized by a weak stress gradient; this aspect complicates predictions of crack positions and number. A stochastic method is used to take into account the spatial variability in concrete properties and create a spatially correlated random property field. Copyright © 2010 John Wiley & Sons, Ltd.     

Computational mechanics of the steel–concrete interface

Concrete cracking in reinforced concrete structures is governed by two mechanisms: the activation of bond forces at the steel–concrete interface and the bridge effects of the reinforcement crossing a macro‐crack. The computational modelling of these two mechanisms, acting at different scales, is the main objective of this paper. The starting point is the analysis of the micro‐mechanisms, leading to an appropriate choice of (measurable) state variables describing the energy state in the surface systems: on the one side the relative displacement between the steel and the concrete, modelling the bond activation; on the other hand, the crack opening governing the bridge effects. These displacement jumps are implemented in the constitutive model using thermodynamics of surfaces of discontinuity. On the computational side, the constitutive model is implemented in a discrete crack approach. A truss element with slip degrees of freedom is developed. This degree of freedom represents the relative displacement due to bond activation. In turn, the bridge effect is numerically taken into account by modifying the post‐cracking behaviour of the contact elements representing discrete concrete cracks crossed by a rebar. First simulation results obtained with this model show a good agreement in crack pattern and steel stress distribution with micro‐mechanical results and experimental results. Copyright © 2001 John Wiley & Sons, Ltd.     

Modelling crack propagation in concrete structures with a two scale approach

A simplified computational technique based on a refined global–local method is applied to the failure analysis of concrete structures. The technique distinguishes the scale of the structure, modelled with large size finite elements, from the scale at which material non‐linearity occurs due to progressive cracking and macro‐crack propagation. The finite element solution is split into two parts: a linear elastic analysis on a coarse mesh over the entire structure and a non‐linear analysis over a small part of the structure where a dense finite element grid is employed. In the non‐linear calculation, a non‐local damage model is implemented. These two computations are coupled with the help of an iterative scheme. The size and location of the region where a non‐linear analysis is performed, is adapted to follow the development of the damage zone. Numerical examples of mode I fracture of concrete specimens with straight and curved cracks are presented. Copyright © 2003 John Wiley & Sons, Ltd.     

An original semi‐discrete approach to assess gas conductivity of concrete structures

For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Damage‐viscoplastic consistency model for rock fracture in heterogeneous rocks under dynamic loading

This paper presents a damage‐viscoplastic consistency model for numerical simulation of brittle fracture in heterogeneous rocks. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept with separate damage variables in tension and compression. This approach does not suffer from ill‐posedness, caused by strain softening, of the underlying boundary/initial value problem since viscoplasticity provides the regularization by introducing a length scale effect under dynamic loading conditions. The model uses the Mohr–Coulomb yield criterion with the Rankine criterion as a tensile cut‐off. The damage law in compression is calibrated via the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition occurring in confined compression at a certain level of confinement. The heterogeneity of rock is accounted for by the statistical approach based on the Weibull distribution. Numerical simulations of confined compression test in plane strain conditions demonstrate a good agreement with the experiments at both the material point and structural levels as the fracture modes are realistically predicted. Copyright © 2009 John Wiley & Sons, Ltd.     

An anisotropic damage–plasticity model for saturated quasi‐brittle materials

This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd.     

Monotonic and cyclic modeling of interface between geotextile and gravelly soil

This paper describes a modified elasto‐plasticity damage model to capture monotonic and cyclic behavior of the interface between a geotextile and gravelly soil. New damage variable and shear strength criterion are introduced on the basis of test observations. The formulations of the modified model are obtained by extending those of the original interface model. The model parameters with physical meaning are easily determined from a group of cyclic shear tests and a confining compression test. The model predictions are compared with the results of a series of direct shear tests and large‐scale pullout tests. The comparison results demonstrate that the model accurately describes the monotonic and cyclic stress–strain relationship of the interface between a geotextile and gravelly soil while capturing new characteristics: (1) the strength that is nonlinearly dependent on the normal stress; (2) significant shear strain‐softening; (3) the comprehensive volumetric strain response with dependency on the shear direction; and (4) the evolution of behavior associated with the changes in the physical state that includes the geotextile damage. This model is used in a finite element analysis of pullout tests, indicating that the tensile modulus of a geotextile has a significant effect on the response of the geotextile–gravel system. Copyright © 2009 John Wiley & Sons, Ltd.     

考虑非线性徐变的混凝土温度裂缝扩展过程模拟

混凝土温度应力超过相应龄期的抗拉强度一半时,混凝土的徐变变形不再与应力呈线性关系。采用一个混凝土非线性徐变模型来进行混凝土的温度应力计算,并给出了相应的算法。同时,吸取无网格方法在裂缝扩展模拟方面的优势,采用无网格方法对混凝土的温度裂缝过程进行了仿真计算。数值算例计算结果表明,所编制的无网格法程序具有较好的收敛性和可靠性,能够有效地模拟混凝土块在温度应力作用下的开裂发展过程,考虑混凝土非线性徐变效应后,在开裂时刻混凝土急剧增加的非线性徐变变形缓和了缝端的温度拉应力,使得相同龄期混凝土的裂缝深度和扩展速率要比只考虑混凝土线性徐变要小一些。     

A multi‐asperity plastic‐contact crack plane model for geomaterials

A mechanistic constitutive model for fully formed cracks in geomaterials, such as concrete and rock, is presented. A three‐dimensional characterisation of the crack morphology is employed in which the crack surface is idealised as a series of conical teeth and corresponding recesses of variable height and slope. Based on this geometrical characterisation, an effective contact function is derived to relate the contact stresses that develop on the sides of the teeth to the net stresses on a crack plane. Plastic embedment and frictional sliding are simulated using a local plasticity model in which the plastic surfaces are expressed in terms of the contact surface function in cylindrical relative displacement space. Finally, the performance of the model is assessed against several sets of experimental data from direct shear tests, and it is concluded that the model is able to capture key characteristics of the behaviour of fully formed cracks in geomaterials. Copyright © 2012 John Wiley & Sons, Ltd.     

Rotational kinematic hardening model for three‐dimensional stress reversals in sand

A kinematic hardening mechanism has previously been proposed to capture the behavior of soil during large stress reversals in the triaxial plane. This mechanism is now extended to the principal stress space. It incorporates rotation and intersection of yield surfaces to achieve a consistent and physically rational fit with experimentally observed soil behavior during large three‐dimensional stress reversals. An existing elasto‐plastic model with isotropic hardening is used as the basic framework to which the rotational kinematic hardening mechanism has been added. The new combined model preserves the behavior of the isotropic hardening model under monotonic loading conditions, and the extension from isotropic to rotational kinematic hardening under three‐dimensional conditions is accomplished without introducing new material parameters. The framework of the model is presented here with some comparisons between theoretical and experimental directions of strain increment vectors to indicate the potential of the model. Copyright © 2008 John Wiley & Sons, Ltd.