The combined finite-discrete element method (FEM/DEM) has been used to simulate processes of brittle fracturing and associated seismicity. With the newly extended FEM/DEM algorithm, two topics involving rock mechanics and geophysics are investigated. In the first topic, a velocity-weakening law is implemented to investigate the initiation of frictional slip, and an innovative method that incorporates surface roughness with varying friction coefficients is introduced to examine the influences of surface roughness. Simulated results revealed detailed responses of stresses to the propagation of the slip front. In the second topic, acoustic activities induced in confined compression tests are simulated and quantitatively studied using the internal monitoring algorithm in FEM/DEM. It is shown that with increasing confinement, AE events are spatially more concentrated and temporally more separated, accompanied by a decreasing b-value. Moreover, interesting correlation between orientations of cracks and the mechanical behavior of the rock was observed.

The combined finite-discrete element method (FEM/DEM) has been used to simulate processes of brittle fracturing and associated seismicity. With the newly extended FEM/DEM algorithm, two topics involving rock mechanics and geophysics are investigated. In the first topic, a velocity-weakening law is implemented to investigate the initiation of frictional slip, and an innovative method that incorporates surface roughness with varying friction coefficients is introduced to examine the influences of surface roughness. Simulated results revealed detailed responses of stresses to the propagation of the slip front. In the second topic, acoustic activities induced in confined compression tests are simulated and quantitatively studied using the internal monitoring algorithm in FEM/DEM. It is shown that with increasing confinement, AE events are spatially more concentrated and temporally more separated, accompanied by a decreasing b-value. Moreover, interesting correlation between orientations of cracks and the mechanical behavior of the rock was observed.

Undulatory locomotion is a common and powerful strategy used in nature at different biological scales by a broad range of living organisms, from flagellated bacteria to prehistoric snakes, which have overcome the complexity of living in ”flowable” media. By taking inspiration from this evolution-induced strategy, we aim at modelling the locomotion in a granular and viscous environment with the objective to provide more insights for designing robots for soil-like media exploration. Moreover, in contrast to common types of movement, the granular locomotion is still not well understood and is an open and challenging field. We approached this phenomenon with several tools: (i.) numerically, via coupling the Finite Element Method (FEM) with the Discrete Element Method (DEM) using ABAQUS; (ii.) analytically, by employing the Lagrangian formalism to derive the equations of motion of a discrete and continuous system subject to non-conservative forces, and (iii.) experimentally, by creating an ad-hoc set up in order to observe the migration of microfibres used for the treatment of spinal cord injuries. The computational attempts to model the motion in a granular medium involved the simulation of the dynamics of an elastic beam (FEM) surrounded by rigid spherical particles (DEM). A propulsion mechanism was introduced by sinusoidally forcing the beam’s tip normally to the longitudinal axis, while the performance of the locomotion was evaluated by means of a parametric study. Depending on the parameters of the external excitation, after a transient phase, the slender body reached a steady-state with a constant translational velocity. In order to gain physical insights, we studied a simplified version of the previous continuous beam by introducing a discrete multi-bar system. The dynamics of the latter was analytically derived, by taking into account the forces exchanged between the locomotor and the environment, according to the Resistive Force Theory. By numerically solving the equations of motion and evaluating the input energy and dissipations, we were able to define the efficiency and thus provide an effective tool to optimise the locomotion. It is worth mentioning that the two approaches, despite the different physical hypothesis, show a qualitatively and quantitatively good accordance. The numerical and analytical models previously analysed have shown promising results for the interpretation of "ad-hoc" experiments that demonstrate the migration of a microfibre embedded in a spinal cord-like matrix. This migration needs to be avoided, once the regenerative microfibre is implanted in the lesioned spinal cord, for the sake of the patients health.

Clay shales are currently being assessed as possible host rock formations for the deep geological disposal of radioactive waste. However, one main concern is that the favourable long-term isolation properties of the intact rock mass could be negatively affected by the formation of an excavation damaged zone (EDZ) around the underground openings. This thesis investigated the deformation and failure process of a clay shale, namely Opalinus Clay, with particular focus on the influence of anisotropy on the short-term response of circular tunnels. To achieve this goal, a hybrid continuum-discontinuum numerical approach was used in combination with new field measurements from the Mont Terri underground research laboratory. The response of Opalinus Clay during the excavation of a full-scale emplacement (FE) test tunnel was characterized by geodetic monitoring of wall displacements, radial extensometers and longitudinal inclinometers. The deformation measurements indicated strong directionality induced by the combined effect of in situ stress field and presence of bedding planes striking parallel to the tunnel axis, with the most severe deformation occurring in the direction approximately perpendicular to the material layering. Computer simulations were conducted using a newly-extended combined finite-discrete element method (FEM/DEM), a numerical technique which allows the explicit simulation of brittle fracturing and associated seismicity. The numerical experimentation firstly focused on the laboratory-scale analysis of failure processes (e.g., acoustic activity) in brittle rocks, and on the role of strength and modulus anisotropy in the failure behaviour of Opalinus Clay in tension and compression. The fracturing behaviour of unsupported circular excavations in laminated rock masses was then analyzed under different in situ stress conditions. Lastly, the modelling methodology was applied to the aforementioned FE tunnel to obtain original insights into the possible EDZ formation process around emplacement tunnels for nuclear waste. The calibrated numerical model suggested delamination along bedding planes and subsequent extensional fracturing as key mechanisms of the damage process potentially leading to buckling and spalling phenomena. Overall, the research findings may have a potential impact on the constructability and support design of an underground repository as well as implications for its long-term safety assessment procedure.

