strain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum princip

Similar documents
Milan Jirasek 1 Introduction In the context of standard continuum theories, strain-softening constitutive models typically lead to ill-posed boundary

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics

EVALUATION OF NONLOCAL APPROACHES FOR MODELLING FRACTURE IN NOTCHED CONCRETE SPECIMENS

Fluid driven cohesive crack propagation in quasi-brittle materials

NUMERICAL MODELLING AND DETERMINATION OF FRACTURE MECHANICS PARAMETERS FOR CONCRETE AND ROCK: PROBABILISTIC ASPECTS

ESCOLA POLITÉCNICA DA UNIVERSIDADE DE SÃO PAULO BOLETIM TÉCNICO PEF-EPUSP. Título:

Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach

INVESTIGATION OF MESH DEPENDENCE OF STOCHASTIC SIMULATIONS OF QUASIBRITTLE FRACTURE

BRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS

Finite element analysis of diagonal tension failure in RC beams

Abstract. 1 Introduction

Contemporary Research in Engineering Science

Modelling of ductile failure in metal forming

A DISCRETE NONLOCAL FORMULATION USING LOCAL CONSTITUTIVE LAWS

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

Cracked concrete structures under cyclic load

3D Finite Element analysis of stud anchors with large head and embedment depth

arxiv: v2 [cond-mat.mtrl-sci] 2 Jan 2009

Gradient-enhanced damage modelling of high-cycle fatigue

Numerical Characterization of Concrete Heterogeneity

Nonlocal model for size effect in quasibrittle failure based on extreme value statistics

MATHEMATICAL AND NUMERICAL MODEL OF ROCK/CONCRETE MECHANICAL BEHAVIOR IN A MULTI-PLANE FRAMEWORK

TENSILE CRACKING VIEWED AS BIFURCATION AND INSTABILITY IN A DISCRETE INTERFACE MODEL

A Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK

Size effect in the strength of concrete structures

Determination of size-independent specific fracture energy of concrete from three-point bend and wedge splitting tests

Equivalent localization element for crack band approach to mesh-sensitivity in microplane model

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS

MODELING OF THE WEDGE SPLITTING TEST USING AN EXTENDED CRACKED HINGE MODEL

Numerical modeling of standard rock mechanics laboratory tests using a finite/discrete element approach

Fracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design

Finite Element Solution of Nonlinear Transient Rock Damage with Application in Geomechanics of Oil and Gas Reservoirs

3D dynamic crack propagation analysis with PDS-FEM

Computational Modelling of Concrete Structures - Meschke, de Borst, Mang & Bieani Taylor & Francis Group, London, ISBN

A RATE-DEPENDENT MULTI-SCALE CRACK MODEL FOR CONCRETE

A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES

An anisotropic multi-plane elastic-damage model with axial and shear damage and its application to concrete behavior

International Journal Of Scientific & Engineering Research Volume 4, Issue 5, May-2013 ISSN

An orthotropic damage model for crash simulation of composites

5 ADVANCED FRACTURE MODELS

Microplane Model formulation ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary

CRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th

Mechanics (American Academy of Mechanics), Vol. 36, number 5-6 (May-June 2007), pp CONTRIBUTED ARTICLE

MODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN LAYERED STRUCTURES USING GEOMETRICALLY NONLINEAR BEAM FINITE ELEMENTS

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

Nonlocal computational methods applied to composites structures

NUMERICAL MODELLING AND DETERMINATION OF FRACTURE MECHANICS PARAMETERS FOR CONCRETE ROCK: INTRODUCTION

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface

MICROMECHANICAL MODELS FOR CONCRETE

CDPM2: A damage-plasticity approach to modelling the failure of concrete

Powerful Modelling Techniques in Abaqus to Simulate

Multiscale modeling of failure in ABS materials

TREBALL FINAL DE MÀSTER

MESO-MECHANICAL MODELING OF ULTRA-LIGHTWEIGHT CONCRETES

On the visco-elastic properties of open-cell polyurethane foams in uniaxial compression

Anisotropic Damage Mechanics Modeling of Concrete under Biaxial Fatigue Loading

A Method for Gradient Enhancement of Continuum Damage Models

Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering

Damage and crack modeling in single-edge and doubleedge notched concrete beams

Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

FRACTURING IN CONCRETE VIA LATTICE-PARTICLE MODEL

On the linearization of the microplane model

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations

SIZE EFFECTS IN THE BIAXIAL TENSILE-COMPRESSIVE BEHAVIOUR OF CONCRETE: PHYSICAL MECHANISMS AND MODELLING

An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation

