Sharif Digital Repository

In this paper, the 3D cohesive crack propagation is presented in partially saturated porous media. The double-nodded zero-thickness cohesive interface elements are employed to capture the mixed mode fracture behavior. In order to describe the behavior of fractured media, two balance equations are applied similar to those employed for the mixture of solid-fluid phase in semi-saturated media, including: the momentum balance of fractured media, and the balance of fluid mass within the fracture. Crack permeability is modified based on the data obtained from experimental results to consider the roughness of fracture walls effect  

An enriched-FEM technique is presented for the crack growth simulation in large deformation ductile fracture problems using a non-local damage-plasticity model in the framework of eXtended Finite Element Method (X-FEM). The Lemaitre damage-plasticity model is used to capture the material degradation effect, in which the non-locality is enforced by solving a Helmholtz type equation in combination with the governing equation of the system based on an operator-split technique. A convergence study is performed to investigate the performance of X-FEM technique in plasticity problems. The accuracy and effectiveness of proposed X-FEM damage-plasticity model are verified through several numerical... 

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... 

In this paper, a polygonal-FEM technique is presented in modeling of arbitrary interfaces in large deformations. The method is used to model the internal interfaces and arbitrary geometries using a uniform non-conformal mesh. The technique is applied to capture discontinuous deformations in the non-conformal elements, which are cut by the interface in a uniform regular mesh. In this approach, a uniform non-conformal mesh is decomposed into sub-elements that conform to the internal interfaces. The geometry of interface is used to produce various triangular, quadrilateral and pentagonal elements at the intersection of interface with regular FE mesh, in which the extra degrees-of-freedom are... 

In the present paper, a fully coupled numerical model is developed for the hydromechanical analysis of deforming, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non-wetting pore fluids. The governing equations involving the coupled two-phase fluid flow and deformation processes in partially saturated porous media containing cohesive cracks are derived within the framework of the generalized Biot theory. The displacement of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the three-phase formulation. A softening cohesive law is employed to describe the nonlinear... 

In this paper, the continuum damage mechanics model originally proposed by Lemaitre [1] is presented through an adaptive finite element method for three-dimensional ductile materials. The macro-crack initiation-propagation criterion is used based on the distribution of damage variable in the continuum damage model. The micro-crack closure effect is incorporated to simulate the damage evolution more realistic. The Zienkiewicz-Zhu posteriori error estimator is employed in conjunction with a weighted superconvergence patch recovery (SPR) technique at each patch to improve the accuracy of error estimation and data transfer process. Finally, the robustness and accuracy of proposed computational... 

In this paper, a generalized-FEM technique is presented in modeling of arbitrary interfaces in large deformations. The method is used to model the internal interfaces and arbitrary geometries using a uniform non-conformal mesh. The technique is applied to capture independent deformations at both sides of separated element cut by the interface in a uniform regular mesh. In this approach, a uniform non-conformal mesh is decomposed into subelements that conform to the internal interfaces. The geometry of interface is used to produce various triangular, quadrilateral and pentagonal elements at the intersection of interface with regular FE mesh, in which the extra degrees-of-freedom are defined... 

In this paper, an automated adaptive remeshing procedure is presented for simulation of arbitrary shape crack growth in multiple cracked bodies. The Zienkiewicz-Zhu error estimator is employed in conjunction with the modified superconvergent patch recovery (SPR) technique to obtain more accurate estimation of error. A stability analysis is performed to determine active cracks from a set of competitive cracks. Various crack growth criteria together with the respective crack trajectory prediction are compared. Several numerical examples are illustrated to demonstrate the efficiency, robustness and accuracy of computational algorithm in the simulation of multiple crack growth  

In this paper, the dynamic large deformation X-FEM method is presented for modeling the full process of dynamic ductile fracture based on a nonlocal damage visco-plasticity model. The effect of inertia is modeled using an explicit central difference scheme which is enhanced through the use of mass lumping, reduced integration with hourglass control, and numerical damping. The material nonlinearity and the flow stress dependency on strain rate, hardening and temperature are modeled with the Johnson-Cook visco-plastic model. The micro-void nucleation, growth and coalescence are modeled macroscopically with an isotropic damage model. The localization phenomenon due to the damage and thermal... 

In this paper, the continuum damage mechanics model originally proposed by Lemaitre (Journal of Engineering Materials and Technology. 1985; 107: 83-89) is presented through an adaptive finite element method for three-dimensional ductile materials. The macro-crack initiation-propagation criterion is used based on the distribution of damage variable in the continuum damage model. The microcrack closure effect is incorporated to simulate the damage evolution more realistic. The Zienkiewicz-Zhu posteriori error estimator is employed in conjunction with a weighted Superconvergence Patch Recovery (SPR) technique at each patch to improve the accuracy of error estimation and data transfer process.... 

