Sharif Digital Repository

Liquefaction is a phenomenon through which saturated sandy soil loses its shear strength and turns into a liquefied state. One of the most detrimental consequences of liquefaction is the reconsolidation volumetric settlements after the earthquakes, which is due to the dissipation of excess pore pressure caused by earthquakes. Severe floods can follow these settlements in free fields such as grounds close to the sea or rivers. Several researchers studied this phenomenon using data obtained from experiments in the lab or observations in the fields. Previous works were mainly based on a limited number of experimental observations and considered loadings and boundary conditions that were... 

The multiphysics process of the electromagnetic (EM) radiation induced from an embedded nanofiber scatterer subjected to SH-waves is of interest. By discarding the commonly employed electroquasistatic approximation, the fully-coupled elastodynamics and Maxwell's equations have been solved simultaneously within the mathematical framework of piezoelectric surface elasticity theory. Certain subtleties regarding the introduced interface characteristic lengths that make the examination of the size effect on the EM radiated power, scattering cross section, fundamental resonance frequency, and distribution of the magnetic field possible will be discussed. The obtained results can be helpful for... 

The coupling of rock and thermal stresses along with fluid pressure are particularly important in fractured rock masses, since stress-induced changes in permeability can be large and irreversible under perturbations resulting from various natural and induced activities. A new method is presented to model fracture permeability changes during drilling in fractured rocks. The approach includes finite element method (FEM) for fully coupled thermo-poroelastic analysis of stress distribution around borehole and displacement discontinuity method (DDM) to model fracture deformation. Three cases of overbalanced, underbalanced, and balanced drilling fluid pressure conditions are employed. The... 

The lateral movement of a liquefiable soil layer on gentle slopes is the most visible and devastating type of liquefaction-induced ground failure. Recent earthquakes have shown that this phenomenon causes severe damages to coastal structures, pier of the bridges and life-lines by exerting large lateral forces on the structures. In this paper coupled dynamic field equations of extended Biot's theory with u-p formulation are used for simulating the phenomenon and the soil behavior is modeled by a critical state two-surface plasticity model for sands. Furthermore, in this study variation of permeability coefficient during liquefaction is taken into account. The permeability coefficient is... 

Abstract Lateral Spreading, which usually occurs as a consequence of liquefaction in gently sloped loose saturated sand layers, is known to be a major source of earthquake-induced damages to structures such as quay walls, bridge piers, pipelines, and highway/railways. Therefore evaluation of the liquefaction potential and using appropriate methods for prediction of the adverse consequences of lateral spreading is of great importance. In this study, numerical modeling has been used to study lateral spreading phenomenon behind rigid waterfront structures. Coupled dynamic field equations of the extended Biot's theory with u-P formulation are used for simulating the phenomenon. A fully coupled... 

In this paper a locally mass conservative finite volume method is employed to model the one-dimensional, two-phase immiscible flow in a poroelastic media. Since, an appropriate choice of primary variables is critical in simulating multiphase subsurface flow, depending on such a choice, the governing equations can be expressed in different forms. By implementing Picard iteration to a highly nonlinear system of equations, three numerical models including pressure form, mixed form and mixed form with a modified Picard linearization are developed in this study. These models have been evaluated in terms of stability, convergence and mass conservation in various one-dimensional test cases.... 

SUMMARY: In this paper, a numerical model is developed for the fully coupled hydro-mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non-wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two-phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid... 

Airflow simulation of the whole respiratory system is still unfeasible due to the geometrical complexity of the lung airways and the diversity of the length scales involved in the problem. Even the new CT imaging system is not capable of providing accurate 3D geometries for smaller tubes, and a complete 3D simulation is impeded by the limited computational resources available. The aim of this study is to develop a fully coupled 3D-1D model to make accurate prediction of airflow and particle deposition in the whole respiratory track, with reasonable computational cost and efficiency. In the new proposed method, the respiratory tree is divided into three parts to be dealt with using different... 

Meshless methods are a relatively new type of numerical methods that have attracted the attention of many researchers over the past years. So far, a number of meshless methods have been developed and applied to solve problems in various fields of engineering, including solid mechanics and geotechnical problems. The Element-Free Galerkin (EFG) method is adopted in this study for solving the governing partial differential equations of equilibrium and continuity of pore fluid flow for numerical simulation of coupled hydro-mechanical problems. For this purpose, the weak form of the governing equations is derived by applying the weighted residual method and Galerkin technique. The penalty method... 

The soil permeability coefficient plays a key role in the process of numerical simulation of the liquefaction phenomenon. Liquefaction causes a considerable increase in soil permeability, due to the creation of easier paths for water flow. The work presented in this paper tries to investigate the effects of permeability coefficient on the results of numerical modeling of the liquefaction phenomenon. To do this, a fully coupled (u-P) formulation is employed to analyze soil displacements and pore water pressures. Two different versions of a well-calibrated critical state bounding surface plasticity model, which possesses the capability to utilize a single set of material parameters for a wide... 

