EARTHQUAKE ANALYSIS OF GRAVITY DAMS BASED ON DAMAGE MECHANICS CONCEPT

In this paper, the seismic response analysis of concrete gravity dams is presented using the concept of Continuum Damage Mechanics. The analysis is performed using the finite element technique and a proper material degradation/damage model. The damage criterion used here is a second order tensor model based on elastic-brittle characterization and on a power function of the principal tensile stress. The methodology employed is shown to be computationally efficient and consistent in its treatment of both damage growth and propagation. Other important features considered in the analysis are: (1) dam–foundation interaction (2) appropriate modelling of joined rock mass using continuum damage mechanics, and (3) proper modelling of unbounded domain of foundation rock. The infinite media representation of the foundation material has been achieved by using doubly asymptotic approximation. The results of the analysis indicate that the seismic response of a damaged concrete dam could be significantly different from that of an undamaged one. In particular, the analysis shows that during a seismic event, the microstructure of a damaged zone can significantly change due to growth and propagation of microcracks.     

A technique to test finite difference schemes to model some geophysical processes in a geological structure

A preliminary problem to solve in the set-up of a mathematical model simulating a geophysical process is the choice of a suitable discrete scheme to approximate the governing differential equations. This paper presents a simple technique to test finite difference schemes used in the modeling of geophysical processes occurring in a geological structure. This technique consists in generating analytical solutions similar to the ones characterizing a geophysical process, given general information on some relevant parameters. Useful information for the choice of the discrete scheme to employ in the mathematical model simulating the original geophysical process can be obtained from the comparison between the analytical solution and the approximated numerical solutions generated by means of different discrete schemes. Two classes of numerical examples approximating the differential equation that governs the steady state earth's heat flow have been treated using three different finite differences schemes. The first class of examples deals with media whose phenomenological parameters vary as continuous space functions; the second class, instead, deals with media whose phenomenological parameters vary as discontinuous space functions. The finite difference schemes that have been utilized are: Centered Finite Difference Scheme (CDS), Arithmetic Mean Scheme (AMS), and Harmonic Mean Scheme (HMS).The numerical simulations showed that: the CDS may yield physically inconsistent solutions if the lattice internodal distance is too large, but in case of phenomenological parameters varying as a continuous function, this pitfall can be avoided increasing the lattice node refinement. In case of phenomenological parameters varying as a discontinuous function, instead, the CDS may yield physically inconsistent solutions for any lattice-node refinement. The HMS produced good results for both classes of examples showing to be a scheme suitable to model situations like these.     

Coupling between progressive damage and permeability of concrete: analysis with a discrete model

A significant increase of the permeability of concrete upon micro‐cracking and a good correlation between the evolution of damage (material stiffness) and permeability are observed experimentally. The present contribution investigates this correlation theoretically, with the help of lattice analyses. Scaling analysis of lattices which contain elastic brittle bonds has shown that the material degradation should be described by the evolution of the material stiffness, or compliance, in a continuum setting (damage models). This result is reviewed and further documented in the first part of the paper. In the second part, hydro‐mechanical problems are considered with the construction of a hydraulic lattice, dual to the mechanical one. We observe that the average permeability upon micro‐cracking is the lattice scale‐independent controlling variable in the hydraulic problem. Additionally, results show that in a continuum poro‐mechanical approach, the evolution of the material permeability ought to be related to the elastic unloading stiffness, described e.g. with the help of continuum damage variables. Copyright © 2005 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.     

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.     

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.     

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.     

On incremental non‐linearity in granular media: phenomenological and multi‐scale views (Part I)

On the basis of fundamental constitutive laws such as elasticity, perfect plasticity, and pure viscosity, many elasto‐viscoplastic constitutive relations have been developed since the 1970s through phenomenological approaches. In addition, a few more recent micro‐mechanical models based on multi‐scale approaches are now able to describe the main macroscopic features of the mechanical behaviour of granular media. The purpose of this paper is to compare a phenomenological constitutive relation and a micro‐mechanical model with respect to a basic issue regularly raised about granular assemblies: the incrementally non‐linear character of their behaviour. It is shown that both phenomenological and micro‐mechanical models exhibit an incremental non‐linearity. In addition, the multi‐scale approach reveals that the macroscopic incremental non‐linearity could stem from the change in the regime of local contacts between particles (from plastic regime to elastic regime) in terms of the incremental macroscopic loading direction. Copyright © 2005 John Wiley & Sons, Ltd.     

Modelling of anisotropic damage in brittle rocks under compression dominated stresses

