Numerical analysis of anisotropic stiffness and strength for geomaterials

In numerical modelling,selection of the constitutive model is a critical factor in predicting the actual response of a geomaterial.The use of oversimplified or inadequate models may not be sufficient to reproduce the actual geomaterial behaviour.That selection is especially relevant in the case of aniso-tropic rocks,and particularly for shales and slates,whose behaviour may be affected,e.g.well stability in geothermal or oil and gas production operations.In this paper,an alternative anisotropic constitutive model has been implemented in the finite element method software CODE_BRIGHT,which is able to account for the anisotropy of shales and slates in terms of both deformability and strength.For this purpose,a transversely isotropic version of the generalised Hooke’s law is adopted to represent the stiffness anisotropy,while a nonuniform scaling of the stress tensor is introduced in the plastic model to represent the strength anisotropy.Furthermore,a detailed approach has been proposed to determine the model parameters based on the stressestrain results of laboratory tests.Moreover,numerical analyses are performed to model uniaxial and triaxial tests on Vaca Muerta shale,Bossier shale and slate from the northwest of Spain(NW Spain slate).The experimental data have been recovered from the literature in the case of the shale and,in the case of the slate,performed by the authors in terms of stress-strain curves and strengths.A good agreement can be generally observed between numerical and experi-mental results,hence showing the potential applicability of the approach to actual case studies.Therefore,the presented constitutive model may be a promising approach for analysing the anisotropic behaviour of rocks and its impact on well stability or other relevant geomechanical problems in aniso-tropic rocks.

An Educational GUI-Based Software for Dynamic Analysis of Framed Structural Models

This paper introduces a new version of the open-source educational software, LESM (Linear Elements Structure Model), developed in MATLAB for structural analysis of one-dimensional models such as frames, trusses, and grillages. The updated program includes dynamic analysis, which incorporates inertial and damping effects, time-dependent load conditions, and a transient solver with multiple time integration schemes. The software assumes small displacements and linear-elastic material behavior. The paper briefly explains the theoretical basis for these developments and the reorganization of the source code using Object-Oriented Programming (OOP). The updated Graphical User Interface (GUI) allows interactive use of dynamic analysis features and displays new results such as animations, envelope diagrams of internal forces, phase portraits, and the response of degrees-of-freedom in time and frequency domain. The new version was used in a structural dynamics course, and new assignments were elaborated to improve students’ understanding of the subject.

A block particle coupled model and its application to landslides

To simulate the progressive failure of slope,a block particle coupled model is introduced.Particle oriented cell mapping(POCM)algorithm is used to enhance the search efficiency,and particlepoint,particle-edge,particle-face contact detecting method is adopted to establish contact pair between particles and blocks precisely.Strain softening Mohr Coulomb model with tensile cutoff is adopted for blocks,and brittle Mohr Coulomb model is used for particles.The particle-block replacement approach is used to describe the fracture and fragmentation process of continuum media.Once the cohesion or tensile strength of one block reaches zero,the block will be deleted,and particles are generated at the same place with all information inherited from the deleted block.Some numerical cases related to landslides demonstrate the precision and rationality of the coupled model.