Associate Professor Klaus Thoeni

This work presents a novel scheme to couple the Boundary Element Method (BEM) and the Discrete Element Method (DEM) in the time domain. The DEM captures discontinuous material behaviour, such as fractured and granular media. However, applying the method to real-life applications embedded into infinite domains is challenging. The authors propose a solution to this challenge by coupling the DEM with the BEM. The capability of the BEM to model infinite domains accurately and efficiently, without the need for numerical artifices, makes it the perfect complement to the DEM. This study proposes a direct monolithic interface-based coupling method that resolves any incompatibilities between the two methods in two dimensions. The benchmark results show that the proposed methodology consistently produces results that align with analytical solutions. The final example in the paper showcases the full potential of this innovative methodology, where the DEM models a fracturing process, and the BEM evaluates its far-field effect.

This work presents a novel scheme to couple the Discrete Element Method (DEM) and the Boundary Element Method (BEM) for the multi-scale modelling in the time domain. The DEM can model discontinuous material at micro scale very well, but it cannot represent infinite domains. Hence, coupling with the BEM is proposed. A formulation employing the DEM and BEM in different subdomains of the same body is presented. There is no overlap between the sub-domains, and the system of equations is derived based on strong equilibrium and compatibility conditions at the interface. The proposed coupling scheme is based on monolithic time integration. The conducted numerical experiments of one-dimensional wave propagation show excellent agreement with the analytical solution. Some spurious wave reflections are observed at the interface, but their effect is quantified and deemed not critical for infinite domains, which are of main interest. Even though the applications for one-dimensional wave propagation are of limited practical engineering interest, this work represents a significant theoretical breakthrough. This paper establishes a reference for future coupling schemes for two- and three-dimensional multi-scale analysis.

This work presents an efficient probabilistic framework for the Bayesian calibration of micro-mechanical parameters for Discrete Element Method (DEM) modelling. Firstly, the superior behaviour of the iterative Bayesian filter over the sequential Monte Carlo filter for calibrating micro-mechanical parameters is shown. The linear contact model with rolling resistance is used for simulating the triaxial responses of Toyoura sand under different confining pressures. Secondly, synthetic data from DEM simulations of triaxial compression are used to assess the reliability of iterative Bayesian filtering with respect to the user-defined parameters, such as the number of samples and predefined parameter ranges. Excellent calibration results with errors between 1 and 2% are obtained when the number of samples is chosen high enough. It is crucial that the sample size is representative for the distribution of individual parameters within the predefined parameter ranges. The wider the ranges, the more samples are required. The investigation also shows the necessity of including both stress and strain histories, at certain confidence levels, for estimation of the correct mechanical responses, especially the correct fabric responses. Finally, based on the findings of this work a fully-automated open-source calibration tool is developed and demonstrated for selected stress paths.

Appropriate rock characterization is beneficial in providing a reliable judgment on rock properties which is crucial for the design process of rock engineering applications. However, it can be difficult to obtain consistent mechanical parameters due to substantial variations in rock properties. In this research, uniaxial compression tests on dolerite specimens collected from a gold mine in Western Australia showed substantial scatter in the results. Rock categorization based on the P-wave velocities is as accurate as the thin section analysis, which suggests that they can be used together to gain a more accurate initial understanding of the rock types before any laboratory testing. The quality of specimen preparation and rock¿machine interaction greatly affect the test results. For instance, non-parallelness of loading platens can lead to considerable scatter of the testing results, which would be perceived as rock variability. It is suggested that the current testing standards should be modified towards a better control of the loading machine performance and equipment precision. Finally, the possibility of pre-existing microcracks in rock, neither detected by the thin section analysis nor by the ultrasonic measurement, must be examined by computed tomography (CT) scanning as they can affect the test results. This study will enhance our knowledge about the sources of variability in laboratory test results of rock which is essential for obtaining reliable data.

This work presents a 3D extension of the deformable discrete element method (DDEM) developed previously for 2D problems. The 3D formulation employs spherical particles. The particle deformation is made up of a global and local deformation mode. The global mode is assumed to be produced by uniform stress due to the contact forces. Particle deformability yields a nonlocal contact model, in which one contact between particles is influenced by contacts with other particles. It also leads to the formation of new contacts in the particle assembly. The DDEM affects the behavior of the granular material at the macroscopic level and gives new possibilities in material modeling by the discrete element method (DEM). The new algorithm is verified on a unconfined uniaxial compression test of a cuboid specimen discretized with equal-size bonded particles aligned in a simple cubic pattern using an analytical solution. Enhanced modeling capabilities are presented by simulating cylindrical specimens discretized with a nonuniform size of bonded particles. The micro¿macro relationships for elastic parameters are obtained. It is shown that the DDEM extends the range of the Poisson's ratio achievable with the DEM. Additional simulations are performed to determine the stability limits of the DDEM.

Digital surface models (DSM) have become one of the main sources of geometrical information for a broad range of applications. Image-based systems typically rely on passive sensors which can represent a strong limitation in several survey activities (e.g., night-time monitoring, underground survey and night surveillance). However, recent progresses in sensor technology allow very high sensitivity which drastically improves low-light image quality by applying innovative noise reduction techniques. This work focuses on the performances of night-time photogrammetric systems devoted to the monitoring of rock slopes. The study investigates the application of different camera settings and their reliability to produce accurate DSM. A total of 672 stereo-pairs acquired with high-sensitivity cameras (Nikon D800 and D810) at three different testing sites were considered. The dataset includes different camera configurations (ISO speed, shutter speed, aperture and image under-/over-exposure). The use of image quality assessment (IQA) methods to evaluate the quality of the images prior to the 3D reconstruction is investigated. The results show that modern highsensitivity cameras allow the reconstruction of accurate DSM in an extreme low-light environment and, exploiting the correct camera setup, achieving comparable results to daylight acquisitions. This makes imaging sensors extremely versatile for monitoring applications at generally low costs.