NUMERICAL SIMULATION OF CONCRETE EXPOSED TO HIGH TEMPERATURE DAMAGE AND EXPLOSIVE SPALLING

Direct Comparison of Anisotropic Damage Mechanics to Fracture Mechanics of Explicit Cracks

The cohesive-frictional crack model applied to the bi-material interface between a dam and the foundation rock

Nigerian Journal of Technology, Vol. 26, No. 2, June 2007 Edelugo 37

FRACTURE IN HIGH PERFORMANCE FIBRE REINFORCED CONCRETE PAVEMENT MATERIALS

Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams

Configurational Forces as Basic Concepts of Continuum Physics

MESOSCOPIC MODELLING OF MASONRY USING GFEM: A COMPARISON OF STRONG AND WEAK DISCONTINUITY MODELS B. Vandoren 1,2, K. De Proft 2

NONLINEAR ANALYSIS OF A DAM-RESERVOIR-FOUNDATION SYSTEM UNDER SPATIALLY VARIABLE SEISMIC EXCITATIONS

Fracture Mechanics of Non-Shear Reinforced R/C Beams

Modeling issues of the FRP detachment phenomenon

HYDRAULIC FRACTURE PROPAGATION NEAR A NATURAL DISCONTINUITY

Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams

Transactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X

Heterogeneous structures studied by interphase elasto-damaging model.

SDM 2013 Student Papers Competition Modeling fiber-matrix splitting failure through a mesh-objective continuum-decohesive finite element method

EFFECT OF SHEAR REINFORCEMENT ON FAILURE MODE OF RC BRIDGE PIERS SUBJECTED TO STRONG EARTHQUAKE MOTIONS

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations

Numerical simulation of delamination onset and growth in laminated composites

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

ISSUES IN MATHEMATICAL MODELING OF STATIC AND DYNAMIC LIQUEFACTION AS A NON-LOCAL INSTABILITY PROBLEM. Ronaldo I. Borja Stanford University ABSTRACT

Loading σ Stress. Strain

436 A. Barani and G.H. Rahimi assessment models have been employed to investigate the LBB of cracked pipes that are not for combined load [8]. Yun-Jae

Fracture Scaling and Safety of Quasibrittle Structures: Atomistic Basis, Computational Challenges and New Advances

Module 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression

PERIDYNAMICS WITH ADAPTIVE GRID REFINEMENT

Coupling of plasticity and damage in glass fibre reinforced polymer composites

Fault Representation Methods for Spontaneous Dynamic Rupture Simulation. Luis A. Dalguer

THE MUTUAL EFFECTS OF SHEAR AND TRANSVERSE DAMAGE IN POLYMERIC COMPOSITES

Discrete Element Modelling of a Reinforced Concrete Structure

A Mesh Objective Algorithm for Modeling Mode-I Cracks Using a Standard Finite Element Formulation

Elsevier Editorial System(tm) for Computer Methods in Applied Mechanics and Engineering Manuscript Draft

MODELING OF NONLINEAR BEHAVIOR OF RC SHEAR WALLS UNDER COMBINED AXIAL, SHEAR AND FLEXURAL LOADING

Transcription:

Nonlocal damage models: Practical aspects and open issues Milan Jirasek LSC-DGC, Swiss Federal Institute of Technology at Lausanne (EPFL), Switzerland Milan.Jirasek@ep.ch Abstract: The purpose of this paper is to discuss certain diculties encountered when nonlocal constitutive models are applied in practical engineering analyses. Some of the problems remain unnoticed in theoretical studies but they may invalidate the results of large-scale simulations. 1. Introduction Quasibrittle materials, such as concrete, rock, tough ceramics, or ice, are characterized by the development of large nonlinear fracture process zones. The progressive growth of microcracks and their gradual coalescence can be described by constitutive laws with softening, i.e., with a descending branch of the stress-strain diagram. In the context of standard continuum mechanics, softening leads to serious mathematical and numerical diculties. The boundary value problem (BVP) becomes ill-posed, and numerical solutions exhibit a pathological sensitivity to the computational grid (nite-element mesh). Well-posedness of the BVP and objectivity of the numerical results can be restored by enriching the continuum theory. Regularization techniques that supply additional information on the internal structure of the material can enforce a realistic and mesh-independent size of the region of localized strain; therefore they are called localization limiters. A wide class of localization limiters is based on the concept of a nonlocal continuum. Generally speaking, the nonlocal approach consists in replacing a certain state variable by its nonlocal counterpart obtained by weighted averaging over a spatial neighborhood of each point under consideration. This idea originally appeared in elasticity (Eringen, 1966). Its early extensions to strain-softening materials, leading to the so-called imbricate continuum (Bazant, 1984), were later improved by the nonlocal damage theory (Pijaudier-Cabot and Bazant, 1987) and adapted for concrete (Saouridis and Mazars, 1992). Nonlocal formulations have been elaborated for a wide spectrum of models. Their detailed overview is out of scope of the present paper; the reader can nd a number of references in Bazant and Planas (1998). 2. Nonlocal Rotating Crack Model In the present study, we shall use the nonlocal version of the rotating crack model as formulated by Jirasek and Zimmermann (1998). This model can capture the cracking-induced anisotropic degradation of the material stiness while remaining relatively simple compared to general anisotropic damage formulations. The starting point is the standard smeared-crack concept based on a decomposition of the total strain into an elastic part and an inelastic part due to cracking. The cracking

