– In this paper, the influence of waste shear strength parameters on landfill slope stability is studied. Namely, using the limit equilibrium method in GeoStudio 2018, stability analyses of a typical landfill slope are performed using semi-probabilistic and probabilistic approaches. The semi-probabilistic computations of slope stability are performed based on the available literature recommendations for the waste shear strength parameters. The results obtained for different recommendations are compared and discussed. In the probabilistic computations, the waste shear strength parameters are treated as random variables with Gaussian random distribution, where the parameters of the distribution are again selected based on the literature recommendations. Here, the influence of dispersion in the values of waste shear strength parameters is also examined. In addition, sensitivity analyses are also performed to gain insights into the relative importance of parameters. With the aforementioned analyses, an effort was made to investigate and derive conclusions about how the selection of waste shear strength parameters affects landfill slope stability.

In this paper, we present a 3D thermo-hydro-mechanical coupled discrete beam lattice model of structure built of the nonisothermal saturated poro-plastic medium subjected to mechanical loads and nonstationary heat transfer conditions. The proposed model is based on Voronoi cell representation of the domain with cohesive links represented as inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities in axial and both transverse directions. The enhanced Timoshenko beam finite element is capable of modeling crack formation in mode I, mode II and mode III. Mode I relates to crack opening, mode II relates to in-plane crack sliding, and mode III relates to the out-of-plane shear sliding. The pore fluid flow and heat flow in the proposed model are governed by Darcy\'s law and Fourier\'s law for heat conduction, respectively. The pore pressure field and temperature field are approximated with linear tetrahedral finite elements. By exploiting nodal point quadrature rule for numerical integration on tetrahedral finite elements and duality property between Voronoi diagram and Delaunay tetrahedralization, the numerical implementation of the coupling results with additional pore pressure and temperature degrees of freedom placed at each node of a Timoshenko beam finite element. The results of several numerical simulations are presented and discussed.

This thesis studies the issue of the overall safety of structures built of heterogeneous and pore-saturated materials under extreme loads in application to fluid-structure interaction problems, such as the dam-reservoir interaction. We propose a numerical model of interaction capable of predicting main tendencies and overall behavior of pore-saturated dam structure interacting with the reservoir in failure analyses of practical interest. The proposed numerical model is first presented in two-dimensional (2D) framework and later extended to three-dimensional (3D) framework. We consider the structure built of porous cohesive material. We assume that the external fluid in interaction with the structure acts as a source of pore saturation. We model the response of the pore-saturated structure with the coupled discrete beam lattice model based on Voronoi cell representation of domain with inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities acting as cohesive links. The coupling between the solid phase and the pore fluid is handled with Biot’s porous media theory, and Darcy’s law governing the pore fluid flow. The numerical consideration of internal coupling results with an additional pressure-type degree of freedom placed at each node of the Timoshenko beam finite element, which is later used at the fluidstructure interface. The confined conditions met for external fluid placed in the reservoir enable the modeling of external fluid motion with the acoustic wave theory. For the numerical representation of the external fluid limited to small (irrotational) motion, we choose a Lagrangian formulation and the mixed displacement/pressure based finite element approximation. The end result are the displacement and pressure degrees of freedom per node of external fluid finite elements, which allows for the issue of the fluid-structure interface to be solved in an efficient and straightforward manner by directly connecting the structure and external fluid finite elements at common nodes. As a result, all computations can be performed in a fully monolithic manner. All numerical implementations and computations are performed with the research version of the computer code FEAP (Finite Element Analysis Program). The proposed numerical models of structure, external fluid and ultimately numerical model of interaction are validated in the linear elastic regime of structure response by comparing computed results against reference values obtained either with analytical solutions or continuum models. The numerical simulations in the nonlinear regime of structure response are performed with the aim to demonstrate the proposed coupled discrete beam lattice model capabilities to capture complete macro-scale response and failure mechanisms in pore-saturated structures. Finally, the proposed numerical model of interaction ability to deal with the progressive localized failure of a dam structure built of porous cohesive material under damreservoir interaction for a particular loading program was tested. To account for the temperature effects, the thermal coupling is introduced in the numerical model of the structure.

Adnan Ibrahimbegovic1, Emina Hadzalic , Emir Karavelic and Mijo Nikolic 3 1 Université de Technologie de Compiègne/Sorbonne Universités, Chair for Computational Mechanics, Centre de Recherche Royallieu, 60200 Compiègne, France, E-mail address: adnan.ibrahimbegovic@utc.fr and URL: http://roberval.utc.fr/Ibrahimbegovic-Adnan 2 Faculty of Civil Engineering, University of Sarajevo, 71000 Sarajevo, Bosnia and Herzegovina 3 Faculty of Civil Engineering, University of Split, Matice Hrvatske 15, 21000, Split, Croatia

