This paper addresses fluid-driven crack propagation in a porous medium. Cohesive interface elements are employed to model the behaviour of the crack. To simulate hydraulic fracturing, a fluid pressure degree of freedom is introduced inside the crack, separate from the fluid degrees of freedom in the bulk. Powell-Sabin B-splines, which are based on triangles, are employed to describe the geometry of the domain and to interpolate the field variables: displacements and interstitial fluid pressure. Due to their -continuity, the stress and pressure gradient are smooth throughout the whole domain, enabling a direct assessment of the fracture criterion at the crack tip and ensuring local mass conservation. Due to the use of triangles, crack insertion and remeshing are straightforward and can be done directly in the physical domain. During remeshing a mapping of the state vector (displacement and interstitial fluid pressure) is required. For this, a new methodology is exploited based on a least-square fit with the energy balance and mass conservation as constraints. The accuracy to model free crack propagation is demonstrated by two numerical examples, including crack propagation in a plate with two notches.

The intergranular strain concept was originally developed to capture the small-strain behaviour of the soil with hypoplastic models. A change of the deformation direction leads to an increase of the material stiffness. To obtain elastic behaviour for smallstrains, only the elastic part of the material stiffness matrix is used. Two different approaches for an application of this concept to nonhypoplastic models are presented in this article. These approaches differ in the determination of the elastic stress response, which is used for reversible deformations. The first approach determines an elastic response from the original material model, and the second one uses an additional elastic model. Both approaches are applied on barodesy. The simulations are compared with experimental results and with simulations using hypoplastic models with the original intergranular strain concept.

A nonlinear hybrid method is developed for multiscale analysis of a bearing-capacity test of a real-scale segmental tunnel ring subjected to point loads. The structural analysis consists of two parts. Part I refers to modeling of bending-induced tensile cracking of the segments, resulting from the external loading. The segments are subdivided into elements, according to the crack spacing. Each element is either intact or contains one central crack band, flanked by lateral undamaged domains. A multiscale model for tensile softening of concrete is used to describe the progressive deterioration of the crack bands. After iterative determination of their state of damage, the effective bending and extensional stiffnesses of the corresponding elements are quantified by means of Voigt-Reuss-Hill estimates. The effective stiffnesses are used for linear-elastic simulations of the segmental tunnel ring. Part II refers to the relative rotation angles at the joints, which are estimated from monitoring data, using the Bernoulli-Euler hypothesis. Since the validity of this hypothesis is questionable for neck-like joints, the relative rotation angles are post-processed such that they refer to rigid body displacements of the segments. The following conclusions are drawn: The presented approach yields good estimates of crack widths. Relative rotation angles at the joints mainly result in rigid body displacements of the segments, governing the convergences. Because realistic interface models are lacking, hybrid analysis based on displacement-monitoring data allows for performing ultimate-load analysis of segmental tunnel rings.

Powell-Sabin B-splines, which are based on triangles, are employed to model cohesive crack propagation without a predefined interface. The method removes limitations that adhere to isogeometric analysis regarding discrete crack analysis. Isogeometric analysis requires that the initial mesh be aligned a priori with the final crack path to a certain extent. These restrictions are partly related to the fact that in isogeometric analysis, the crack is introduced in the parameter domain by meshline insertions. Herein, the crack is introduced directly in the physical domain. Because of the use of triangles, remeshing and tracking the real crack path in the physical domain is relatively standard. The method can be implemented in existing finite element programmes in a straightforward manner through the use of Bézier extraction. The accuracy of the approach to model free crack propagation is demonstrated by several numerical examples, including discrete crack modelling in an L-shaped beam and the Nooru-Mohamed tension-shear test.

Discontinuity layout optimization (DLO) is a recently presented topology optimization method for determining the critical layout of discontinuities and the associated upper bound limit load for plane two-dimensional and three-dimensional (3D) problems. The modelling process (pre-processing) for DLO includes defining the discontinuities inside a specified domain and building the target function and the global constraint matrix for the optimization solver, which has great influence on the the efficiency of the computation processes and the reliability of the final results. This paper focuses on efficient and reliable pre-processing of the discontinuities within the 3D DLO and presents a multi-slicing strategy, which naturally avoids the overlapping and crossing of different discontinuities. Furthermore, the formulation of the 3D discontinuity considering a shape of an arbitrary convex polygon is introduced, permitting the efficient assembly of the global constraint matrix. The proposed method eliminates unnecessary discontinuities in 3D DLO, making it possible to apply 3D DLO for solving large-scale engineering problems such as those involving landslides. Numerical examples including a footing test, a 3D landslide and a punch indentation are considered, illustrating the effectiveness of the presented method. © 2016 The Authors. International published by John Wiley & Sons Ltd.

The purpose of this paper is to analyse the development and the evolution of the fracture process zone during fracture and damage in quasi-brittle materials. A model taking into account the material details at the mesoscale is used to describe the failure process at the scale of the heterogeneities. This model is used to compute histograms of the relative distances between damaged points. These numerical results are compared with experimental data, where the damage evolution is monitored using acoustic emissions. Histograms of the relative distances between damage events in the numerical calculations and acoustic events in the experiments exhibit good agreement. It is shown that the mesoscale model provides relevant information from the point of view of both global responses and the local failure process. © 2015 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons Ltd.

This paper describes a new method for representing concave polyhedral particles in a discrete element method as unions of convex dilated polyhedra. This method offers an efficient way to simulate systems with a large number of (generally concave) polyhedral particles. The method also allows spheres, capsules, and dilated triangles to be combined with polyhedra using the same approach. The computational efficiency of the method is tested in two different simulation setups using different efficiency metrics for seven particle types: spheres, clusters of three spheres, clusters of four spheres, tetrahedra, cubes, unions of two octahedra (concave), and a model of a computer tomography scan of a lunar simulant GRC-3 particle. It is shown that the computational efficiency of the simulations degrades much slower than the increase in complexity of the particles in the system. The efficiency of the method is based on the time coherence of the system, and an efficient and robust distance computation method between polyhedra as particles never intersect for dilated particles.