Clay shales are currently being assessed as possible host rock formations for the deep geological disposal of radioactive waste. However, one main concern is that the favourable long-term isolation properties of the intact rock mass could be negatively affected by the formation of an excavation damaged zone (EDZ) around the underground openings. This thesis investigated the deformation and failure process of a clay shale, namely Opalinus Clay, with particular focus on the influence of anisotropy on the short-term response of circular tunnels. To achieve this goal, a hybrid continuum-discontinuum numerical approach was used in combination with new field measurements from the Mont Terri underground research laboratory. The response of Opalinus Clay during the excavation of a full-scale emplacement (FE) test tunnel was characterized by geodetic monitoring of wall displacements, radial extensometers and longitudinal inclinometers. The deformation measurements indicated strong directionality induced by the combined effect of in situ stress field and presence of bedding planes striking parallel to the tunnel axis, with the most severe deformation occurring in the direction approximately perpendicular to the material layering. Computer simulations were conducted using a newly-extended combined finite-discrete element method (FEM/DEM), a numerical technique which allows the explicit simulation of brittle fracturing and associated seismicity. The numerical experimentation firstly focused on the laboratory-scale analysis of failure processes (e.g., acoustic activity) in brittle rocks, and on the role of strength and modulus anisotropy in the failure behaviour of Opalinus Clay in tension and compression. The fracturing behaviour of unsupported circular excavations in laminated rock masses was then analyzed under different in situ stress conditions. Lastly, the modelling methodology was applied to the aforementioned FE tunnel to obtain original insights into the possible EDZ formation process around emplacement tunnels for nuclear waste. The calibrated numerical model suggested delamination along bedding planes and subsequent extensional fracturing as key mechanisms of the damage process potentially leading to buckling and spalling phenomena. Overall, the research findings may have a potential impact on the constructability and support design of an underground repository as well as implications for its long-term safety assessment procedure.

This work proposes the computational modeling of two-dimensional media mechanical behavior with a continuous approach, related to the Finite Element Method (FEM) associated with Cohesive Fracture Mechanics (CFM) models, and a discrete approach, using the Discrete Element Method (DEM). The FEM consists in a numerical tool widely used to achieve approximate solutions of Continuum Mechanics problems, involving physical and geometrical nonlinearities phenomena with quasi-static or dynamic behaviors, having already established practical applications in many areas of science and industry. On the other hand, DEM has more recent development and has been increasingly used to model discrete nature problems involving contact, impact and fragmentation phenomena and flow of particulate systems. Focused on analysis of problems with interactions between these natures, a FEMDEM coupling code is developed to solve the problem by a sub-region scheme where the FEM is used on modeling of nucleation process and crack propagation in continuous media, and DEM is employed to model granular media, whether due its nature or its conception, in a transient behavior. The possibility of opening and propagation of cracks is considered by using CFM models, intrinsically incorporated into the FEM formulation through interfaces inserted into the inner edges of the finite element mesh. Illustrative examples are presented and discussed in order to validate the proposed formulation and implementation. / FUNDEPES - Fundação Universitária de Desenvolvimento de extensão e Pesquisa / PRH-ANP - Programa de Recursos Humanos da Agência Nacional do Petróleo / Este trabalho propõe a modelagem computacional do comportamento mecânico bidimensional de meios com abordagens contínua, relacionada ao Método dos Elementos Finitos (MEF) associado a modelos da Mecânica da Fratura Coesiva (MFC) e discreta, através do Método dos Elementos Discretos (MED). O MEF consiste em uma ferramenta numérica bastante utilizada na determinação de soluções aproximadas para problemas da Mecânica do Contínuo, envolvendo fenômenos com não linearidades físicas e geométricas associadas e com comportamento quase-estático ou dinâmico, possuindo aplicações práticas já consagradas em diversas áreas do campo científico e industrial. Por outro lado, o MED tem desenvolvimento mais recente e vem sendo cada vez mais utilizado no tratamento de problemas de natureza discreta envolvendo fenômenos de contato, impacto, fragmentação e fluxo de sistemas particulados. Com foco na análise de problemas que envolvem interações entre tais naturezas, implementa-se uma estratégia de acoplamento MEF-MED para solução do problema em subregiões, onde o MEF é utilizado na modelagem de processos de nucleação e propagação de fraturas em meios contínuos e o MED é empregado na modelagem de meios granulares por natureza, ou assim concebidos, em comportamento transiente. A possibilidade de abertura e propagação de fraturas é considerada através da utilização de modelos da MFC, incorporados intrinsecamente na formulação do MEF através de interfaces inseridas nas arestas internas da malha de elementos finitos. Exemplos ilustrativos são apresentados e discutidos visando-se validar a formulação e a implementação propostas.