strain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum principal stress. Crack opening and sliding contribute only to the strain normal to the crack and to the shear strain but not to the strain parallel to the crack. The normal and shear cracking strains are linked to the normal and tangential components of the traction transmitted by the crack. It is assumed either that the crack direction remains constant (xed-crack concept) or that it can vary (rotating-crack concept). Here we consider the latter case, in which the principal axes of stress and strain are kept aligned. The tangential components of cracking strain and crack traction remain zero, and the softening process is described by a simple relation between the normal traction across the crack and the inelastic strain due to crack opening. The crack opening strain is an internal variable that can be eliminated, either analytically (for a linear softening law), or numerically (for more general laws). The total stress-strain law can be written in the secant format resembling damage mechanics, = D(") " (1) where = stress, " = strain, and D(") = secant (damaged) material stiness. The nonlocal version of the model replaces (1) by where "(x) = = D(") " (2) Z V (x;)"() (3) is the nonlocal strain, and (x; ) = nonlocal weight function, decaying with the growing distance between points x and. The subsequent sections shall discuss certain diculties encountered when the nonlocal approach is applied to structural models with a general geometry and loading. The examples have been calculated with the nonlocal rotating crack model but many conclusions are valid for other damage-type models as well. 3. Spurious Shifting of Process Zone due to Gravity Forces The growth and coalescence of microcracks in the fracture process zone eventually lead to the formation of a stress-free macroscopic crack. Traditionally, this has been taken into account by softening laws for which the stress across the crack tends to zero, either at a nite value of cracking strain, or asymptotically. The latter case is more convenient from the numerical point of view, since the material keeps a nonzero stiness and the secant stiness matrix of the structure remains regular. A typical example is the exponential softening law = f t exp, " c " f where = normal stress transmitted by the crack, f t = uniaxial tensile strength, " c = cracking strain, and " f = material parameter related to the fracture energy.! (4)

The structure and evolution of the fracture process zone simulated with various nonlocal damage models were studied by Jirasek (1998a) on a simple one-dimensional test problem. It was shown that the zone of increasing local strain gets thinner as the loading process continues while the zone of increasing nonlocal strain keeps an approximately constant width. The secant material stiness and the residual strength tend to zero in a nite zone whose thickness corresponds to the support diameter of the nonlocal weight function. This zone must always contain several nite elements, otherwise the nonlocal interaction would not be properly captured by the numerical model. The degradation of the residual strength in a band containing several elements across its thickness causes serious numerical problems if the structure is subjected to distributed body forces. This is the case, e.g., in the analysis of gravity dams, where the eect of gravity forces is essential and cannot be neglected. The presence of gravity forces inside the fracture process zone leads to spurious shifting of the center of the zone and nally to the divergence of the equilibrium iteration process. This is documented in Fig. 1. (a) (b) 14 6 36 66 pressure full due to reservoir overflow pressure 72 (c) (d) Figure 1: Analysis of a gravity dam: (a) Geometry and loading, (b) nite element mesh, (c-d) evolution of the crack pattern (close-up area is marked in (a) by a dotted rectangle)

A gravity dam with dimensions corresponding to the well-known Koyna Dam is loaded by its own weight, full reservoir pressure, and an increasing hydrostatic pressure due to reservoir overow (Fig. 1a). The computational mesh, consisting of 4295 linear triangular elements, is shown in Fig. 1b. In a simulation using the nonlocal rotating crack model with exponential softening, the initially diuse cracking zone at the upstream face eventually localizes and propagates into the interior of the dam (Fig. 1c). Blue and red rectangles represent opening and closing smeared cracks, resp. When a critical degradation level is reached at the mouth of the macroscopic crack, the (constant) body forces start pulling the material down, which results into a shifting of the numerical crack (Fig. 1d) and later into the loss of convergence. Figure 2: Fracture process zone development and propagation of a discrete crack modeled as a displacement discontinuity embedded in nite elements (close-up area is marked in Fig. 1a by a dashed rectangle) An elegant remedy is a gradual transition from a smeared to a discrete description of material failure. Here we exploit the recently proposed combination of the nonlocal model with strong discontinuities embedded into the interior of nite elements (Jirasek, 1998b). This approach is appealing from the physical point of view. It is intuitively clear that diuse damage at early stages of material degradation is adequately described by a model dealing with inelastic strain while highly localized fracture is better modeled by a displacement discontinuity. It can also be argued that, as fracture progresses, long-range interaction between material points becomes more dicult and nally impossible. This would be best reected by a nonlocal

