In its numerical implementation, the variational approach to brittle fracture approximates the crack evolution in an elastic solid through the use of gradient damage models. In the present paper, we first formulate the quasi-static evolution problem for a general class of such damage models. Then, we introduce a stability criterion in terms of the positivity of the second derivative of the total energy under the unilateral constraint induced by the irreversibility of damage. These concepts are applied in the particular setting of a one-dimensional traction test. That allows us to construct homogeneous as well as localized damage solutions in a closed form and to illustrate the concepts of loss of stability, of scale effects, of damage localization, and of structural failure. Considering several specific constitutive models, stress-displacement curves, stability diagrams, and energy dissipation provide identification criteria for the relevant material parameters, such as limit stress and internal length. Finally, the one-dimensional analytical results are compared with the numerical solution of the evolution problem in a two dimensional setting.
Phase-field models, sometimes refered to as gradient damage or smeared crack models, are widely used methods for the numerical simulation of crack propagation in brittle materials. Theoretical results and numerical evidences show that they can predict the propagation of a pre-existing crack according to Griffith' criterion. For a one-dimensional problem, it has been shown that they can predict nucleation upon a critical stress, provided that regularization parameter be identified with the material's internal or characteristic length. In this article, we draw on numerical simulations to study crack nucleation in commonly encountered geometries for which closed-form solutions are not available. We use U-and V-notches to show that the nucleation load varies smoothly from that predicted by a strength criterion to that of a toughness criterion, when the strength of the stress concentration or singularity varies. We present validation and verifications numerical simulations for both types of geometries. We consider the problem of an elliptic cavity in an infinite or elongated domain to show that variational phase field models properly account for structural and material size effects. We conclude that variational phase-field models can accurately predict crack nucleation through energy minimization in a nonlinear damage model instead of introducing ad-hoc criteria.
We study the genesis and the selective propagation of complex crack networks induced by thermal shock or drying of brittle materials. We use a quasistatic gradient damage model to perform large-scale numerical simulations showing that the propagation of fully developed cracks follows Griffith criterion and depends only on the fracture toughness, while crack morphogenesis is driven by the material's internal length. Our numerical simulations feature networks of parallel cracks and selective arrest in two dimensions and hexagonal columnar joints in three dimensions, without any hypotheses on cracks geometry, and are in good agreement with available experimental results.
Collagen fibres play an important role in the mechanical behaviour of many soft tissues. Modelling of such tissues now often incorporates a collagen fibre distribution. However, the availability of accurate structural data has so far lagged behind the progress of anisotropic constitutive modelling. Here, an automated process is developed to identify the orientation of collagen fibres using inexpensive and relatively simple techniques. The method uses established histological techniques and an algorithm implemented in the MATLAB image processing toolbox. It takes an average of 15 s to evaluate one image, compared to several hours if assessed visually. The technique was applied to histological sections of human skin with different Langer line orientations and a definite correlation between the orientation of Langer lines and the preferred orientation of collagen fibres in the dermis (p < 0.001, R(2) = 0.95) was observed. The structural parameters of the Gasser-Ogden-Holzapfel (GOH) model were all successfully evaluated. The mean dispersion factor for the dermis was κ = 0.1404±0.0028. The constitutive parameters μ, k(1) and k(2) were evaluated through physically-based, least squares curve-fitting of experimental test data. The values found for μ, k(1) and k(2) were 0.2014 MPa, 243.6 and 0.1327, respectively. Finally, the above model was implemented in ABAQUS/Standard and a finite element (FE) computation was performed of uniaxial extension tests on human skin. It is expected that the results of this study will assist those wishing to model skin, and that the algorithm described will be of benefit to those who wish to evaluate the collagen dispersion of other soft tissues.
a b s t r a c tWe consider a wide class of gradient damage models which are characterized by two constitutive functions after a normalization of the scalar damage parameter. The evolution problem is formulated following a variational approach based on the principles of irreversibility, stability and energy balance. Applied to a monotonically increasing traction test of a one-dimensional bar, we consider the homogeneous response where both the strain and the damage fields are uniform in space. In the case of a softening behavior, we show that the homogeneous state of the bar at a given time is stable provided that the length of the bar is less than a state dependent critical value and unstable otherwise. However, we also show that bifurcations can appear even if the homogeneous state is stable. All these results are obtained in a closed form. Finally, we propose a practical method to identify the two constitutive functions. This method is based on the measure of the homogeneous response in a situation where this response is stable without possibility of bifurcation, and on a procedure which gives the opportunity to detect its loss of stability. All the theoretical analyses are illustrated by examples.
The paper is devoted to gradient damage models which allow us to describe all the process of degradation of a body including the nucleation of cracks and their propagation. The construction of such model follows the variational approach to fracture and proceeds into two stages: (i) definition of the energy; (ii) formulation of the damage evolution problem. The total energy of the body is defined in terms of the state variables which are the displacement field and the damage field in the case of quasi-brittle materials. That energy contains in particular gradient damage terms in order to avoid too strong damage localizations. The formulation of the damage evolution problem is then based on the concepts of irreversibility, stability and energy balance. That allows us to construct homogeneous as well as localized damage solutions in a closed form and to illustrate the concepts of loss of stability, of scale effects, of damage localization, and of structural failure. Moreover, the variational formulation leads to a natural numerical method based on an alternate minimization algorithm. Several numerical examples illustrate the ability of this approach to account for all the process of fracture including a 3D thermal shock problem where the crack evolution is very complex.
SUMMARYThe variational approach to fracture is effective for simulating the nucleation and propagation of complex crack patterns but is computationally demanding. The model is a strongly nonlinear non-convex variational inequality that demands the resolution of small length scales. The current standard algorithm for its solution, alternate minimization, is robust but converges slowly and demands the solution of large, ill-conditioned linear subproblems. In this paper, we propose several advances in the numerical solution of this model that improve its computational efficiency. We reformulate alternate minimization as a nonlinear Gauss-Seidel iteration and employ over-relaxation to accelerate its convergence; we compose this accelerated alternate minimization with Newton's method, to further reduce the time to solution, and we formulate efficient preconditioners for the solution of the linear subproblems arising in both alternate minimization and in Newton's method. We investigate the improvements in efficiency on several examples from the literature; the new solver is five to six times faster on a majority of the test cases considered.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.