Liquefaction phenomenon is usually accompanied by large amounts of settlement owing to disruption of soil structure. In addition to that, large settlement also occurs by a significant increase in soil permeability during seismic excitation. To properly simulate the post-liquefaction settlement, it is important to take the compressibility properties of the liquefied sand as well as the permeability increase into account. Using initial permeability coefficient in the course of simulation of liquefaction leads to underestimation of settlement. In addition to that, using unrealistic values for permeability may cause erroneous predictions of other aspects of soil behavior. Therefore, an accurate... 

Arbitrary Lagrangian Eulerian (ALE) method is widely used for simulation of large deformation problems, such as metal forming. However, in many such applications, modeling of the heat generation and transfer in conjunction with the stress analysis is necessary. In this work, a fully coupled dynamic ALE formulation is developed. The ALE form of energy balance equation is derived, and is coupled with the dynamic, rate dependent ALE stress analysis. The proposed formulation is used for simulation of a few thermo-mechanical problems. The effectiveness and efficiency of the ALE method is verified by comparing the results of this simulation with available experimental and numerical results  

Occurrence of liquefaction in saturated sand deposits underlying foundation of structure can cause a wide range of structural damages starting from minor settlement, and ending to general failure due to loss of bearing capacity. If the bearing capacity failure is not the problem, reliable estimation of the liquefaction-induced settlement will be of prime importance in assessment of the overall performance of the structure. Currently, there are few procedures with limited application in practice for estimation of settlement of foundations on liquefied ground. Therefore, development of a general relationship is important from the practical viewpoint. In this paper, the dynamic response of... 

A fully coupled formulation of thermo-fluid shape design problems has recently been developed in which the unknown nodal coordinates appear explicitly in the formulation of the problem. This "direct design" approach is, in principle, generally applicable and has been successfully applied in the context of potential and Euler flow models. This paper focuses on the direct design of ducts using the ideal flow model and may be considered as an addendum to the paper entitled "Direct Design of Ducts" [1]. However, a cell-vertex finite volume method is used and a different boundary condition implementation technique is applied, as compared to the method presented in the previous paper. The other... 

In this paper, formulation and implementation of finite element analysis within an Arbitrary Lagrangian Eulerian (ALE) description is presented for large deformation analysis of elasto-viscoplastic materials. The rate effects are included using a consistent procedure. An implicit algorithm with backward Euler integration scheme is used to integrate the elasto-viscoplastic constitutive equations. Also, the closed form of the consistent tangent operator is derived using the momentum balance equation to reduce the computation time. A fully coupled ALE procedure is used which includes dynamic effects. The proposed algorithm is implemented in an ALE code and its effectiveness and efficiency is... 

Mathematical simulation of non-isothermal multiphase flow in deformable unsaturated porous media is a complicated issue because of the need to employ multiple partial differential equations, the need to take into account mass and energy transfer between phases and because of the non-linear nature of the governing partial differential equations. In this paper, an analytical solution for analyzing a fully coupled problem is presented for the one-dimensional case where the coefficients of the system of equations are assumed to be constant for the entire domain. A major issue is the non-linearity of the governing equations, which is not considered in the analytical solution. In order to... 

The liquefaction phenomenon is usually accompanied by a large amount of settlement. Based on the observations made in past earthquakes, ground improvement by densification is one of the most useful approaches to reduce the liquefaction-induced settlement. Currently, there is no analytical solution for evaluation of the amount of settlement and tilting of footings that are constructed on densified ground surrounded by liquefiable soil. A number of factors, such as underlying soil properties, dimensions of the footing and earthquake loading characteristics, cause the problem to become complicated. In this paper, the dynamic response of shallow foundations on both liquefiable and... 

Numerical analyses of liquefiable sand are presented in this paper. Liquefaction phenomenon is an undrained response of saturated sandy soils when they are subjected to static or dynamic loads. A fully coupled dynamic computer code is developed to predict the liquefaction potential of a saturated sandy layer. Coupled dynamic field equations of extended Biot's theory with u-P formulation are used to determine the responses of pore fluid and soil skeleton. Generalized Newmark method is employed for integration in time. The soil behavior is modelled by two constitutive models; a critical state two-surface plasticity model, and a densification model. A class 'B' analysis of a centrifuge... 

Dynamic compaction is a widely used soil improvement method in dry and/or saturated soils. Despite its vast application, its design basis is still empirical and the mechanisms that are involved in the procedure are not fully understood. A fully coupled dynamic finite element code has been developed in order to clarify the ambiguities in the process and predict the strain/displacement field in the ground, determine depth and degree of improvement, and also calculate the pore pressure variation during the process. This model can be used as a rational design tool for dynamic compaction projects