Search

Fracture growth in a material is strongly influenced by the presence of inhomogeneities, which deviate crack trajectories from rectilinearity and deeply affect failure. Increasing crack tortuosity is connected to enhancement of fracture toughness, while often a crack may even be stopped when it impinges a void, which releases the stress concentration. Therefore, the determination of crack trajectories is important in the design against failure of materials and mechanical pieces. The recently developed phase-field approach (AT1 and AT2 models), based on a variational approach to damage localization, is believed to be particularly suited to describe complex crack trajectories. This belief is examined through a comparison between simulations and photoelastic experiments on PMMA plates, which have been designed in a new way, to highlight the effects of notches and circular holes on fracture propagation. The latter is shown to initiate from a notch and to be strongly attracted by voids. When a void is hit, fracture is arrested, unless the void contains a notch on its internal surface, from which a new crack nucleates and propagates. Different mechanical models are tested where fracture initiates and grows (i.) under Mode I compact tension, (ii.) four-point bending and (iii.) a tensile stress indirectly generated during compression of samples containing a circular hole. The experiments show that the fracture propagation may be designed to develop in different tortuous paths, involving multiple arrests and secondary nucleation. Simulations performed with an ad hoc implemented version of the AT1 and AT2 phase-field methods (equipped with spectral decomposition, in which a crack is simulated as a highly localized zone of damage accumulation) are shown to be in close agreement with experiments and therefore confirm the validity of the approach and its potentialities for mechanical design., In press on the International Journal of Solids and Structures

In this paper, the counter – current imbibition phenomena in a heterogeneous porous media is studied with the consideration of two types of porous materials like volcanic and fine sand and Adomian decomposition method is applied to find the saturation of wetting phase and the recovery rate of the reservoir. A simulation result is developed here to study the effect of heterogeneity, capillarity and relative permeability on saturation rate and to obtain an optimum recovery rate of the reservoir with the choices of some interesting parametric value.

In this paper, the counter-current imbibition phenomenon in two phase fluid through fracture porous media is discussed and Adomian decomposition method is applied to find the saturation of wetting phase and the recovery rate of the reservoir. A simulation result is developed for the saturation of wetting phase in fracture matrix and in porous matrix for some interesting choices of parametric value to study the recovery rate of the oil reservoir with dimensionless time. This problem has a great importance in the oil recovery process.

