Search

Screw piles have been used to support a variety of structures due to their ease of installation and high axial capacity. Recently, screw piles have been proposed as an alternative foundation solution for offshore renewable structures due to their quiet or silent installation. Due to their variable geometry, design and prediction of installation requirements and its effect on in-service capacity may be challenging. In this research study, the discrete-element method (DEM) is used to numerically recreate a series of onshore field tests. The aim of the study is to investigate the ability of DEM to be used as a practical design tool for the design and deployment of screw piles. In this case study, the effect of the geometric helix pitch on the installation torque and tensile capacity of screw piles installed into sand is investigated. DEM results show that the geometric pitch of a screw pile appears to have little effect on the installation torque. The results show that DEM has the potential to be used as a practical design procedure for complex foundation installation where the simulation needs to capture installation effects.

The final publication is available at Springer via http://dx.doi.org/10.1007/s00603-020-02142-y This paper addresses the problem of the stability of structures on calcareous rocks due to long-term weathering processes. The case study consists of a building resting on a calcarenite rock formation where two abandoned man-made caves exist directly under the structure. The boundaries of the caves were exposed to a slightly acidic environment inducing time-dependent weathering. Analyses were performed following a semi-decoupled approach, where the weathering process, driven by a reactive transport mechanism, was first solved and its results were fed to the mechanical problem which hence accounted for the spatial and temporal evolution or rock damage. For the mechanical problem, a nonlocal constitutive model was employed for the objective simulation of localised deformations. Relevant outcomes are obtained regarding the evolution of the structure’s stability and about the importance of regularising the finite element solution in the presence of brittle materials.

Biaxial shearing tests on dual-sized, 2d particle assemblies are conducted at several confining pressures. The effect of particle angularity, an important mesoscale shape descriptor, is investigated at the macro and micro levels. Macroscopically, it is observed that assemblies composed of angular particles exhibit higher strengths and dilations. The difference observed in bulk behavior due to particle angularity can be explained reasonably by considering particle-level mechanisms. A novel 2D image analysis technique is employed to estimate particle kinematics. Particle rotation results to be a key mechanism strongly influenced by particle shape determining the overall granular behavior. Unlike circular particles, angular ones are more resistant to rotations due to stronger interlocking and consequently exhibit higher strengths.

This paper aims to present a few aspects of the development of a plastic-hardening macro-element model for steel wires in flexible protection systems. First, the material behaviour is obtained using uniaxial tensile tests. Successively, the evolution of the elastic and plastic domain is obtained using a combination of physical tests, analytical models, and numerical simulations. Finally, the results obtained with the macro-element model are compared to those obtained using other approaches found in literature.

Computational load of discrete element modelling (DEM) simulations is known to increase with the number of particles. To improve the computational efficiency hybrid methods using continuous elements in the far-field, have been developed to decrease the number of discrete particles required for the model. In the present work, the performance of using such coupling methods is investigated. In particular, the coupled wall method, known as the “wall-zone” method when coupling DEM and the continuum Finite Differences Method (FDM) using the Itasca commercial codes PFC and FLAC respectively, is here analysed. To determine the accuracy and the efficiency of such a coupling approach, 3-point bending tests of cemented materials are simulated numerically. To validate the coupling accuracy first the elastic response of the beam is considered. The advantage of employing such a coupling method is then investigated by loading the beam until failure. Finally, comparing the results between DEM, DEM-FDM coupled and FDM models, the advantages and disadvantages of each method are outlined.

In this paper an isotropic hardening elastoplastic constitutive model for structured soils is applied to the simulation of a standard CPTu test in a saturated soft structured clay. To allow for the extreme deformations experienced by the soil during the penetration process, the model is formulated in a fully geometric non-linear setting, based on: i) the multiplicative decomposition of the deformation gradient into an elastic and a plastic part; and, ii) on the existence of a free energy function to define the elastic behaviour of the soil. The model is equipped with two bonding-related internal variables which provide a macroscopic description of the effects of clay structure. Suitable hardening laws are employed to describe the structure degradation associated to plastic deformations. The strain-softening associated to bond degradation usually leads to strain localization and consequent formation of shear bands, whose thickness is dependent on the characteristics of the microstructure (e.g, the average grain size). Standard local constitutive models are incapable of correctly capturing this phenomenon due to the lack of an internal length scale. To overcome this limitation, the model is framed using a non-local approach by adopting volume averaged values for the internal state variables. The size of the neighbourhood over which the averaging is performed (characteristic length) is a material constant related to the microstructure which controls the shear band thickness. This extension of the model has proven effective in regularizing the pathological mesh dependence of classical finite element solutions in the post-localization regime. The results of numerical simulations, conducted for different soil permeabilities and bond strengths, show that the model captures the development of plastic deformations induced by the advancement of the cone tip; the destructuration of the clay associated with such plastic deformations; the space and time evolution of pore water pressure as the cone tip advances. The possibility of modelling the CPTu tests in a rational and computationally efficient way opens a promising new perspective for their interpretation in geotechnical site investigations.

