Professor Anna Giacomini

This study investigates the dynamic behavior between a rockfall and a new anchor cable ribbed rockfall retaining walls on a rock shed, which is designed to meet project requirements due to the tremendous space occupation of the protection structure, via numerical simulations by the coupling between PFC3D (discrete element method, DEM) and FLAC3D (finite element method, FEM). In the model, the slope model was imported after the point cloud of the slope was processed. The novel structure is modeled by the finite element method (FEM) through the zone element, while the rockfall and the buffer layer are modeled by the discrete element method (DEM) through the ball element. The dynamic movement of rockfall was traced, and the impact position and velocity were obtained on the structure. The numerical results show that three stages of the rockfall movement were modeled, namely, movement, impact, and stagnation. The rockfall falls on the slope, impacting the buffer layer on the top of the rock shed, before rebounding to the anchored rockfall retaining walls (with a height lower than 2 m). Then, the stress and deflection can be unified and related to the impact velocity of rockfall to examine the stability of the structure. In addition, the rockfall radius is the dominant parameter in the three parameters (rockfall shape, rockfall radius, and impact velocity). While the study focuses on a specific case study, the results provide valuable guidelines for future applications of the proposed combined structure for railway transportation protection.

Reliable shear strength determination of large in situ discontinuities is still a challenge faced by the rock mechanics field. This is principally due to the limited availability of surface roughness and morphology information of in situ discontinuities and the unresolved management of the ¿scale effect¿ phenomenon. Recently, a stochastic approach for predicting the shear strength of large-scale discontinuities was established, encompassing random field theory, a semi-analytical shear strength model, and a stochastic analysis framework. A key aspect of the new approach is the application at field scale, thereby minimising or bypassing the scale effect. The approach has been validated at laboratory scale and an initial large-scale deterministic-based validation showed promising results. However, to date, no large-scale experimental-based validation has been undertaken. This paper presents the first rigorous application of the employed semi-analytical shear strength model and the stochastic approach on a 2¿m-by-2¿m discontinuity surface, with comparison of prediction to experimental shear strength data. The shear strength model was found to generally produce peak and residual predictions within a ± 10% relative error range, with good agreement between predicted and observed damage areas. It was observed that, applying the stochastic approach to seed traces with gradient statistics equivalent to that of the surface, produced predictions that closely resemble the experimental results. Whereas, predicting shear strength from different seed traces results in more variability of predictions, with many falling within ± 20% of the experimental data. The predictions of residual shear strength tended to be more accurate than peak shear strength.

In the Load and Resistance Factor Design (LRFD) of piles, several design codes recommend higher resistance factors if load tests are conducted. However, no information is provided on how these resistance factors are determined. In this paper, a probabilistic approach based on Bayes' theorem and the First Order Reliability Method (FORM) is proposed to calibrate resistance factors for different numbers of load tests and the corresponding test results. In addition, within-site variability, design methods, types of piles, and ground conditions can also be considered. The proposed approach applied to H-piles under axial load tests shows consistent results with current design codes. Results show that resistance factors are significantly increased even if only one positive test is observed among all the tests. For low variability sites, the differences of resistance factors between various design methods are significantly reduced if one or more tests are positive, while for high variability sites, the differences of resistance factors are only slightly decreased, indicating that design methods should be considered in the latter case. Most of the increase in resistance factors is achieved with a small number of tests. For ß-Method used in clay sites, 80% of the increase in resistance factors is achieved with two, four and five consecutive positive tests are observed for low, medium and high variability sites, respectively.

The number of women employed in STEM in Australia is increasing, however, they continue to remain underrepresented in most industries. A significant corpus of literature on female underrepresentation has emerged in the past 20 years, however, many of those studies focus on educational access and retention and not many look at the lived experiences of women after they have left higher education. In this article, we take a different stance and explore the heterogeneous experiences of female STEM professionals in regional Australia. Through the qualitative analysis of 25 interviews, we learn what women have endured, accepted, and valued on their individual STEM journeys. While these journeys are often quite different, our interviewees independently reported having experienced similar societal prejudices and possessing similar personality traits. Our data reveals that resilience and determination proved vital for these women, as did a strong early interest in STEM. Our interviews also unearth issues in which women¿s opinions are fiercely divided, such as whether positive discrimination has been a barrier or an enabler for their careers. Based on what we have learnt from their accounts, we argue that these women have ¿survived¿ their work environments despite structural barriers, only due to their determination, resilience and fervent interest.