A new model for describing induced anisotropic damage in brittle rocks is proposed. Although phenomenological, the model is based on physical grounds of micromechanical analysis. Induced damage is represented by a second rank tensor, which is related to the density and orientation of microcracks. Damage evolution is related to the propagation condition of microcracks. The onset of microcrack coalescence leading to softening behaviour is also considered. The effective elastic compliance of the damaged material is obtained from a specific form of Gibbs potential. Irreversible damage‐related strains due to residual opening of microcracks after unloading are also captured. All the model's parameters could be determined from conventional triaxial compression tests. The proposed model is applied to a typical brittle rock. Comparison between test data and numerical simulations shows an overall good agreement. The proposed model is able to describe the main features related to induced microcracks in brittle geomaterials. Copyright © 2002 John Wiley & Sons, Ltd.     

A Strain Rate Dependent Constitutive Model for the Lower Silurian Longmaxi Formation Shale in the Fuling Gas Field of the Sichuan Basin,China

Shale, as a kind of brittle rock, often exhibits different nonlinear stress-strain behavior, failure and time-dependent behavior under different strain rates. To capture these features, this work conducted triaxial compression tests under axial strain rates ranging from 5×10?6 s?1 to 1×10?3 s?1. The results show that both elastic modulus and peak strength have a positive correlation relationship with strain rates. These strain rate-dependent mechanical behaviors of shale are originated from damage growth, which is described by a damage parameter. When axial strain is the same, the damage parameter is positively correlated with strain rate. When strain rate is the same, with an increase of axial strain, the damage parameter decreases firstly from an initial value (about 0.1 to 0.2), soon reaches its minimum (about 0.1), and then increases to an asymptotic value of 0.8. Based on the experimental results, taking yield stress as the cut-off point and considering damage variable evolution, a new measure of micro-mechanical strength is proposed. Based on the Lemaitre’s equivalent strain assumption and the new measure of micro-mechanical strength, a statistical strain-rate dependent damage constitutive model for shale that couples physically meaningful model parameters was established. Numerical back-calculations of these triaxial compression tests results demonstrate the ability of the model to reproduce the primary features of the strain rate dependent mechanical behavior of shale.     

Propagation of a hydraulic fracture parallel to a free surface

This paper analyses the plane strain problem of a fracture, driven by injection of an incompressible viscous Newtonian fluid, which propagates parallel to the free surface of an elastic half‐plane. The problem is governed by a hyper‐singular integral equation, which relates crack opening to net pressure according to elasticity, and by the lubrication equations which describe the laminar fluid flow inside the fracture. The challenge in solving this problem results from the changing nature of the elasticity operator with growth of the fracture, and from the existence of a lag zone of a priori unknown length between the crack tip and the fluid front. Scaling of the governing equations indicates that the evolution problem depends in general on two numbers, one which can be interpreted as a dimensionless toughness and the other as a dimensionless confining stress. The numerical method adopted to solve this non‐linear evolution problem combines the displacement discontinuity method and a finite difference scheme on a fixed grid, together with a technique to track both crack and fluid fronts. It is shown that the solution evolves in time between two asymptotic similarity solutions. The small time asymptotic solution corresponding to the solution of a hydraulic fracture in an infinite medium under zero confining stress, and the large time to a solution where the aperture of the fracture is similar to the transverse deflection of a beam clamped at both ends and subjected to a uniformly distributed load. It is shown that the size of the lag decreases (to eventually vanish) with increasing toughness and compressive confining stress. Copyright © 2005 John Wiley & Sons, Ltd.     

Structural similarity and variety at the tips in a wide range of strike–slip faults: a review

Strike–slip faults are often accompanied by a variety of structures, particularly at their tips. The zones of additional fracturing are classified as tip‐damage zones. These zones can be subdivided into several different damage patterns based on the nature and orientation of faults and fractures developed. Damage zones at the ends of small strike–slip faults (mode II tips) develop wing cracks, horsetail splays, antithetic faults, synthetic branch faults and solution surfaces. Similar tip‐damage patterns are also commonly observed at larger (regional) scales, but with a dominance of faulting over tensile cracks and solution surfaces. Wing cracks and horsetail splays developed at small‐scale faults are replaced by normal faults in large‐scale faults. Antithetic faults and synthetic branch faults are observed at small and large scales. Thrust faults are developed at large scales, in a similar pattern to solution surfaces at a small scale. All these structures may show slightly different angular relationships to the master fault at small and large scale, but develop in response similar stress distribution and mechanics around the fault. Thus, mode II tip‐damage zones show similar patterns over a wide range of fault scales.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