. Behavior of soil is usually described with continuum type of failure models such as Mohr-Coulomb or Drucker-Prager model. The main advantage of these models is in a relatively simple and efficient way of predicting the main tendencies and overall behavior of soil in failure analysis of interest for engineering practice. However, the main shortcoming of these models is that they are not able to capture post-peak behavior of soil nor the corresponding failure modes under extreme loading. In this paper we will significantly improve on this state-of-the-art. In particular, we propose the use of a discrete beam lattice model to provide a sharp prediction of inelastic response and failure mechanisms in coupled soil-foundation systems. In the discrete beam lattice model used in this paper, soil is meshed with one-dimensional Timoshenko beam finite elements with embedded strong discontinuities in axial and transverse direction capable of representing crack propagation in mode I and mode II. Mode I relates to crack opening, and mode II relates to crack sliding. To take into account material heterogeneities, we determine fracture limits for each Timoshenko beam with Gaussian random distribution. We compare the results obtained using the discrete beam lattice model against those obtained using the modified three-surface elasto-plastic cap model.

The presence of the pore fluid strongly influences the reponse of the soil subjected to external loading and in many cases increases the risk of final failure. In this paper, we propose the use of a discrete beam lattice model with the aim to investigate the coupling effects of the solid and fluid phase on the response and failure mechanisms in the saturated soil. The discrete cohesive link lattice model used in this paper, is based on inelastic Timoshenko beam finite elements with enhanced kinematics in axial and transverse direction. The coupling equations for the soil-pore fluid interaction are derived from Terzaghi\'s principle of effective stresses, Biot\'s porous media theory and Darcy\'s law for fluid flow through porous media. The application of the model in soil mechanics is illustrated through several numerical simulations.

. In this paper, we present a numerical model for fluid-structure interaction between structure built of porous media and acoustic fluid, which provides both pore pressure inside porous media and hydrodynamic pressures and hydrodynamic forces exerted on the upstream face of the structure in an unified manner and simplifies fluid-structure interaction problems. The first original feature of the proposed model concerns the structure built of saturated porous medium whose response is obtained with coupled discrete beam lattice model, which is based on Voronoi cell representation with cohesive links as linear elastic Timoshenko beam finite elements. The motion of the pore fluid is governed by Darcy’s law, and the coupling between the solid phase and the pore fluid is introduced in the model through Biot’s porous media theory. The pore pressure field is discretized with CST (Constant Strain Triangle) finite elements, which coincide with Delaunay triangles. By exploiting Hammer quadrature rule for numerical integration on CST elements, and duality property between Voronoi diagram and Delaunay triangulation, the numerical implementation of the coupling results with an additional pore pressure degree of freedom placed at each node of a Timoshenko beam finite element. The second original point of the model concerns the motion of the outside fluid which is modeled with mixed displacement/pressure based formulation. The chosen finite element representations of the structure response and the outside fluid motion ensures for the structure and fluid finite elements to be connected directly at the common nodes at the fluid-structure interface, because they share both the displacement and the pressure degrees of freedom. Numerical simulations presented in this paper show an excellent agreement between the numerically obtained results and the analytical solutions.

The text analyses stresses in composite beam, as a result of time dependent concrete slab shrinkage, where in addition to analytical solution of stress distribution, numerical solution is given as well, using structural software SAP2000 v 16. For linear behavior of the headed stud anchors, whose stiffness for varying diameters is determined using separate finite element model, a comparison between numerical and analytical solutions of stress distribution is shown. Shrinkage strain is determined using EC2. For purposes of the finite element model, this strain is converted into equivalent temperature load.

. In this paper, we present the parameter identification for inelastic and multi-scale problems. First, the theoretical background of several fundamental methods used in the upscaling process is reviewed. Several key definitions including random field, Bayesian theorem, Polynomial chaos expansion (PCE), and Gauss-Markov-Kalman filter are briefly summarized. An illustrative example is given to assimilate fracture energy in a simple inelastic problem with linear hardening and softening phases. Second, the parameter identification using the Gauss-Markov-Kalman filter is employed for a multi-scale problem to identify bulk and shear moduli and other material properties in a macro-scale with the data from a micro-scale as quantities of interest (QoI). The problem can also be viewed as upscaling homogenization.

. In this paper, soil-structure interaction effects on the seismic response of multistory frame structure on raft foundation are numerically analyzed. The foundation soil profile is assumed to consists of a clay layer of variable thickness resting on bedrock. A modified plane-strain numerical model is formed in the software Plaxis, and both free vibration analysis, and earthquake analysis for a selected ground motion accelerogram are performed. The behavior of the structure is assumed to be linear elastic with Rayleigh viscous damping included. The behavior of the clay layer is modeled with a Hardening soil model with small strain stiffness. The computed results in terms of fundamental period and structural horizontal displacements for the case of fixed base and for different thicknesses of clay layer are presented, compared, and discussed