In this paper, an analysis of kinematics of the isotropic elastic Cosserat continuum is presented in infinitesimal and finite deformations. Emphasis is given on the applicability of corotational stress rates for hypoelasticity in micro-polar continua. A non-linear finite element analysis is performed with an explicit formulation of tangent stiffness matrices in the case of Truesdell stress and couple stress rates. A comprehensive path-dependent procedure is employed based on the arc-length method to calculate the stability points and handle the snap-back problem. Finally, the accuracy and efficiency of method are illustrated by numerical examples  

In this paper, a new computational technique is presented for the modeling of moving boundaries in large plastic deformations based on an enriched arbitrary Lagrangian-Eulerian finite element method. An Arbitrary Lagrangian-Eulerian (ALE) technique is employed to capture the advantages of both Lagrangian and Eulerian methods and alleviate the drawbacks of mesh distortion in Lagrangian formulation. An enriched finite element method is implemented based on the extended FEM technique to capture the arbitrary interfaces independent of element boundaries. The process is accomplished by performing a splitting operator to separate the material (Lagrangian) phase from the convective (Eulerian)... 

In this paper, a fully coupled numerical model is developed for the modeling of the hydraulic fracture propagation in porous media using the extended finite element method in conjunction with the cohesive crack model. The governing equations, which account for the coupling between various physical phenomena, are derived within the framework of the generalized Biot theory. The fluid flow within the fracture is modeled using the Darcy law, in which the fracture permeability is assumed according to the well-known cubic law. By taking the advantage of the cohesive crack model, the nonlinear fracture processes developing along the fracture process zone are simulated. The spatial discretization... 

In this paper, a numerical model is developed for the fully coupled analysis of deforming porous media containing weak discontinuities which interact with the flow of two immiscible, compressible wetting and non-wetting pore fluids. The governing equations involving the coupled solid skeleton deformation and two-phase fluid flow in partially saturated porous media are derived within the framework of the generalized Biot theory. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three-phase formulation. The other variables are incorporated into the model via the experimentally determined functions that specify the... 

This paper presents a general framework for modeling dynamic large deformation contact-impact problems based on the phantom-node extended finite element method. The large sliding penalty contact formulation is presented based on a master-slave approach which is implemented within the phantom-node X-FEM and an explicit central difference scheme is used to model the inertial effects. The method is compared with conventional contact X-FEM; advantages, limitations and implementational aspects are also addressed. Several numerical examples are presented to show the robustness and accuracy of the proposed method. © 2017 Springer-Verlag GmbH Germany  

The main goal of the current study is developing an advanced and robust numerical tool for accurate capturing heat front propagation. In some applications such as impermeable medium, Heat transfer in the surrounding domain of fracture acts just as a conduction process but the heat transfer through the fractures appears as a convection process. From a mathematical point of view, a parabolic partial differential equation (PDE) should be solved in the surrounding domain whereas a hyperbolic PDE should be solved in the domain of fractures. In fact, they have completely different treatments and this is one of the complicated problems in this area. In this paper, the presence of fractures and... 

This paper presents a general framework for modeling dynamic large deformation contact-impact problems based on the phantom-node extended finite element method. The large sliding penalty contact formulation is presented based on a master-slave approach which is implemented within the phantom-node X-FEM and an explicit central difference scheme is used to model the inertial effects. The method is compared with conventional contact X-FEM; advantages, limitations and implementational aspects are also addressed. Several numerical examples are presented to show the robustness and accuracy of the proposed method. © 2017, Springer-Verlag GmbH Germany  

In this paper, a multiscale homogenization approach is developed for fully coupled saturated porous media to represent the idealized sugar cube model, which is generally employed in fractured porous media on the basis of dual porosity models. In this manner, an extended version of the Hill-Mandel theory that incorporates the microdynamic effects into the multiscale analysis is presented, and the concept of the deformable dual porosity model is demonstrated. Numerical simulations are performed employing the multiscale analysis and dual porosity model, and the results are compared with the direct numerical simulation through 2 numerical examples. Finally, a combined multiscale-dual porosity... 

In this paper, an adaptive finite element analysis is presented for 3D modeling of non-planar curved crack growth. The fracture mechanical evaluation is performed based on a general technique for non-planar curved cracks. The Schollmann's crack kinking criterion is used for the process of crack propagation in 3D problems. The Zienkiewicz-Zhu error estimator is employed in conjunction with a weighted SPR technique at each patch to improve the accuracy of error estimation. Applying the proposed technique to 3D non-planar curved crack growth problems shows significant improvements particularly at the boundaries and near crack tip regions. Several numerical examples are presented to illustrate... 

In this paper, a multiscale homogenization approach is developed for fully coupled saturated porous media to represent the idealized sugar cube model, which is generally employed in fractured porous media on the basis of dual porosity models. In this manner, an extended version of the Hill-Mandel theory that incorporates the microdynamic effects into the multiscale analysis is presented, and the concept of the deformable dual porosity model is demonstrated. Numerical simulations are performed employing the multiscale analysis and dual porosity model, and the results are compared with the direct numerical simulation through 2 numerical examples. Finally, a combined multiscale-dual porosity...