When testing an 1:N geotechnical structure in the centrifuge, it is desirable to choose a large scale factor (N) that can fit the small-scale model in a model container and avoid unwanted boundary effects, however, this in turn may cause scale effects when the structure is overscaled. This is more significant when it comes to small-scale modelling of sinker root-soil interaction, where root-particle size ratio is much lower. In this study the Distinct Element Method (DEM) is used to investigate this problem. The sinker root of a model root system under axial loading was analysed, with both upward and downward behaviour compared with the Finite Element Method (FEM), where the soil is modelled as a continuum in which case particle-size effects are not taken into consideration. Based on the scaling law, with the same prototype scale and particle size distribution, different scale factors/g-levels were applied to quantify effects of the ratio of root diameter (𝑑𝑟) to mean particle size (𝐷50) on the root rootsoil interaction.

Contact rolling resistance is the most widely used method to incorporate particle shape effects in the discrete element method (DEM). The main reason for this is that such approach allows for using spherical particles hence offering substantial computational benefits compared to non-spherical DEM models. This paper shows how rolling resistance parameters for 3D DEM models can be easily calibrated with 2D sand grain images.

Screw piles need to be upscaled for offshore use e.g. being an alternative foundation and anchor form for offshore floating wind turbines, although the high demand of vertical installation forces could prevent its application if conventional pitch-matched installation is used. Recent studies, using numerical and centrifuge physical tests, indicated that the vertical installation force can be reduced by adopting over-flighting which also improved axial uplift capacity of the screw pile. The current study extends the scope to axial cyclic performance with respect to the installation approach. Using quasi-static discrete element method (DEM) simulation it was found that the over-flighted screw pile showed a lower displacement accumulation rate, compared to a pitch-matched installed pile, in terms of load-controlled cyclic tests. Sensitivity analysis of the setup of the cyclic loading servo shows the maximum velocity during the tests should be limited to avoid significant exaggeration of the pile displacement accumulation but this may lead to very high run durations.

Screw piles are efficient anchors to sustain large uplift loads and can be installed with low noise or vibration. Screw piles dimensions are currently increasing, renewing research interest to reduce the installation requirements (torque and crowd or vertical force). The Discrete Element Method (DEM) is an ideal technique to investigate the complex soil behaviour during screw pile installation. Different techniques such as particle upscaling or increase of pile penetration rate have been used to reduce the CPU time to more acceptable durations (e.g. few days or weeks). This paper investigates how such techniques can affect the accuracy of the results and change the installation mechanisms. Results show that maintaining a low particle scaling factor is essential to reproduce the correct mechanism at low pile advancement ratio (AR, defined as the vertical displacement per rotation divided by the helix pitch). The pile overflighting (AR ≤ 1) creates an upwards movement of particles, which in turn creates some tension in the pile. Smaller advancement ratios require smaller particles to accurately capture this effect. Results also show that the pile penetration rate must be maintained relatively low to avoid spurious inertial effects.

Due to their low-noise installation and relatively large axial capacity, screw piles have been proposed as an alternative foundation solution in dense sand for offshore renewable energy applications in deeper water. For this to occur, significant upscaling of onshore dimensions is required. Furthermore, the effects of certain geometric features on installation requirements are still not well understood. In this work, using the three-dimensional discrete-element method, the effects of base geometry, shaft diameter and helix pitch were investigated by simulating the full installation process prior to conducting axial compression and tension tests. The results of the investigation showed it is possible to optimise the geometry of the screw pile to reduce installation requirements, in terms of both vertical installation force (up to 61%) and installation torque (up to 39%), without reducing the axial capacity of the pile significantly.

Existing guidance on the installation of screw piles suggest that they should be installed in a pitch-matched manner to avoid disturbance to the soil that may have a detrimental effect on the in-service performance of the pile. Recent insights from centrifuge modelling have shown that installing screw piles in this way requires large vertical compressive (or crowd) forces, which is inconsistent with the common assumption that screw piles pull themselves into the ground requiring minimal vertical compressive force. In this paper, through the use of the discrete element method (DEM), the effects of advancement ratio, i.e., the ratio between the vertical displacement per rotation to the geometric pitch of the helix of the screw pile helix, on the installation resistance and in-service capacity of a screw pile is investigated. The findings are further used to assess the applicability of empirical torque capacity correlation factors for large diameter screw piles. The results of the investigation show that it is possible to reduce the required vertical compressive installation force by 96% by reducing the advancement ratio and that although over-flighting a screw pile can decrease the subsequent compressive capacity, it appears to increase the tensile capacity significantly.

The virtual calibration chamber technique, based on the discrete element method, is here applied to study the standard penetration test (SPT). A macro-element approach is used to represent a rod driven with an impact like those applied to perform SPT. The rod is driven into a chamber filled with a scaled discrete analogue of a quartz sand. The contact properties of the discrete analogue are calibrated simulating two low-pressure triaxial tests. The rod is driven changing input energy and controlling initial density and confinement stress. Energy-based blowcount normalization is shown to be effective. Results obtained are in good quantitative agreement with well-accepted experimentally-based relations between blowcount, density and overburden. It is also shown that the tip resistance measured under impact dynamic penetration conditions is close to that under constant velocity conditions, hence supporting recent proposals to relate CPT and SPT results.

