SUMMARYNew enrichment functions are proposed for crack modelling in orthotropic media using the extended finite element method (XFEM). In this method, Heaviside and near-tip functions are utilized in the framework of the partition of unity method for modelling discontinuities in the classical finite element method. In this procedure, by using meshless based ideas, elements containing a crack are not required to conform to crack edges. Therefore, mesh generation is directly performed ignoring the existence of any crack while the method remains capable of extending the crack without any remeshing requirement. Furthermore, the type of elements around the crack-tip remains the same as other parts of the finite element model and the number of nodes and consequently degrees of freedom are reduced considerably in comparison to the classical finite element method. Mixed-mode stress intensity factors (SIFs) are evaluated to determine the fracture properties of domain and to compare the proposed approach with other available methods. In this paper, the interaction integral (M-integral) is adopted, which is considered as one of the most accurate numerical methods for calculating stress intensity factors.

SUMMARYA novel approach based on a combination of isogeometric analysis (IGA) and extended FEM is presented for fracture analysis of structures. The extended isogeometric analysis is capable of an efficient analysis of general crack problems using nonuniform rational B-splines as basis functions for both the solution field approximation and the geometric description, and it can reproduce crack tip singular fields and discontinuity across a crack. IGA has attracted a lot of interest for solving different types of engineering problems and is now further extended for the analysis of crack stability and propagation in two-dimensional isotropic media. Concepts of the extended FEM are used in IGA to avoid the necessity of remeshing in crack propagation problems and to increase the solution accuracy around the crack tip. Crack discontinuity is represented by the Heaviside function and isotropic analytical displacement fields near a crack tip are reproduced by means of the crack tip enrichment functions. Also, the Lagrange multiplier method is used to impose essential boundary conditions. Moreover, the subtriangles technique is utilized for improving the accuracy of integration by the Gauss quadrature rule. Several two-dimensional static and quasi-static crack propagation problems are solved to demonstrate the efficiency of the proposed method and the results of mixed-mode stress intensity factors are compared with analytical and extended FEM results.

Macroscopic properties of granular materials are important in modeling a variety of flow and transport phenomena in many fields of science. Determination of these parameters has always been an issue among both researchers and engineers, mainly in view of tortuosity and permeability. This paper presents analytical functions for the tortuosity and permeability of monosized sphere arrays based on a volume averaging approach and eliminates some ambiguities by modification of the original representative elementary volume model. Veracity of the proposed formulations has been illustrated through comparisons with the latest available results on the subject. Good agreement is found.

No abstract

No abstract