A method to estimate blast-induced damage in rock considering both stress waves and gas expansion phases is presented. The method was developed by assuming a strong correlation between blast-induced damage and stress wave amplitudes, and also by adapting a 2D numerical method to estimate damage in a 3D real case. The numerical method is used to determine stress wave damage and provides an indication of zones prone to suffer greater damage by gas expansion. The specific steps carried out in this study are: i) extensive blast monitoring in hard rock at surface and underground test sites; ii) analysis of seismic waveforms in terms of amplitude and frequency and their azimuthal distribution with respect to borehole axis, iii) measurement of blast-induced damage from single-hole blasts; iv) assessment and implementation of method to utilize 2D numerical model to predict blast damage in 3D situation; v) use of experimental and numerical results to estimate relative contribution of stress waves and gas penetration to damage, and vi) monitoring and modeling of full-scale production blasts to apply developed method to estimate blast-induced damage from stress waves. The main findings in this study are: i) both P and S-waves are generated and show comparable amplitudes by blasting in boreholes; ii) amplitude and frequency of seismic waves are strongly dependent on initiation mode and direction of propagation of explosive reaction in borehole; iii) in-situ measurements indicate strongly non-symmetrical damage dependent on confinement conditions and initiation mode, and existing rock structure, and iv) gas penetration seems to be mainly responsible for damage (significant damage extension 2-4 borehole diameters from stress waves; > 22 from gas expansion). The method has the potential for application in regular production blasts for control of over-breaks and dilution in operating mines. The main areas proposed for future work are: i) verification of seismic velocity changes in rock by blast-induced damage from controlled experiments; ii) incorporation of gas expansion phase into numerical models; iii) use of 3D numerical model and verification of crack distribution prediction; iv) further studies on strain rate dependency of material strength parameters, and v) accurate measurements of in-hole pressure function considering various confinement conditions.

A method to estimate blast-induced damage in rock considering both stress waves and gas expansion phases is presented. The method was developed by assuming a strong correlation between blast-induced damage and stress wave amplitudes, and also by adapting a 2D numerical method to estimate damage in a 3D real case. The numerical method is used to determine stress wave damage and provides an indication of zones prone to suffer greater damage by gas expansion. The specific steps carried out in this study are: i) extensive blast monitoring in hard rock at surface and underground test sites; ii) analysis of seismic waveforms in terms of amplitude and frequency and their azimuthal distribution with respect to borehole axis, iii) measurement of blast-induced damage from single-hole blasts; iv) assessment and implementation of method to utilize 2D numerical model to predict blast damage in 3D situation; v) use of experimental and numerical results to estimate relative contribution of stress waves and gas penetration to damage, and vi) monitoring and modeling of full-scale production blasts to apply developed method to estimate blast-induced damage from stress waves. The main findings in this study are: i) both P and S-waves are generated and show comparable amplitudes by blasting in boreholes; ii) amplitude and frequency of seismic waves are strongly dependent on initiation mode and direction of propagation of explosive reaction in borehole; iii) in-situ measurements indicate strongly non-symmetrical damage dependent on confinement conditions and initiation mode, and existing rock structure, and iv) gas penetration seems to be mainly responsible for damage (significant damage extension 2-4 borehole diameters from stress waves; > 22 from gas expansion). The method has the potential for application in regular production blasts for control of over-breaks and dilution in operating mines. The main areas proposed for future work are: i) verification of seismic velocity changes in rock by blast-induced damage from controlled experiments; ii) incorporation of gas expansion phase into numerical models; iii) use of 3D numerical model and verification of crack distribution prediction; iv) further studies on strain rate dependency of material strength parameters, and v) accurate measurements of in-hole pressure function considering various confinement conditions.