In granular soils grain crushing reduces dilatancy and stress obliquity enhances crushability. These are well-supported specimen-scale experimental observations. In principle, those observations should reflect some peculiar micromechanism associated with crushing, but which is it? To answer that question the nature of crushing-induced particle-scale interactions is here investigated using an efficient DEM model of crushable soil. Microstructural measures such as the mechanical coordination number and fabric are examined while performing systematic stress probing on the triaxial plane. Numerical techniques such as parallel and the newly introduced sequential probing enable clear separation of the micromechanical mechanisms associated with crushing. Particle crushing is shown to reduce fabric anisotropy during incremental loading and to slow fabric change during continuous shearing. On the other hand, increased fabric anisotropy does take more particles closer to breakage. Shear-enhanced breakage appears then to be a natural consequence of shear-enhanced fabric anisotropy. The particle crushing model employed here makes crushing dependent only on particle and contact properties, without any pre-established influence of particle connectivity. That influence does not emerge, and it is shown how particle connectivity, per se, is not a good indicator of crushing likelihood.

The Standard Penetration Test (SPT) is the most widely used method for dynamic testing of soils. The test is simple and robust but difficult to control and not fully standardized. As a result, experimental results typically show large variations and poor repeatability. To mitigate that correction factors such as energy normalization and rod length have been introduced in SPT practice. This study provides an examination of the two correction factors using models based on the discrete element method (DEM).

A 3D discrete element model is used to investigate the axial cyclic response of a small-scale displacement piles installed in Fontainebleau sand. Calibration chamber experimental results from literature are used to validate the pile penetration phase of the DEM model which is then employed to simulate stress controlled vertical cyclic loading. The crushable DEM particle model is calibrated using high pressure element test data for the same sand. The model predicts the experimental stress measurements surrounding the jacked pile in both penetrating and unloaded conditions. The DEM model is used to assess micromechanical features hard to detect using experimental and continuum numerical methods. Grain crushing within the soil is observed to occur only below the cone during pile penetration. The analysis of particle stresses and force chains highlight how arching develops around the shaft. These arching effects create a sort of shield around the shaft causing the radial stresses to be lower. After pile installation is completed, a numerical parametric study of stress controlled cyclic axial loading of the pile is performed. The results show that depending on the magnitude of the cyclic load stable or metastable pile cyclic response is attained. The cyclic load amplitude also influences in different ways both stress and density profiles around the pile. These results may serve as a step forward to the understanding of installation effects on axial cyclic performance of jacked piles in sand.

Herein a discrete-element model of a full-scale compound mesh barrier, formed by interweaved double-twist hexagonal meshes and strand rope square meshes is presented. A plastic hardening model is employed to investigate the energy dissipation terms within the barrier. The computational effort required for model initialisation is reduced through a novel generation procedure. The effect of different impact positions on the energy dissipation within the barrier and the way the impact load is transmitted to the main structural elements of the barrier are investigated through a sensitivity analysis. The results show that the impact position has both qualitative and quantitative effects on the barrier response, which significantly modify the peak load acting on the structural elements and the combination of shear, tensile and moment loading on the fence posts.

Compound mesh panels are structures in which two different nets geometries are employed: a main mesh that provides the bearing capacity and a weaker mesh with a thin sieve size to catch smaller blocks that can pass through otherwise. Typically, only the effect of the main mesh is investigated, and the weaker mesh is considered to provide negligible structural resistance. In this paper, after a calibration procedure, numerical simulations of quasi-static punch tests and a dynamic block impact on a composite double-twist and strand rope mesh are performed. The results show that, under dynamic conditions, the presence of the finer mesh lowers the peak force acting on the main mesh. This effect is not found under quasi-static conditions and has important repercussions on the overall structural resistance as the energy dissipation mechanism reduces the stress on the mesh fence posts.

This work presents a robust and mesh-independent implementation of an elasto-plastic constitutive model at large strains, appropriate for structured soils, into a Particle Finite Element code specially developed for geotechnical simulations. The constitutive response of structured soils is characterized by softening and, thus, leading to strain localization. Strain localization poses two numerical challenges: mesh dependence of the solution and computability of the solution. The former is mitigated by employing a non-local integral type regularization whereas an Implicit-Explicit integration scheme is used to enhance the computability. The good performance of these techniques is highlighted in the simulation of the cone penetration test in undrained conditions.