Rockfall is a natural hazard that needs to be rigorously managed along all the major road and railways transport networks by identifying the most appropriate mitigation measures. There has been significant progress in rockfall modelling and rockfall protection systems in recent years but there remains one aspect that is not very well understood and quite challenging to account for in the design of rockfall protection structures, namely the fragmentation of falling blocks upon impact. Rocks often break up upon impact, which leads to a change in size, shape and energy of falling blocks, parameters that affect the design of the protective structures. Before being able to incorporate fragmentation into predictive trajectory models, it is required to better understand the fragmentation process and its likely outcome (number, volume of fragments and their trajectories). To that aim, an innovative experimental setup was developed at the University of Newcastle (Australia) to study rock fragmentation upon impact. The setup was designed to perform controlled vertical drop tests and record the following impact parameters: impact force, impulse, impact duration, velocities (of the block before impact and its fragment after impact) and all components of energy, pre and post impact. Six views (four high-speed cameras and two mirrors) are used for an accurate reconstruction of the 3D trajectory of blocks and fragments, in translation and rotation. This paper presents the validation of the setup via two series of drop tests using mortar spheres. Attention was focused on the evaluation of impact force and impulse from load cells placed under the impacted surface, tracking of translational and rotational velocity and the computation of total kinetic energy (before and after impact) and all components of energy dissipation. The results confirm that the experimental setup and the approach developed can be used to obtain impact force, impulse and to compute the energy balance during the impact and fragmentation and conduct advanced fragmentation testing.

Fragmentation of rocks upon impact during rockfall is a phenomenon that is poorly understood, scarcely researched and difficult to predict. However, to adequately predict the outcomes of rockfall events, it is essential to know whether a given block is likely to fragment given the impact conditions and what will be the outcome of the fragmentation process; that is, the number, size and trajectory of fragments. To date, there is no model or data that can be used to fully answer these questions. This paper presents the first theoretical model that can predict the fragmentation survival probability of brittle spherical blocks upon dynamic impact (i.e. drop tests) from the statistical distribution of material properties, determined from a range of standard quasi-static tests. Considering that survival probabilities tend to follow a Weibull distribution, the model focuses on predicting the two Weibull parameters, commonly known as the shape parameter (m) and the scale parameter (here, the critical kinetic energy). The model is based on theoretically-derived conversion factors used to turn the critical work required to fail disc samples in quasi-static indirect tension into the critical kinetic energy to cause failure of spheres at impact in drop tests. The mechanistic conversion factors specifically account for the shape and size of the specimens tested and the increase of strength under dynamic loading (strain rate effect). Three series of drop tests were conducted (on spheres of three different diameters) and complemented by extensive material characterisation testing in order to validate the novel predictive model. The variability of material properties was characterised, and it was found that the material strength found by the characterisation tests generally follows a Weibull form, but the survival probability distribution of the drop tests seems to be linear. The predicted conversion factors were first compared against their experimental counterparts before validating the prediction of survival probability of the spheres upon dynamic impact (in drop tests). It was found that it is possible to predict the survival probability of artificial rock of three different diameters (50 mm, 75 mm, 100 mm) and two different strengths upon impact solely from the statistical information coming from Brazilian tests and with an average relative error of less than 9%.

This paper presents a new methodology for estimating the expected energies and first impact distances at the base of a rock cliff, subject to the geometry and properties of the cliff and the representative block being known. The method is based on a sensitivity analysis, conducted by means of kinematic simulations and carried out for a large range of input parameters and their combinations, taking into account the uncertainty associated with their estimate. The proposed approach is validated by comparing predictions to experimental data and shows great potential for a quick qualitative hazard assessment.

This paper investigates the factors affecting the shear strength of granular materials by means of the tilt box test. The test involves the titling of a specifically designed split-box and measuring the angle at which the upper half of the box slides off. This angle can be interpreted as the peak friction angle of the sample for a low stress state. Experimental tests and numerical analyses using the Discrete Element Method (DEM) are carried out and compared. The study is carried out for two different materials, a rounded and an angular aggregate. It is shown that rolling resistance is needed to match the experimental results if only spheres are used but this is not the case for non-spherical particles. © 2015 Taylor & Francis Group, London.

This work consists of a preliminary numerical study on the impact behaviour of rock blocks on muckpiles using the Discrete Element Method (DEM). A contact law using a linear hysteresis model for normal contacts and a Mohr-Coulomb model for tangential contacts is used. In addition, two dissipative laws are introduced consisting of a moment transfer law and normal viscous damping law. One of the aims of this study is to investigate if the contact law enables a realistic modelling of the rebound characteristics of a block during impact on a muckpile. A parametric study is carried out in order to investigate the influence of the constitutive model parameters and the presence of clumps. Numerical predictions are compared to experimental results. It is shown that clumps in the muckpile are essential to reproduce realistic results. © 2015 Taylor & Francis Group, London.

This paper presents a discrete modelling approach which allows the simulation of rockfalls behind drapery systems. The falling rock is represented by a rigid assembly of spheres whereas the slope is represented by triangular elements. Energy dissipation during impact on the slope is considered via friction and viscous damping. The drapery is represented by a set of spherical particles which interact remotely. The numerical model is used to investigate the efficiency of a drapery system. Various simulations with two different block sizes are performed and the numerical predictions of simulations with drapery are compared to simulations without drapery in order to assess the performance of the protection system.