model with an evolving (decreasing) characteristic length. However, such a model would be computationally very expensive, since the interaction weights for all interacting pairs of integration points would have to be continuously recomputed. The approach proposed in Jirasek (1998b) can be considered as a reasonable approximation, with a constant interaction length in the early stage of degradation and a zero interaction length in the nal stage. The simulation of a gravity dam has been repeated using the model with transition from nonlocal rotating cracks to embedded displacement discontinuities (discrete cracks). The transition takes place when the cracking strain reaches a certain critical value, and the softening laws describing the smeared and discrete parts of the model are matched such that a smooth transition and correct energy dissipation are ensured. The initial cracking pattern is the same as before, but soon after localization the rst discrete crack appears and propagates along the center of the fracture process zone (see the white lines in Fig. 2). The partially damaged material in the wake of the discrete crack unloads, and its residual strength does not diminish anymore. Therefore, the correct crack trajectory can be captured without any numerical instabilities. 4. Pathological Interaction Between Tensile and Compressive Regions Another interesting issue related to nonlocal damage models is the spurious interaction between tensile and compressive regions. This phenomenon is illustrated by the example in Fig. 3. A rectangular sheet with a pre-existing slit-like crack is loaded by pressure applied on the crack faces. Only one fourth of the specimen is simulated; see Fig. 3a. Due to the specic type of loading, the region of high tensile stresses around the crack tip is very close to the compressed region around the crack faces. If the usual nonlocal formulation is used, two spurious interaction eects appear. Firstly, the nonlocal strains just at the crack tip can be smaller than at a certain distance from the crack tip, and so the rst numerical crack does not necessarily ap- (a) (b) (c) Figure 3: Cracked plate loaded on the crack faces: (a) geometry and loading (1/4 model), (b) initial stage of simulation (displacements exaggerated 500), (c) - nal stage of simulation (displacements exaggerated 100)

pear in the element closest to the tip (Fig. 3b). Secondly, as the macroscopic crack opens wide, the nonlocal strains in the compressed region can become positive (even though the local strains are of course negative), and cracking appears below the compressive load (Fig. 3c). It should be emphasized that this cracking is not primarily caused by the lateral expansion of the compressed material but by the proximity of a region that experiences large tensile strains. This spurious eect nally leads to a completely wrong failure mechanism instead of splitting of the specimen we obtain compressive failure at the faces of the initial crack. 5. Concluding Remarks We have briey described and illustrated two aspects of nonlocal damage models that are not apparent in simple academic examples, namely spurious shifting due to gravity forces and pathological interaction between tensile and compressive regions. Due to the limited size of this paper, other related issues are left for the oral presentation. Acknowledgment Financial support of the Swiss Committee for Technology and Innovation under project CTI.3201.1 is gratefully acknowledged. References Bazant, Z.P. (1984). "Imbricate continuum and its variational derivation." J. Eng. Mech. ASCE 110, 1693{1712. Bazant, Z.P., and Planas, J. (1998). Fracture and Size Eect in Concrete and Other Quasibrittle Materials. CRC Press, Boca Raton, FL. Eringen, A.C. (1966). "A unied theory of thermomechanical materials." Int. J. Eng. Sci. 4, 179{202. Jirasek, M. (1998a). "Comparison of nonlocal models for damage and fracture." LSC Internal Report 98/02, Swiss Federal Institute of Technology at Lausanne, Switzerland. Jirasek, M. (1998b). "Embedded crack models for concrete fracture." In Computational Modelling of Concrete Structures, ed. by R. de Borst et al., Balkema, Rotterdam, pp. 291{300. Jirasek, M., and Zimmermann, T. (1998). "Rotating crack model with transition to scalar damage." J. Eng. Mech. ASCE 124, 277{284. Pijaudier-Cabot, G., and Bazant, Z.P. (1987). "Nonlocal damage theory." J. Eng. Mech. ASCE 113, 1512{1533. Saouridis, C., and Mazars, J. (1992). "Prediction of the failure and size eect in concrete via a bi-scale damage approach." Eng. Computations 9, 329{344.