The stability of six regional dyke cross-sections in the Netherlands was re-assessed using the random finite element method (RFEM), which explicitly accounts for the spatial variability of strength parameters. The RFEM assessments of the cross-sections were shown to result in significantly narrower response distributions than those obtained by ignoring the spatial variability, and therefore would result in more economical designs. Given the complexity of RFEM for applications in daily engineering practice, the results obtained from the re-assessments of the six dyke cross-sections were used to propose partial factors that can be used in practice to achieve the desired reliability levels for regional dykes. When applied in a conventional semi-probabilistic assessment of a dyke cross-section, these partial factors would result in the same level of reliability as would have been obtained by carrying out an RFEM analysis of the same cross-section.

Soil liquefaction is investigated considering a saturated soil deposit and by implementing standard techniques of random field theory to distribute initial void ratio values and assess liquefaction risk. The soil domain is represented in a 2-dimensional (2D) random finite element model for the dynamic analysis of coupled behavior. Multiple Monte Carlo realizations are subjected to a base acceleration, while cyclic and small strain soil behaviours are achieved through a hypoplastic constitutive model. This investigation demonstrates that 2D stochastic simulations converge to 2D deterministic simulations when small standard deviations and/or small scales of fluctuation are used. However, large standard deviations combined with relatively large scales of fluctuation may cause significant uncertainty in the response of the soil deposit. Finally, common techniques employed to assess soil liquefaction are evaluated based on the results of the deterministic and random field analyses.

Clay rocks are multiphase porous media having a complex structure and behaviour characterised by heterogeneity, damage and viscosity, existing on a wide range of scales. The mesoscopic scale of mineral inclusions embedded in a clay matrix has an important role in the mechanisms of deformation under mechanical loading by cracking and creeping. This study introduces a micromechanical approach to model the time-dependent mechanical behaviour of clay rocks. A heterogeneous clay rock is represented at the mesoscopic scale as a composite material consisting of rigid elastic mineral inclusions (quartz, calcite and pyrite) embedded in a clay matrix. To describe the damageable rock behaviour and its failure modes at the small scale, interfaces between different mineral phases and within the clay matrix are considered. Viscous effects are incorporated inside the clay aggregates, with intergranular microfractures propagating in the clay matrix, in order to investigate their contribution to the creep behaviour of clay rock at the macroscale. The mesostructure of the clay rock is represented in digital 2D Representative Elementary Areas (REAs). The overall mesoscale behaviour of the clay rock under mechanical solicitation is numerically obtained from the REA by computational homogenisation within a two-scale finite element squared framework. Then, the model is validated at mesoscale against experimental data. The variability of the material response and the time evolution of the mineral interfacial damage state are investigated in relation to the small-scale properties and failure, while considering mesostructure variability. The results can give some valuable insights into creep behaviour of the clay rock from a small-scale perspective.

Uncertainty is inevitable in the characterisation of a geotechnical site, especially due to the inherently heterogeneous nature of the ground. In this paper, a method for characterising a subsurface with limited cone penetration test (CPT) data is proposed. The method is based on integrating predictions of CPT parameters with a probabilistic approach for subsoil classification at the CPTs. The predicted stratigraphy is able to capture the spatial variability of soil measured via CPTs and takes account of the uncertainties that arise from transforming CPT measurements into soil units as well as errors in the measurements themselves. The applicability of the proposed method is demonstrated for a site in the Netherlands. The results show that the proposed approach can identify the most likely classification in the domain with good accuracy. Furthermore, the significance of considering the uncertainties in predicting the most likely classification is illustrated via finite element stability analyses of a slope cut-out in the domain.

A kriging-based metamodelling approach for analysing the structural reliability of a sheetpile wall in a dyke is formulated. This specific problem is characterised by high target reliabilities ((Formula presented.)) in combination with a noisy and incomplete numerical model response. Starting from the original formulation of active learning kriging-based Monte Carlo simulation (AK-MCS), a robust two-stage metamodel framework is formulated in combination with adaptive multiple importance sampling, Gaussian process classification and kernel enhancements. Learning functions and convergence criteria are revised to maintain consistency with the metamodel enhancements. The developed metamodel is applied in the reliability analysis of a soil-structure interaction problem involving a sheetpile wall in a dyke body, which is representative of a class of problems encountered in engineering practice. Low dimensional example studies demonstrate the workings of the model and give insight into the model response. Full probabilistic analyses are then performed to estimate the probabilities of structural failure in a reliability updating context. The results show that after several necessary enhancements of the classical formulations, metamodelling approaches can be used successfully in combination with noisy and incomplete computational models as are often encountered in geotechnical engineering practice.

In the context of underground exploitation, the behaviour of rocks near galleries and tunnels conditions their stability. Underground drilling generates deformations, damage, fracturing, and significant modification of flow characteristics in the surrounding rock. However, the influence of small-scale characteristics and behaviour on the rock deformations and damage at engineering scale remains a complex issue. Consequently, the multiscale behaviour of a clay rock is modelled starting from the large scale of the excavation damaged zone around galleries and then enriching the approach by considering microstructural characteristics from the scale of mineral inclusions. Lastly, a double-scale numerical framework is considered. It allows to relate small- to large-scale rock behaviour in terms of deformations and material rupture. In fact, the development of damage and cracking at microscale allows to predict large-scale fracturing. The developed method focuses on a claystone in the particular context of long-term management of high-level nuclear wastes by deep geological repository. The results highlight the possibilities of double-scale computing in the prediction of the behaviour of underground engineering structures.