Estimating the shear strength of large in situ rock discontinuities is non-trivial because of the multiscale nature of roughness and the fact that only a very limited extent of discontinuity morphology is visible along traces. Recently, a novel stochastic method based on random field theory and Monte Carlo semi-analytical estimation of shear strength was proposed. The method was validated at laboratory scale and its application to one large natural surface showed that it has the potential to bypass the scale effect. However, a critical issue was reported by the authors of the study: it was found that the random field approach used could not generate the correct distribution of gradients on the simulated surfaces, which translates into an inaccurate prediction of shear strength. The authors had to manually adjust the input of the model to achieve a satisfactory prediction. This paper presents a new multiscale approach using random field theory, which now allows a rigorous generation of large synthetic 3D surfaces with controlled distribution of asperity heights and gradients from the profile of a 2D fracture trace (as might be visible in a rock face) referred to as a seed trace. Each seed trace is first decomposed into three daughter profiles corresponding to three levels of roughness. The statistics of each daughter profile form the input of the random field model at each scale level, allowing synthetic daughter surfaces to be created at each scale level. The synthetic daughter surfaces are then superimposed to obtain a composite rough surface comprising three levels of roughness. The approach was successfully validated with 25 input seed traces, coming from 5 different natural surfaces. This rigorous multiscale approach is essential to apply the stochastic method that was recently developed to predict the shear strength of large in situ discontinuities.

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.

In rockfall science, the bullet effect refers to the perforation of a rockfall mesh by a small block traveling at high speed. To date, there is still no comprehensive experimental data set investigating the underlying mechanisms of such effect. The bullet effect illustrates the fact that the capacity of a rockfall mesh depends on the size and speed of the impacting block. This paper presents the results of an experimental study on the effect of block size and mesh geometry (aperture and wire diameter) on the mesh performance. The results clearly show that the amount of energy required to perforate the mesh drops as the blocks get smaller. They also suggest that the mesh performance reaches a maximum and reduces to zero when the mesh cannot sustain the static load imposed by very large blocks. The outcome of the first series validates an analytical model for mesh perforation, making it the first simple model capturing the bullet effect. A second series of tests focused on the effect of mesh geometry and it was found that decreasing the mesh aperture by 19¿% improves the performance by 50¿% while only an extra 30¿% could be gained by increasing the wire diameter by 33¿%. The outcomes of the second series were used to discuss and redefine a dimensionless geometrical parameter G* and to validate a simple power type equation relating the mesh characteristics and the mesh performance.

From a consideration of the concepts of geological weathering and structure, it can be expected that rockfall hazards should be characteristically different in different geological environments. This paper tests this idea by looking at the geometric characteristics of rock fragments formed on natural slopes in four different geological environments in Eastern Australia, where rockfall phenomena are often characterised by rolling of pre-detached debris. By measuring the three principle dimensions and making a systematic assessment of the shape characteristics of samples of rock debris in significant geological environments, it is found that the distributions of size and shape for the surface debris are statistically different. From the results, it is shown that the size and shape of debris is directly controlled by the rock type, its weathering characteristics and the structure of the parent rock mass. The severity of rockfall hazards is shown to be relatively lower in areas of Tertiary basalt, as the size of rolling fragments is limited by closely spaced fracturing inherited from its formation and the tendency to deteriorate further as it weathers deeply and rapidly. It is also lower in areas of Palaeozoic volcanics, since these tend to produce relatively angular fragments with higher proportions of fragments that are inherently more resistant to rolling. By contrast, thickly bedded sandstones form larger blocks with a larger proportion of shapes that are more prone to rolling. The size distribution of fragments is shown to be well approximated by a log-normal statistical distribution, and using the data provided in this study, it is possible to generate the size and shape data needed to undertake a stochastic assessment of rockfall trajectories in different geological environments. © 2013.

This paper reports an analysis procedure for the evaluation of the features of the motion of blocks detaching from a steep rock wall and traveling down the slope below. Starting from the execution of real scale rock fall tests, carried out on two slopes having different morphology and lithology, the paper describes the methodology used for test interpretation and a procedure for the evaluation of the parameters best suited to the description of rock fall motion. The influence of the parameters assessed on the prediction of the rock fall trajectory was also investigated using two-dimensional and three-dimensional numerical models. These models were calibrated by means of a back analysis of the in situ tests, which also allowed the evaluation of the uncertainties involved in the parameters experimentally estimated.