The goal of this work is to analyze a model for the rate-independent evolution of sets with finite perimeter. The evolution of the admissible sets is driven by that of a given time-dependent set, which has to include the admissible sets and hence is to be understood as an external loading. The process is driven by the competition between perimeter minimization and minimization of volume changes. \par In the mathematical modeling of this process, we distinguish the adhesive case, in which the constraint that the (complement of) the `external load' contains the evolving sets is penalized by a term contributing to the driving energy functional, from the brittle case, enforcing this constraint. The existence of Energetic solutions for the adhesive system is proved by passing to the limit in the associated time-incremental minimization scheme. In the brittle case, this time-discretization procedure gives rise to evolving sets satisfying the stability condition, but it remains an open problem to additionally deduce energy-dissipation balance in the time-continuous limit. This can be obtained under some suitable quantification of data. The properties of the brittle evolution law are illustrated by numerical examples in two space dimensions., Comment: Dedicated to Alexander Mielke on the occasion of his 60th birthday

We revisit the weak, energetic-type existence results obtained in (Rossi and Thomas, ESAIM Control Optim. Calc. Var. 21, 1–59, (2015)) for a system for rate-independent, brittle delamination between two visco-elastic, physically nonlinear bulk materials and explain how to rigorously extend such results to the case of visco-elastic, linearly elastic bulk materials. Our approximation result is essentially based on deducing the Mosco-convergence of the functionals involved in the energetic formulation of the system. We apply this approximation result in two different situations at small strains: Firstly, to pass from a nonlinearly elastic to a linearly elastic, brittle model on the time-continuous level, and secondly, to pass from a time-discrete to a time-continuous model using an adhesive contact approximation of the brittle model, in combination with a vanishing, super-quadratic regularization of the bulk energy. The latter approach is beneficial if the model also accounts for the evolution of temperature.

Recent experimental evidence on rubber has revealed that the internal cracks that arise out of the process, often referred to as cavitation, can actually heal. We demonstrate that crack healing can be incorporated into the variational framework for quasistatic brittle fracture evolution that has been developed in the last twenty years. This will be achieved for two-dimensional linearized elasticity in a topological setting, that is, when the putative cracks are closed sets with a preset maximum number of connected components. Other important features of cavitation in rubber, such as near incompressibility and the evolution of the fracture toughness as a function of the cumulative history of fracture and healing, have yet to be addressed even in the proposed topological setting.

In this paper, the counter - current imbibition phenomenon is discussed in an inclined heterogeneous porous media with the consideration of two types of porous materials like volcanic sand and fine sand. Adomian decomposition method is applied to find the saturation of wetting phase and the recovery rate of the reservoir. Finally, a simulation result is developed to study the saturation of wetting phase and the optimum recovery rate of reservoir with the choices of some interesting parametric values. This problem has a great importance in the field of oil recovery process.

This contribution presents recent results in the modeling and the analysis of delamination problems. It addresses adhesive contact, brittle, and cohesive zone models both in a quasistatic and a viscous, dynamic setting for the bulk part. Also different evolution laws for the delaminating surface are discussed.

Phase-field models have already been proven to predict complex fracture patterns in two and three dimensions for brittle fracture at small strains. In this paper we discuss a model for phase-field fracture at finite deformations in more detail. Among the identification of crack location and projection of crack growth the numerical stability is one of the main challenges in solid mechanics. We here present a phase-field model at finite strains, which takes into account the anisotropy of damage by applying an anisotropic split and the modified invariants of the right Cauchy-Green strain tensor. We introduce a suitable weak notion of solution that also allows for a spatial and temporal discretization of the model. In this framework we study the existence of solutions %Second the mathematical background of the approach is examined and and we show that the time-discrete solutions converge in a weak sense to a solution of the time-continuous formulation of the model. Numerical examples in two and three space dimensions are carried out in the range of validity of the analytical results.

We consider the weak solution of the Laplace equation in a planar domain with a straight crack, prescribing a homogeneous Neumann condition on the crack and a nonhomogeneous Dirichlet condition on the rest of the boundary. For every k we express the k-th derivative of the energy with respect to the crack length in terms of a finite number of coefficients of the asymptotic expansion of the solution near the crack tip and of a finite number of other parameters, which only depend on the shape of the domain.

The focus of this note lies on the numerical analysis of models describing the propagation of a single crack in a linearly elastic material. The evolution of the crack is modeled as a rate-independent process based on the Griffith criterion. We follow two different approaches for setting up mathematically well defined models: the global energetic approach and an approach based on a viscous regularization. &nbsp We prove the convergence of solutions of fully discretized models (i.e. with respect to time and space) and derive relations between the discretization parameters (mesh size, time step size, viscosity parameter, crack increment) which guarantee the convergence of the schemes. Further, convergence rates are provided for the approximation of energy release rates by certain discrete energy release rates. Thereby we discuss both, models with self-contact conditions on the crack faces as well as models with pure Neumann conditions on the crack faces. The convergence proofs rely on regularity estimates for the elastic fields close to the crack tip and local and global finite element error estimates. Finally the theoretical results are illustrated with some numerical calculations.

In this paper, we analyze a system for brittle delamination between two visco-elastic bodies, also subject to inertia, which can be interpreted as a model for dynamic fracture. The rate-independent flow rule for the delamination parameter is coupled with the momentum balance for the displacement, including inertia. This model features a nonsmooth constraint ensuring the continuity of the displacements outside the crack set, which is marked by the support of the delamination parameter. A weak solvabi- lity concept, generalizing the notion of energetic solution for rate-independent systems to the present mixed rate-dependent/rate-independent frame, is proposed. Via refined variational convergence techniques, existence of solutions is proved by passing to the limit in approximating systems which regularize the nonsmooth constraint by conditions for adhesive contact. The presence of the inertial term requires the design of suitable recovery spaces small enough to provide compactness but large enough to recover the information on the crack set in the limit.

Brittle Griffith-type delamination of compounds is deduced by means of Gamma-convergence from partial, isotropic damage of three-specimen-sandwich-structures by flattening the middle component to the thickness 0. The models used here allow for nonlinearly elastic materials at small strains and consider the processes to be unidirectional and rate-independent. The limit passage is performed via a double limit: first, we gain a delamination model involving the gradient of the delamination variable, which is essential to overcome the lack of a uniform coercivity arising from the passage from partial damage to delamination. Second, the delamination-gradient is supressed. Noninterpenetration- and transmission-conditions along the interface are obtained.

We discuss a model for crack propagation in an elastic body, where the crack path is described a-priori. In particular, we develop in the framework of finite-strain elasticity a rate-independent model for crack evolution which is based on the Griffith fracture criterion. Due to the nonuniqueness of minimizing deformations, the energy-release rate is no longer continuous with respect to time and the position of the crack tip. Thus, the model is formulated in terms of the Clarke differential of the energy, generalizing the classical crack evolution models for elasticity with strictly convex energies. We prove the existence of solutions for our model and also the existence of special solutions, where only certain extremal points of the Clarke differential are allowed.

We address the analysis of a model for brittle delamination of two visco-elastic bodies, bonded along a prescribed surface. The model also encompasses thermal effects in the bulk. The related PDE system for the displacements, the absolute temperature, and the delamination variable has a highly nonlinear character. On the contact surface, it features frictionless Signorini conditions and a nonconvex, brittle constraint acting as a transmission condition for the displacements. We prove the existence of (weak/energetic) solutions to the associated Cauchy problem, by approximating it in two steps with suitably regularized problems. We perform the two consecutive passages to the limit via refined variational convergence techniques.

A unilateral contact problem between elastic bodies at small strains glued by a brittle adhesive is addressed in the quasistatic rate-independent setting. The delamination process is modeled as governed by stresses rather than by energies. This leads to a specific scaling of an approximating elastic adhesive contact problem, discretized by a semi-implicit scheme and regularized by a BV-type gradient term. An analytical zero-dimensional example motivates the model and a specific local-solution concept. Two-dimensional numerical simulations performed on an engineering benchmark problem of debonding a fiber in an elastic matrix further illustrate the validity of the model, convergence, and algorithmical efficiency even for very rigid adhesives with high elastic moduli.

We consider the nonlocal formulation of continuum mechanics described by peridynamics. We provide a link between peridynamic evolution and brittle fracture evolution for a broad class of peridynamic potentials associated with unstable peridynamic constitutive laws. Distinguished limits of peridynamic evolutions are identified that correspond to vanishing peridynamic horizon. The limit evolution is associated with dynamic brittle fracture and satisfies a dynamic energy inequality expressed in terms of the kinetic energy of the motion together with a bulk elastic energy and a Griffith surface energy. It corresponds to the simultaneous evolution of elastic displacement and brittle fracture with displacement fields satisfying the wave equation inside the cracking domain. The wave equation provides the dynamic coupling between elastic waves and the evolving fracture path inside the media. The elastic moduli, wave speed and energy release rate for the evolution are explicitly determined by moments of the peridynamic influence function and the peridynamic potential energy., Comment: 26 pages, 2 figures, typos and references added

The nonlinear evolution problem for a crack with a kink in elastic body is considered. This nonlinear formulation accounts the condition of mutual non-penetration between the crack faces. The kinking crack is presented with the help of two unknown shape parameters of the kink angle and of the crack length, which minimize an energy due to the Griffith hypothesis. Based on the obtained results of the shape sensitivity analysis, solvability of the evolutionary minimization problem is proved, and the necessary conditions for the optimal crack are derived.

International audience; This work aims to take into account the influence of boundary defects on the behaviour till rupture of structures without any fine geometrical description of the domain. This is achieved by appealing to two approaches: an asymptotic analysis of Navier equations and strong discontinuity models. We present in this work a strategy to couple the two approaches in order to provide the analysis till rupture of the structure behavior. The approach is validated on an academic example.

We want to understand he concentration of damage in microfractured elastic media. Due to the different scallings of the volume and area (or area and length in two dimensions) the traditional method of homogenization using periodic arrays of cells seems to fail when applied to the Mumford-Shah functional and to periodically fractured domains. In the present paper we are departing from traditional homogenization. The main result implies the use of Mumford-Shah energies and leads to an explanation of the observed concentration of damage in microfractured elastic bodies.

We study the evolution of a single crack in an elastic body and assume that the crack path is known in advance. The motion of the crack tip is modeled as a rate-independent process on the basis of Griffith's local energy release rate criterion. According to this criterion, the system may stay in a local minimum before it performs a jump. The goal of this paper is to prove the existence of such an evolution and to shed light on the discrepancy between the local energy release rate criterion and models which are based on a global stability criterion (as for example the Francfort/Marigo model). We construct solutions to the local model via the vanishing viscosity method and compare different notions of weak, local and global solutions.

By a combination of geometrical and configurational analysis we study the properties of absolute minimal and equilibrium states of general Mumford-Shah functionals, with applications to models of quasistatic brittle fracture propagation. The main results concern the mathematical relations between physical quantities as energy release rate and energy concentration for 3D cracks with complex shapes, seen as outer measures living on the crack edge., Comment: larger version

Griffith's fracture criterion describes in a quasistatic setting whether or not a pre-existing crack in an elastic body is stationary for given external forces. In terms of the energy release rate (ERR), which is the derivative of the deformation energy of the body with respect to a virtual crack extension, this criterion reads: if the ERR is less than a specific constant, then the crack is stationary, otherwise it will grow. In this paper, we consider geometrically nonlinear elastic models with polyconvex energy densities and prove that the ERR is well defined. Moreover, without making any assumption on the smoothness of minimizers, we rigorously derive the well-known Griffith formula and the J-integral, from which the ERR can be calculated. The proofs are based on a weak convergence result for Eshelby tensors. Copyright © 2007 John Wiley & Sons, Ltd.

It is well known that high stress concentrations can occur in elastic composites in particular due to the interaction of geometrical singularities like corners, edges and cracks and structural singularities like jumping material parameters. In the project C5 "Stress concentrations in heterogeneous materials" of the SFB 404 "Multifield Problems in Solid and Fluid Mechanics" it was mathematically analyzed where and which kind of stress singularities in coupled linear and nonlinear elastic structures occur. In the linear case asymptotic expansions near the geometrical and structural peculiarities are derived, formulae for generalized stress intensity factors included. In the nonlinear case such expansions are unknown in general and regularity results are proved for elastic materials with power-law constitutive equations with the help of the difference quotient technique combined with a quasi-monotone covering condition for the subdomains and the energy densities. Furthermore, some applications of the regularity results to shape and structure optimization and the Griffith fracture criterion in linear and nonlinear elastic structures are discussed. Numerical examples illustrate the results.

The purpose of this paper is to fill the gap between the classical treatment of brittle fracture mechanics and the new idea of considering the crack evolution as a free discontinuity problem. Griffith and Irwin criterions of crack propagation are studied and transformed in order to be no longer dependent on any prescription of the geometry of the crack during its evolution. The inequality contained in theorem 6.1. represents the link between generalized Irwin and Griffith criterions of brittle crack propagation. The physical meaning of this inequality is explained in the last section. Further information available at http://irmi.epfl.ch/cag/buliga_bfrac.html ., Comment: unpublished paper from 1997, circulated in manuscript