Time dependence of axial pile response

This communication describes an approximate method for generalizing homogeneous, linearly elastic, subgrade type solutions of axial pile response to account for soil creep behaviour. The method is very simple to apply and can be used in conjunction with either analytical or numerical elastic solutions. Exact solutions for limiting cases and finite difference solutions in both space and time are presented to confirm that errors introduced by the approximations are small. Because of the theoretical basis of this approach it is anticipated that method can be used with other pile problems, lateral and axial, for both subgrade adn continuum idealizations. Methods for generalizing the results to more complex conditions such as non-homogeneity or time varying loads are given.     

A New Damage Index

A new damage index to estimate damage to buildings relies on construction costs per square metre, and a replacement ratio which approximates costs relative to the cost of replacing a median-sized family home. Building damage is estimated against a five-point scale with Central Damage Values at 0.02, 0.1, 0.4, 0.75 and 1.0 of the replacement cost.Damage is expressed as damage in House Equivalents (HE) = Replacement Ratio × Central Damage Value. The Damage Index = log₂ (HE) provides a simple 0–20 scale covering total damage of less than 1 HE to>1 million HE. For all natural hazard impacts in Australia DI is less than 12.Where the only damage data available are of lesser quality Generic or Qualitative Damage Indices (GDI and QDI) can be used. The various advantages and limitations of the Damage Index are discussed.     

Control‐volume mixed finite element methods

A key ingredient in simulation of flow in porous media is accurate determination of the velocities that drive the flow. Large‐scale irregularities of the geology (faults, fractures, and layers) suggest the use of irregular grids in simulation. This paper presents a control‐volume mixed finite element method that provides a simple, systematic, easily implemented procedure for obtaining accurate velocity approximations on irregular (i.e., distorted logically rectangular) block‐centered quadrilateral grids. The control‐volume formulation of Darcy’s law can be viewed as a discretization into element‐sized “tanks” with imposed pressures at the ends, giving a local discrete Darcy law analogous to the block‐by‐block conservation in the usual mixed discretization of the mass‐conservation equation. Numerical results in two dimensions show second‐order convergence in the velocity, even with discontinuous anisotropic permeability on an irregular grid. The method extends readily to three dimensions.     

Transient solution for a plane‐strain fracture driven by a shear‐thinning,power‐law fluid

This paper analyses the problem of a fluid‐driven fracture propagating in an impermeable, linear elastic rock with finite toughness. The fracture is driven by injection of an incompressible viscous fluid with power‐law rheology. The relation between the fracture opening and the internal fluid pressure and the fracture propagation in mobile equilibrium are described by equations of linear elastic fracture mechanics (LEFM), and the flow of fluid inside the fracture is governed by the lubrication theory. It is shown that for shear‐thinning fracturing fluids, the fracture propagation regime evolves in time from the toughness‐ to the viscosity‐dominated regime. In the former, dissipation in the viscous fluid flow is negligible compared to the dissipation in extending the fracture in the rock, and in the later, the opposite holds. Corresponding self‐similar asymptotic solutions are given by the zero‐viscosity and zero‐toughness (J. Numer. Anal. Meth. Geomech. 2002; 26 :579–604) solutions, respectively. A transient solution in terms of the crack length, the fracture opening, and the net fluid pressure, which describes the fracture evolution from the early‐time (toughness‐dominated) to the large‐time (viscosity‐dominated) asymptote is presented and some of the implications for the practical range of parameters are discussed. Copyright © 2006 John Wiley & Sons, Ltd.     

A parallel, implicit, cell‐centered method for two‐phase flow with a preconditioned Newton–Krylov solver

A new parallel solution technique is developed for the fully implicit three‐dimensional two‐phase flow model. An expandedcell‐centered finite difference scheme which allows for a full permeability tensor is employed for the spatial discretization, and backwardEuler is used for the time discretization. The discrete systems are solved using a novel inexact Newton method that reuses the Krylov information generated by the GMRES linear iterative solver. Fast nonlinear convergence can be achieved by composing inexact Newton steps with quasi‐Newton steps restricted to the underlying Krylov subspace. Furthermore, robustness and efficiency are achieved with a line‐search backtracking globalization strategy for the nonlinear systems and a preconditioner for each coupled linear system to be solved. This inexact Newton method also makes use of forcing terms suggested by Eisenstat and Walker which prevent oversolving of the Jacobian systems. The preconditioner is a new two‐stage method which involves a decoupling strategy plus the separate solutions of both nonwetting‐phase pressure and saturation equations. Numerical results show that these nonlinear and linear solvers are very effective.