Reliable prediction of shear strength for engineering scale in-situ discontinuities is still problematic. There is currently no consensus or a satisfactory method to estimate shear strength, account for surface variability and manage the effects of the recognised ¿scale effect¿ phenomenon. Shear behaviour of a discontinuity rock mass greatly depends upon the rock joint roughness, which is generally concealed within the rock mass. As such, the limited amount of accessible surface information complicates even further the exercise of shear strength prediction. Often, small size specimens (e.g. rock core) are recovered to conduct experimental tests and predictions can be made from analysing traces. In this paper, three different shear strength prediction approaches were followed and their relative performances were compared. The results reveal that the prediction of shear strength from recovered sub samples can be significantly variable. The application of Barton¿s empirical model to four selected traces, produced a large scattering of results, with the prediction highly dependent on the trace used. Conversely, the application of Casagrande and co-workers¿ stochastic approach on the same traces, produced the least scatter and provides statistical data to quantify variability and uncertainty.

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.

Numerical simulation of collapse of tunnels in weak rocks is one of the challenging problems in rock mechanics. Collapse can impose technical and financial difficulties during the construction of underground excavations. This phenomenon, particularly in shallow excavations, may lead to propagation of a caved zone toward the surface and formation of chimney subsidence. Chimney subsidence has numerous catastrophic consequences such as the adverse environmental impacts and damage to surface structures in vicinity of excavations. Due to the importance of subsidence in civil engineering, the main goal of this study is to propose an appropriate method for preventing the formation of subsidence in shallow tunnels in weak rocks. The Dolaei road tunnel is investigated as a case study, where a large chimney subsidence has occurred. The geology of the area and the geostructural characteristics of the rock mass surrounding the tunnel are briefly reviewed in this paper. Two-dimensional as well as three-dimensional numerical simulations are conducted to study the tunnel stability. Moreover, the effects of staged excavation and pre-supporting method on controlling the ground movements are studied. Outcomes of this study indicate that employing forepoling and staged excavation can be a practical remedy to control the ground movements and prevent the formation of chimney subsidence.

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.

Ground subsidence due to mining is a very complicated geomechanics phenomenon, particularly in the longwall mining method. Due to the adverse consequences of subsidence on surface structures, facilities, and environment, subsidence prediction has attracted interest in geotechnical engineering over the past few decades. Previous research has indicated that the thickness and inclination of seam, depth of working panel, and ground characteristics are the main factors which should be considered in estimating the subsidence. Methods for predicting the subsidence are: empirical, analytical and/or numerical. The geotechnical characteristics of ground affect the caving, collapse, and bulking of goaf and the extent of the cavity toward the surface. The main aim of this paper is to investigate the mechanism of ground surface subsidence based on analytical and numerical simulations. Analytical solutions are based on the theory of elasticity. The Distinct Element Method, which can simulate jointed and stratified rock mass, is used for numerical modelling of subsidence due to mining. In this paper, numerical and analytical results are obtained for two different cases of isotropic and transverse isotropic (layered) strata. Outcomes of this study indicate that the weak planes and stratification not only play a significant role in geotechnical behavior of the rock mass, but they affect the mechanism of ground movement. Results also show that considering the real structure and properties of ground will lead to more accurate estimation of measured subsidence in the field.

The paper investigates the soil-structure interaction between a clay embankment of the new high speed railway network Milan-Naples (Italy) and its clay foundation reinforced with different length driven concrete piles and geotextile. The embankment has been placed next to the old railway line. A monitoring system has been installed in a section of the embankment from the beginning of the construction of the new line. Numerical simulations describing both the old and new embankments (for a given section) have been carried out, by using Abaqus computer code, in order to back analyse the monitoring data. The analysis has been carried out in effective stress conditions, referring to long term conditions, to analyse the engineering problem to reduce the final embankment settlement. The results of the numerical simulations put into evidence that the model can reproduce the experimental data regarding the settlements below the new embankment and the interaction between the new embankment supported by piles and the old railway line embankment. © 2007 Taylor & Francis Group.