The state parameter (y) measures the distance from the current state to the critical state line (CSL) in the compression plane. The existence of a correlation between both the peak angle of shearing resistance (�# $ ) and peak dilatancy and y is central to many constitutive models used to predict granular soil behaviour. These correlations do not explicitly consider particle crushing. Crushing induced evolution of the particle size distribution influences the CSL position and recent research supports used of a critical state plane (CSP) to account for changes in grading. This contribution evaluates the whether the CSP can be used to calculate y and thus enable prediction of the peak angle of �# $ and peak dilatancy where crushing takes place. The data considered were generated from a validated DEM model of Fontainebleau sand that considers particle crushing. It is shown that where y is calculated by considering the CSL of the original uncrushed material there can be in a significant error in predicting the material response. Where the CSP is used there is a significant improvement in our ability to predict behaviour whether the CSP is accurately determined using a large number of tests or approximated using crushing yield envelopes. It is shown that the state parameter calculated using the previously available definition can give a false sense of security when assessing liquefaction potential of potentially crushable soils. The contribution also highlights the stress-path dependency of the relationship between �# $ and y whichever approach is used to determine y

In this study oedometric compression tests of hydrocarbon coke, Fontainebleau sand and silica sand are simulated in three dimensions using breakable particles. The method adapts a rigorous breakage criterion for elasto-brittle spheres to represent failure of grains isolated between platens or within granular masses. The breakage criterion allows for the effect of particle bulk and contact properties to be treated separately. A discrete fragmentation multigenerational approach is applied as a spawning procedure. The number of particles quickly increases during the simulation, but is kept manageable by systematic fine exclusion and upscaling. Fine exclusion leads to mass losses between generations, but that loss is accounted for outside the mechanical model. Sensitivity analysis shows that it is enough to keep 53% of the crushed particle mass within the mechanical model to correctly reproduce experimental macroscopic behaviour. Practical upscaling rules are proposed for (a) contact stiffness, (b) breakage criteria and (c) grain size distribution, and validated simulating the same test, reducing by half the initial number of particles. The results are promising as both the mechanical and grading evolution are well captured with two orders of magnitude savings in computing efficiency.

The assessment of susceptibility to failure of soft rock coastal cliffs, along with the associated failure mechanism, is not a simple task. Equilibrium conditions depend on the combination of several factors such as structural setting, rock mechanical strength, weathering processes, the hydro-mechanical action of sea waves, the variation of the rock cliff geometry, to mention some of the most important ones. From a geomechanical perspective, the brittle - strain softening behaviour of the rocks plays a key role in the onset and propagation of failure (Ciantia & Castellanza 2015). In particular, the rapid strength reduction occurring after peak under mechanical loading leading to localised deformations within shear fractures is detrimental for rock cliffs. Taking rock brittleness into account in numerical simulations under the framework of continuum mechanics is not straightforward, due to the problems related to a strong dependence of the numerical results from the adopted mesh when strain-softening laws are implemented (Vermeer and Brinkgreve 1994). Nowadays, several regularization techniques are available to control the size of the localised region and prevent the mesh dependence. Within regularization techniques, the nonlocal integral type solution has the advantage of not changing the field equations which facilitates numerical implementation. In this approach, the chosen nonlocal variables are valuated from spatial averages of the field variables in a neighbourhood, and the constitutive model is updated by replacing a local variable with its nonlocal counterpart. Consequently, the constitutive response of a Gauss point is influenced by all the integration points within a neighbourhood, with the size determined through a characteristic length (Bažant and Jirásek 2002). This contribution addresses the problem of the stability of an ideal 2-D plane strain coastal cliff, 20-m high, by means of the use of a non-local constitutive model implemented in a commercial finite element code (Mánica et al. 2018). The numerical results show insights into the evolution of the strain field and the process of slip surface/fracture propagation in the rock cliff as well as highlight the importance of regularising the finite element solution in the presence of brittle materials.

In recent years, civil engineers have started to use discrete-element modelling to simulate large-scale soil volumes thanks to technological improvements in both hardware and software. However, existing procedures to prepare ‘representative elementary volumes’ are unsatisfactory in terms of computational cost and sample homogeneity. In this work, a simple but efficient procedure to initialise large-scale discrete-element models is presented. Periodic cells are first generated with a sufficient number of particles (enough to consider the cell a representative elementary volume) matching the desired particle size distribution and equilibrated at the desired stress state, porosity and coordination number. When the cell is in equilibrium, it is replicated in space to fill the problem domain. And when the model is filled, only a small number of mechanical cycles is needed to equilibrate a large domain. The result is an equilibrated homogeneous sample at the desired initial state in a large volume.

Soft highly porous carbonate rocks, such as calcarenites, and soluble sulphate rocks, such as gypsum, are very common in the Mediterranean region and, due to their microstructure and chemical composition, are prone to water induced weathering mechanisms. Cliffs, underground cavities and other morphological features in such formations are hence affected by intense erosion phenomena and weathering processes responsible for unexpected collapses and sinkholes. Just considering the Apulian region (Italy), 150 sinkholes have been recorded since 1925, with increasing frequency since 2000 ( Fiore et al., 2018 ). The geosystem’s failure is often the short- or long-term result of a very complex hydro-chemo-mechanical process taking place at the micro-scale which can be detected and analysed by means of field and laboratory experimental test campaigns. Therefore, stability problems are often related to changes of the mechanical properties of the rock forming the cave caused by environmental weathering processes, despite the external boundary conditions are not changing with time. The paper deals with the assessment of hazard associated with the stability of abandoned underground caves, which is nowadays frequently required for land and urban planning activities. A methodological approach for hazard assessment based on a step-by-step procedure is proposed. This includes in-situ surveys, laboratory experimental studies, theoretical analyses and finally numerical investigations. The approach derives from the experience developed from several case studies analysed by the authors. In this work, two of these are presented. The first one concerns the stability of an anthropic cavity in a calcarenite formation which is affected by a water induced short-term and long-term debonding processes. The second one regards the stability of a three-level abandoned gypsum mine, the lowest level being partially flooded by water. The methodological procedure aims to evaluate the factors controlling the change of the mechanical properties of the rock so that efficient remediation measures can be designed in order to avoid any further decay of the rock mass stability with time. The proposed methodological approach, validated on real case studies, shows the convenience of performing advanced experimental, theoretical and numerical studies to properly assess the hazard in space and time and to better design the mitigation measures if they are required. The adoption of the proposed approach reduced the remediation costs of the second case study of one order of magnitude.

The goal of this paper is to show how a recently developed advanced chemo-hydro-mechanical (CHM)-coupled constitutive and numerical model for soft rocks can be applied to predict the temporal evolution of settlement damage to buildings on cavities subject to weathering. In particular, a Building Damage Index (BDI) and its evolution with time is proposed. The definition of the BDI is inspired by the work of Boscardin and Cording (J Geotech Eng 115:1–21, 1989) and uses the surface differential settlements obtained by finite element (FE) analyses to assess how far a building is from a non-acceptable service condition. By modelling the reactive transport of chemical species in 3D and using a coupled CHM constitutive and numerical model, it is possible to simulate weathering scenarios and monitor the temporal evolution of surface settlements making the BDI time dependent. This approach is applied to evaluate the damage evolution of two buildings lying on two anthropic caves in a calcarenite deposit belonging to the Calcarenite di Gravina Formation. Standard and advanced experimental tests are performed on the in situ material, and the results are used to calibrate the constitutive model. The soundness of both constitutive relationship and reactive transport solver is subsequently tested by simulating two laboratory-scale boundary value experiments. The first is a model footing test on dry and wet calcarenite, while the second is a small-scale pillar that, after the saturation-induced short-term water weakening, fails due to a long-term dissolution weathering process. Finally, both 2D and 3D coupled FE analyses simulating different weathering scenarios and corresponding settlements affecting the buildings above the considered cavities are presented. Particular attention is placed on assessing the BDI and its temporal evolution.

This paper details the work undertaken at the University of Dundee in the last 5 years to develop understanding of screw piles to allow them to be deployed offshore as an alternative foundation type to driven piles used in jacket structures. This has been prompted by several UK and European funded research initiatives to develop silent foundation techniques to mitigate affects on marine mammals and other animals. Current mitigation systems also have significant associated costs and questionable environmental credentials. Prior to starting this work, it was recognized that development would not progress unless the ability to predict installation requirements did not form an early part of the research investigation. This then led to investigation of methods to reduce installation requirements and development of new design approaches for the required new screw pile geometries. This paper details the progress of these investigations which are still ongoing.

Interaction diagrams in the generalized 3D loading space of vertical (V), horizontal (H) and moment (M) actions constitute the basis of the design of foundation structures in case of complex loads combinations. The mechanical response of such systems is frequently interpreted in terms of the ‘macroelement’ theory, where a generalized incremental constitutive relationship is introduced, linking the displacements and rotations of the foundation (playing the role of generalized strains) to the histories of applied loading components (i.e. the generalized stresses). In this paper an attempt to extend a classical macroelement framework, to the case of root-soil interaction presented. The model is calibrated on small scale experimental data on 3D printed plastic root systems, subject to combined V-H-M loads, and a parametric analysis on the main governing parameters is discussed. The comparison between numerical and experimental data suggests that the macroelement approach could be an efficient and simple analytical tool for describing the whole moment-rotation curve, overcoming the main simplifying hypotheses currently employed in arboriculture practice.

Aim (1) To understand the tree root-soil interaction under lateral and moment loading using a physical modelling technique; (2) To detect the possible factors (e.g. root architecture, water condition, and stress level) influencing a tree’s push-over behaviour; (3) To identify suitable scaling laws to use in physical modelling. Methods Two 1:20 scaled root models with different architectures (namely, deep and narrow, and shallow and wide) were reconstructed and 3D printed based on the field-surveyed root architecture data. Push-over tests were performed both in elevated-gravity (centrifuge 20-g) and normal-gravity (1-g) conditions. Results The shallow and wide model showed higher anchorage strength than the deep and narrow model. Regardless of the root architecture, the root anchorage strength measured from dry soil was higher than that from saturated soil. However, once the effective stress was the same, regardless of water conditions, the root anchorage strength would be the same. Conclusions The presence of water decreasing the soil effective stress and key lateral roots extending along the wind direction play a significant role on a tree’s push-over resistance. Centrifuge tests showed comparable results to the field pull-over measurements while 1-g model tests overestimated the root-soil interaction, which could be corrected for soil strength by using modified scaling laws.

New technologies are changing the way engineers work within the construction sector. Newly developed software solutions have provided effective methods to explore the design space at the interface between Structural Engineering and Architecture, allowing more efficient design strategies. These technologies are based on the integration of parametric generation and visualisation of geometries with powerful numerical solvers, employing user-customised routines. While the construction industry is rapidly moving the design of new construction towards a fully digitalised process, the assessment and the analysis of existing structures with such tools are still largely unexplored. In this context, a visual script for the structural assessment of out-of-plane mechanisms in historic masonry structures subject to seismic loading has recently been proposed by the authors. This relies on two successive steps of analysis, which are integrated into a digital work-flow. Datasets describing the geometric configuration of masonry structures are employed to automatically generate a non-linear Finite Element (FE) model and investigate possible collapse modes. A preliminary global analysis is performed using the commercial software ABAQUS CAE. This, in combination with the Control Surface Method (CSM), allows identifying the most likely failure mechanisms which are described by the geometry of the macro blocks. The parametric modelling of the macro-blocks geometry allows exploring the domain of possible solutions using the upper bound method of limit analysis. A Genetic Algorithms (GA) solver is used to refine the geometry of the macro-blocks and search the minimum of the upper-bound load multipliers, which guarantees equilibrium. The script is implemented in the visual programming environment offered by Rhino3D+Grasshopper. In this paper, a set of parametric analyses considering various input variables such as friction coefficient and opening incidence are performed to verify both the sensitivity and the accuracy of the proposed met, (undefined)

Gabion-type support is a favoured option to restrain bulking in pillar walls of mine footprint tunnels. It uses closely spaced short reinforcements in tunnel walls (typically fully grouted rebar) in combination with surface support (rock fragment retention systems such as shotcrete, weld wire mesh, straps, etc.). The system is installed while the rock is still mostly intact and is conceived to maintain support capacity even when, the rock attains a fully fragmented state, acting then like a gabion or earth-reinforced type retaining wall. In this paper the interaction between the support system and the highly stressed pillar walls is investigated numerically by means of finite element analyses within the framework of displacement-based design. Because the material response should capture the passage from intact rock to fully fragmented state, an advanced elasto-plastic bonded constitutive model was adopted as a simulation framework. The model is calibrated to replicate the mechanical behaviour of Bursnip Sandstone and Amarelo Pais Granite. These two rocks were selected because of high quality triaxial tests results from the literature. After showing the good performance of the model to reproduce both low and high pressure triaxial compression behaviour an extensive parametric study investigating the effects of bolt types on gabion response is presented. This work was partly supported by the Ministry of Science and Innovation of Spain through the research grant BIA2017-84752-R

The Standard Penetration Test (SPT) is the most popular example of dynamic probing, a large category of soil testing techniques. Understanding and interpretation of these tests is hampered by the difficulties of reproducing them under controlled laboratory conditions. The virtual calibration chamber technique, based on the Discrete Element Method (DEM), may supplement or substitute this complex experimentation. In this paper SPT in sand are analyzed considering the energy transfer involved. Energy balances are written for the penetrating rod and for the material in the chamber. All the terms are computed for a number of cases in which the main variables controlling test response in the field -initial density and stress level- are systematically varied. The analysis confirms previous field observations indicating that, when an energy-based interpretation is used, SPT provides a value of equivalent penetration resistance that is the same that would be obtained with a static cone penetration test. The analyses also provide an unequivocal explanation for this observation: although the impacting rod shows complicated dynamics the response of the sand is quasi-static. The first author was supported by a CSC grant. This work was partly supported by the Ministry of 486 Science and Innovation of Spain through the research grant BIA2017-84752-R.

Grouting and jet grouting are geotechnical consolidation techniques commonly employed to improve the mechanical behaviour of soils. Although these techniques are common, the micromechanical processes taking place at the local level are not yet totally understood and modelled. In this work, such a problem has been approached from a micromechanical perspective via the discrete element method by considering the local interaction among soil grains and pseudo-fluid particles. Homogeneous representative elementary volumes of a virtual analogue of silica sand have been first generated and tiny rigid frictionless particles have been subsequently injected through them, to simulate the grouting in granular materials. Various injection pressures, initial soil pressures and initial soil densities have been considered. The different diffusion patterns, the flow rate, the consequent increase in local stresses and the consequent reduction in local porosity have been discussed. To overcome the DEM computational restrictions and to speed up the injection simulations, a novel procedure based on the replication of pre-equilibrated cells has been adapted for both the initialization and injection phase. Finally, a qualitative laboratory-scale pressure grouting test has been reproduced to validate the results.

Summary The mechanical behavior of granular materials is characterized by strong nonlinearity and irreversibility. These properties have been differently described by a variety of constitutive models. To test any constitutive model, experimental data relative to the nature of the incremental stress–strain response of the material is desirable. However, this type of laboratory data is scarce because of being expensive and difficult to obtain. The discrete element method has been used several times as an alternative to obtain incremental responses of granular materials. Crushable grains add one extra source of irreversibility to granular materials. Crushability has been variously incorporated into different constitutive models. Again, it will be helpful to obtain incremental responses of crushable granular materials to test these models, but the experimental difficulties are increased. Making use of a recently introduced crushing model for discrete element simulation, this paper presents a new procedure to obtain incremental responses in discrete analogs of granular crushable materials. The parallel probe approach, previously used for uncrushable discrete analogs, is here extended to account for the presence of crushable grains. The contribution of grain crushing to the incremental irreversible strain is identified and separately measured. Robustness of the proposed method is examined in detail, paying particular attention to aspects such as dynamic instability or crushing localization. The proposed procedure is later applied to map incremental responses of a discrete analog of Fontainebleau sand on the triaxial plane. The effect of stress ratio and granular state on plastic flow characteristics is highlighted. Copyright © 2016 John Wiley & Sons, Ltd.

© Springer-Verlag Berlin Heidelberg 2016.In this work, an extended theory of plasticity with generalized hardening is proposed to describe the response of geomaterials under both mechanical and environmental processes, which include as special cases several elastoplastic constitutive equations proposed in the literature to model such processes as desaturation or suction hardening, thermal softening, chemo-mechanical coupling effects in fine-grained soils, as well as weathering of soft rocks. In the formulation of the theory, the coupling between mechanical and environmental processes takes place at two levels: first, as an additional direct contribution to the constitutive stress changes, taking place in both elastic and elastoplastic processes; and second, as a result of the evolution of the internal state variables induced by changes in the environmental process variables. This last effect is incorporated through a set of generalized hardening rules. As an example of application, the general formulation is specialized to the particular case of weak calcarenite rocks undergoing degradation processes due to plastic deformations, changes in degree of saturation (short-term debonding) and chemical dissolution of the bond material and the solid grains (long-term debonding). The resulting model is implemented in a FE code by means of an implicit generalized backward Euler algorithm, suitably modified to incorporate the full formalism of plasticity with generalized hardening. Results of numerical simulations carried out at the element level show the accuracy and efficiency properties of the proposed stresspoint algorithm. The simulation of a representative initialboundary value problem demonstrates the practical relevance of environmental degradation effects in practical applications, over periods of time comparable with the life cycle of most geotechnical structures.

Summary The mechanical properties of calcarenites are known to be significantly affected by water saturation: both stiffness and strength decrease for wetting in the short term and for chemical dissolution in the long term. Both processes mainly affect bonds among grains: immediately after inundation depositional bonds fall in suspension, whereas diagenetic bonds dissolve more slowly. In this paper, the authors started from the micro-structural analysis of the weathering processes to conceive a strain hardening hydro-chemo-mechanical coupled elastoplastic constitutive model. The concept of extended hardening rules is here enriched: weathering functions have been determined by employing a micro to macro simplified upscaling procedure. Chemical damage is incorporated into the formulation by means of a scalar damage function. Its evolution is also described by using a multiscale approach. A new term is added to the strain rate tensor in order to incorporate the dissolution induced chemical deformations developing once the soft rock is turned into a granular material. A calibration procedure for the constitutive parameters is suggested, and the model is validated by using both coupled and uncoupled chemo-mechanical experimental test results. Copyright © 2015 John Wiley & Sons, Ltd.

Chalks are highly porous rocks formed of mainly carbonate grains bonded together by calcite bridges. From a microstructural point of view they are very similar to calcarenites. The above characteristics make them prone to water induced weathering, frequently featuring large caverns and inland natural underground cavities. This study is aimed to determine the main physical processes at the base of the short - and long-term weakening experienced by these rocks when interacting with water. We present the results of microscale experimental investigations performed on calcarenites from different sites in Southern Italy. SEM, thin sections, X-ray CT observations and related analyses are used for both the interpretation - definition of the structure changes and the identification - quantification of the degradation mechanisms. Two distinct types of bonding could be identified within the rock: depositional bonds and diagenetic bonds. The diverse mechanisms linked to these two types of bonding explain both the observed fast decrease in rock strength when water fills the pores (shortterm effect of water), identified with a short-term debonding (STD), and a long-term weakening of the material, when the latter is persistently kept in water saturated conditions (long-term effect of water), identified with a long-term debonding (LTD). To highlight the microhydro-chemo mechanical processes of formation and annihilation of the depositional bonds and their role in the evolution of the mechanical strength of the material, mechanical tests on samples prepared by drying partially saturated calcarenite powder, or a mix of glass ballotini and calcarenite powder were conducted. Finally the use of transparent chemical consolidants as a more sustainable alternative to grouting is explored experimentally.

Dynamic penetration has been widely employed in some site investigation procedures and for the installation of driven piles. A 3-dimensional discrete element method model has been developed to simulate dynamic rod penetration into a calibration chamber. The chamber has been filled with a scaled analogue of Fontainebleau sand. A novel method for simulating the physical rod as boundary condition is presented. It is shown that the method leads to a good agreement in static penetration comparisons. A parametric study of involved variables shows that penetration increases with force magnitude and impact time but tip resistance is quite insensitive to them. The tests are then extended to six conditions that combine two density levels and three confining stresses. It is found that increasing the confining stress does increase tip resistance, whereas it reduces penetration distance. As expected, penetration in dense sand results in larger penetration resistances and smaller penetration distances.

A 3-dimensional discrete element method model has been developed to simulate static and dynamic rod penetration test in a calibration chamber. The chamber has been filled with a scaled analogue of Fontainebleau sand. Crushing effects on penetration results have been examined. It has been found that particle crushing reduces penetration resistance for both static and dynamic tests. Microscale observation of crushed particles has been conducted. It is shown that dynamic impact causes more crushing events. The crushed particles are distributed within 2–3 radius from the rod for both tests.

Carbonatic rocks, such as calcarenites, are very often subject to damage processes, causing a progressive degradation of their mechanical properties. In nature, in some cases, this phenomenon can cause the collapse of cliffs and underground cavities, with dangerous consequences for the anthropic environment. In this paper, the results of an experimental campaign, intended to both clarify and quantify the mechanical consequences of this process, are illustrated. To achieve such a goal, suitable physical and geotechnical indices are introduced and different time scales to describe the physical/chemical reactions induced by the water saturation of the material are taken into consideration. In particular, the authors have observed: (1) a short-term marked and instantaneous reduction in strength when water fills the pores of the rock; (2) a long-term dissolution; and (3) a progressive chemically induced reduction in the grain size. To describe the degradation processes induced by the material water saturation, owing to the complexity of the hydro-chemo-mechanical phenomena taking place within the material, suitably designed tests under controlled “weathering” conditions were also performed and discussed.

We describe the large strain implementation of an elasto-plastic model for structured soils into G-PFEM, a code developed for geotechnical simulations using the Particle Finite Element Method. The constitutive model is appropriate for naturally structured clays, cement-improved soils and soft rocks. Structure may result in brittle behavior even in contractive paths; as a result, localized failure modes are expected in most applications. To avoid the pathological mesh-dependence that may accompany strain localization, a nonlocal reformulation of the model is employed. The resulting constitutive model is incorporated into a numerical code by means of a local explicit stress integration technique. To ensure computability this is hosted within a more general Implicit-Explicit integration scheme (IMPLEX). The good performance of these techniques is illustrated by means of element tests and boundary value problems.

Long portions of the Apulian coast are steep cliffs in carbonate soft rocks. These, especially the calcarenite, are affected by weathering processes that markedly alter their mechanical properties with time, potentially leading to instability of coastal geomorphological structures. Such alterations are mainly due to chemical reactions between the solid and fluid phases, and are driven by chemical variables, which are internal variables and hence uncontrollable. In a search for the variables that drive the process of rock weakening, recourse is made to the micro scale, at which most of the chemical processes are observed and quantified. Observations using scanning electron microsope, thin sections and X-ray computed tomography analyses appear to be crucial for the understanding, interpretation and definition of the degradation mechanisms of the material. A chemo-mechanical coupled model at the meso scale of the chemically reactive stressed porous system is presented and framed in the context of a multi-scale scenario of an array of coupled phenomena. An analogous model at the macro scale is developed in parallel together with upscaling and identification procedures for meso-scale and macro-scale material constants. The main outcome of the study is a tool for predicting the progress of time-dependent weathering phenomena, potentially allowing the stability of geological structures to be assessed as it evolves with a progressing chemical degradation in a specific configuration and under a specific set of loads.

In this work acoustic emission (AE) is used as experimental evidence of the progressive nature of grain crushing. Stress controlled high pressure oedometric compression test are carried out on 1.2 mm monodisperse samples of glass beads. It was observed that the granular assembly starts to experience particle breakage at a vertical stress of about 25MPa. When this yield pressure is exceeded the glass beads start to break emitting loud impulsive sound and the vertical displacement increases rapidly. The load was increased beyond the yield stress and at each increment while the vertical stress remained constant the sample continued to emit sound. The emission of sound at a constant vertical stress indicates that crushing is a progressive failure mechanism; once the first crushing event occurs, the structure starts to rearrange causing other crushing events to occur and additional settlement. In particular, two signal processing algorithms are used on the samples of the acoustic signal to obtain two additional metrics of the crushing evolution. The first is the cumulative energy versus time. The second is the number of crushing events versus time, which is based on the automatic detection of the peaks of the sound signal envelope. There is a clear correlation between the cumulative acoustic energy emitted and the observed sample displacement. Using laser scanning, the evolution of the particle size distribution and particle shape are measured in detail so that a link between the acoustic data and the crushing intensity is established. The crushing intensity was controlled using materials with different strengths.