Failure mechanisms in coupled poro-plastic medium

Hadzalic, Emina,Ibrahimbegovic, Adnan,Nikolic, Mijo Techno-Press 2018 Coupled systems mechanics Vol.7 No.1

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.

3D thermo-hydro-mechanical coupled discrete beam lattice model of saturated poro-plastic medium

Hadzalic, Emina,Ibrahimbegovic, Adnan,Dolarevic, Samir Techno-Press 2020 Coupled systems mechanics Vol.9 No.2

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.

Fluid-structure interaction system predicting both internal pore pressure and outside hydrodynamic pressure

Hadzalic, Emina,Ibrahimbegovic, Adnan,Dolarevic, Samir Techno-Press 2018 Coupled systems mechanics Vol.7 No.6

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.

Failure mechanisms in coupled soil-foundation systems

Hadzalic, Emina,Ibrahimbegovic, Adnan,Dolarevic, Samir Techno-Press 2018 Coupled systems mechanics Vol.7 No.1

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.

Numerical analysis of embankment primary consolidation with porosity-dependent and strain-dependent coefficient of permeability

Balic, Anis,Hadzalic, Emina,Dolarevic, Samir Techno-Press 2022 Coupled systems mechanics Vol.11 No.2

The total embankment settlement consists of three stages: the initial settlement, the primary consolidation settlement, and the secondary consolidation settlement. The total embankment settlement is largely controlled by the primary consolidation settlement, which is usually computed with numerical models that implement Biot's theory of consolidation. The key parameter that affects the primary consolidation time is the coefficient of permeability. Due to the complex stress and strain states in the foundation soil under the embankment, to be able to predict the consolidation time more precisely, aside from porosity-dependency, the strain-dependency of the coefficient of permeability should be also taken into account in numerical analyses. In this paper, we propose a two-dimensional plane strain numerical model of embankment primary consolidation, which implements Biot's theory of consolidation with both porosity-dependent and strain-dependent coefficient of permeability. We perform several numerical simulations. First, we demonstrate the influence of the strain-dependent coefficient of permeability on the computed results. Next, we validate our numerical model by comparing computed results against in-situ measurements for two road embankments: one near the city of Saga, and the other near the city of Boston. Finally, we give our concluding remarks.

Identification of the strain-dependent coefficient of permeability by combining the results of experimental and numerical oedometer tests with free lateral movement

Balic, Anis,Hadzalic, Emina,Dolarevic, Samir Techno-Press 2022 Coupled systems mechanics Vol.11 No.1

The key parameter that affects the consolidation process of soil is the coefficient of permeability. The common assumption in the consolidation analysis is that the coefficient of permeability is porosity-dependent. However, various authors suggest that the strain-dependency of the coefficient of permeability should also be taken into account. In this paper, we present results of experimental and numerical analyses, with an aim to determine the strain-dependency of the coefficient of permeability. We present in detail both the experimental procedure and the finite element formulation of the two-dimensional axisymmetric numerical model of the oedometer test (standard and modified). We perform a set of experimental standard and modified oedometer tests. We use these experimental results to validate our numerical model and to define the model input parameter. Finally, by combining the experimental and numerical results, we propose the expression for the strain-dependent coefficient of permeability.

Soil-structure interaction effects on the seismic response of multistory frame structure

Botic, Amina,Hadzalic, Emina,Balic, Anis Techno-Press 2022 Coupled systems mechanics Vol.11 No.5

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 thicknessresting on bedrock. Amodified plane-strain numerical model isformed 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 displacementsfor the case of fixed base and for different thicknesses of clay layer are presented, compared, and discussed.

Viscoplasticity model stochastic parameter identification: Multi-scale approach and Bayesian inference

Nguyen, Cong-Uy,Hoang, Truong-Vinh,Hadzalic, Emina,Dobrilla, Simona,Matthies, Hermann G.,Ibrahimbegovic, Adnan Techno-Press 2022 Coupled systems mechanics Vol.11 No.5

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.