Topological search of the crack pattern from a continuum mechanical computation
Introduction
A correct evaluation of crack path and opening is important in many concrete engineering applications.
For some type of structures, crack opening estimation is required in relation with permeability issues, in order to limit fluid releases. Examples are confinement vessels and cooling towers of nuclear power plants, dams and liquid natural gas reservoirs [41], [16], [38]. In any case, since concrete is exposed to environmental agents, the knowledge of crack characteristics may be exploited in order to assess durability, in the sense of penetration of external chemical aggressive agents [5] or simply for aesthetic reasons. Finally, even design codes prescribe as a service limit state a maximum crack opening which in turn must be addressed carefully.
The assessment of the mechanical state of structures is more and more often obtained by numerical simulations, the finite element method being by far the most employed technique. The mechanical simulation of a degraded or fractured structure can be performed at present days following different approaches.
A first approach, based on fracture mechanics, assumes the existence of a crack or a set of potential cracks. The crack is then modelled as a discontinuity, sometimes by describing it explicitly in the geometry [19], sometimes by enriching with Heaviside functions the interpolation of displacements in the element, as in the X-FEM method [28], [30], [29]. Due to the direct incorporation of the discontinuity, this approach is here referred to as discrete. Using such a method has the clear advantage, that crack properties (crack path and opening) are retrieved directly, because they are main variables of the model. It is also true that X-FEM has allowed to overcome some traditional shortcomings of discrete approaches, such as necessity of incremental remeshing [3] or constraints on crack path direction for cohesive elements, which otherwise must be known a priori [19]. However, this approach has still a limitation when it comes to describe appropriately the dissipative bulk behaviour (i.e. microcracking and kinematic field across the fracture process zone) on one hand, and crack initiation on the other hand.
A completely different approach, called here continuum approach, consists in adopting a nonlinear constitutive law, where the degradation state of the material is described by at least one state variable [2], [21], [22]. Through continuum models the description of both crack initiation and propagation in the same frame is possible. A further benefit of the continuum with respect to discrete approach lies in the correct description of the degradation process for quasi-brittle materials, from diffuse damage to formation of the macro-crack due to coalescence of micro-cracks. In other terms, in quasi-brittle materials energy dissipation takes place in a volume of non negligible size with respect to the characteristic size of the structure; this volume is usually called Fracture Process Zone (FPZ). Unfortunately, continuum approaches do not directly provide the crack path and opening.
A third alternative, more and more explored by researchers, consists in adopting both approaches at the same time. An embedded or explicitly meshed discontinuity models the macro-crack, while a non-linear constitutive law is used to model the diffuse damage [1], [32]. However, intrinsic difficulties may arise by establishing the transition point between the two kind of approaches whatever the mechanical quantity used for the equivalence [8], [42]. Some recent works provide arguments to overcome this problem by proposing a transition between a damage model and a cohesive zone approach [7], [24], but still there are limitations on the crack direction. Finally, the method of thick level sets [31] seems very promising in the next future.
The choice made in this contribution is to perform a continuum nonlinear computation based on classical finite element method and then retrieve the crack characteristics a posteriori in a post-processing phase. The complete procedure is two fold; in a first step the crack path is estimated and in the second step the crack opening along the found crack path is computed. For the latter phase, it is possible to establish a strain equivalence with a strong discontinuity [11]; another method consists in computing a displacement jump by isolating the inelastic part if any of the strain tensor concurring to open the crack, so ignoring the elastic part released upon unloading [26]. In order to obtain the crack path, a method has recently been developed in [10], inspired by the Global Tracking Method used in [33], [32]; the crack path is an isoline of a scalar field obtained from a secondary gradient problem. However, this method is, so far, limited to radial loading and mode I failure, since it is based on the assumption that, at the time the crack path is computed, the principal direction of the largest principal strain is perpendicular to the crack. Furthermore, this is not verified at the intersection between a FPZ and a rebar in concrete. The aim of this new contribution is thus to obtain the crack path in an automatic way for a wider range of mechanical problems.
In the first part of this paper, the method for tracking the crack path is presented. In a second stage, choices of numerical parameters are discussed. Thirdly, two continuum models used for illustrative purposes are briefly described. Finally, the powerfulness of the method is analysed on two testcases.
Section snippets
General discussion
The key idea is to search a crack path as a ridge in the topological three-dimensional space generated by the scalar field , i.e. a state variable field representing material degradation from a 2D simulation. This state variable can be an internal variable of a continuum model, such as isotropic damage or one component of an anisotropic damage model [9]; it can also be a plastic deformation or the hardening/softening variable, when plasticity models are used; eventually, it can also be an
Discussion on parameters range
Overall, the procedure is based on the use of five numerical parameters: the search length a, the smoothing length , the orthogonal line length , the point density on the orthogonal line and the value to arrest the search on the single FPZ.
Due to spatial discretization, the field is usually known only on a set of points on the considered domain Ω. Depending on the finite element solution, the field X has seldom an interpolation of order higher than linear. In this
Continuum models used in the examples
In this section, the continuum models employed for the illustrative applications are summarized. Two isotropic damage models have been used to highlight the versatility of the proposed method. Both are associated to two different regularisation techniques, since conventional continuum damage descriptions suffer from ill-posedness beyond a certain level of damage and lead to spurious localizations. A certain number of regularisation techniques exist in literature, such as the integral [36], [15]
Specimen with fixed circular inclusion
A first example is a squared specimen with a fixed circular inclusion. The specimen is 0.3 m by 0.3 m, the inclusion has a radius equal to 0.05 m and is located in the middle of the specimen. An imposed displacement is applied on the upper horizontal edge, while displacements of the inclusion are set to zero in all directions (see Fig. 9). The finite element code Code_Aster [4] is used. The specimen is computed up to its complete failure by means of the brittle damage model of Section 4.1. Law
Conclusion
In this contribution, an original method to track the crack path a posteriori from the results of continuum computation in the finite element frame is proposed and validated. The continuum model must be able to describe failure; for example, damage and plasticity models can be used for this purpose. A topological space is generated by a scalar field on a bidimensional domain, with X being the variable describing failure. The crack paths are identified as the ridges of this topological
Acknowledgements
Financial support was provided by the R&D department of French company EDF and by the French national “ANR Mefisto” grant. These supports are gratefully acknowledged. 3SR is part of the LabEx Tec 21 (Investissements d’Avenir – Grant agreement no. ANR-11-LABX-0030).
References (42)
- et al.
On the use of embedded discontinuity elements with crack path continuity for mode I and mixed mode fracture
Eng Fract Mech
(2002) - et al.
A general framework for continuum damage models. I. Infinitesimal plastic damage models in stress space
Int J Solids Struct
(2000) - et al.
Remeshing strategies for adaptive ALE analysis of strain localisation
Eur J Mech A/Solids
(2000) - et al.
A cohesive zone model which is energetically equivalent to a gradient-enhanced coupled damage-plasticity model
Eur J Mech A Solids
(2010) - et al.
Nonlocal anisotropic damage model and related computational aspects for quasi-brittle materials
Eng Fract Mech
(2007) - et al.
Damage, gradient of damage and principle of virtual power
Int J Solids Struct
(1996) - et al.
Strain-based transient-gradient damage model for failure analyses
Comput Methods Appl Mech Eng
(1998) - et al.
Modified nonlocal damage model based on stress state influence
Int J Solids Struct
(2011) - et al.
Containment evaluation under severe accidents (CESA): synthesis of the predictive calculations and analysis of the first experimental results obtained on the Civaux mock-up
Nucl Eng Des
(2001) - et al.
Air leakage characteristics in cracked concrete
Nucl Eng Des
(1995)
How to use damage mechanics
Nucl Eng Des
A variational formulation for nonlocal damage models
Comptes Rendus Mécanique
Convergence of a gradient damage model toward a cohesive zone model
Comptes Rendus Mécanique
Gradient damage models: toward full-scale computations
Comput Methods Appl Mech Eng
Extended finite element method for cohesive crack growth
Eng Fract Mech
Continuum approach to material failure in strong discontinuity settings
Comput Methods Appl Mech Eng
Theoretical and computational issues in modelling material failure in strong discontinuity scenarios
Comput Methods Appl Mech Eng
Prediction of air and steam leak rate through cracked reinforced concrete panels
Nucl Eng Des
Airsteam leakage through cracks in concrete walls
Nucl Eng Des
Three-dimensional liquid transport in concrete cracks
Int J Numer Anal Methods Geomech
Cited by (12)
An explicit representation of cracks in the variational phase field method for brittle fractures
2021, Computer Methods in Applied Mechanics and EngineeringCitation Excerpt :Elena Tamayo-Mas et al. [41] proposed a medial-axis-based approach to reconstruct crack path via damage variable. Marina Bottoni et al. [42] retrieve the crack path out of a continuum numerical model. Rangarajan et al. [43] proposed an algorithm for crack path in two-dimensional.
A discrete anisotropic damage constitutive law with an enhanced mixed-mode kinematics: Application to RC shear walls
2017, Engineering Fracture MechanicsCitation Excerpt :Thus, the mixed-mode formulation of the enhanced damage model is relevant. Classical tracking algorithm technics [31,52,53] in order to obtain good orientation of the crack and avoid such mesh bias effect have not been used here considering that welded mesh for concrete reinforcement enforced the direction of cracks. The developed model is tested with respect to its capacity in reproducing the anisotropic cracking pattern under complex loading and also crack closure effects.
Identifying the crack path for the phase field approach to fracture with non-maximum suppression
2016, Computer Methods in Applied Mechanics and EngineeringCitation Excerpt :Then, a number of interpolation points are determined from this bounding area, from which a cubic spline is built as the initial guess for the crack. This step is comparable to the approach introduced in [19] but ours ensures: (1) a more uniform distance between the identified points and also (2) a stable weighted average expression to determine these points rather than local maximum. Afterwards, an optimization step is undertaken to improve the crack path.
Reconstruct lower-dimensional crack paths from phase-field point cloud
2023, International Journal for Numerical Methods in EngineeringNumerical simulations of crack path control using soundless chemical demolition agents and estimation of required pressure for plain concrete demolition
2018, Materials and Structures/Materiaux et Constructions