Search results

The purpose of this paper is to achieve numerical simulation of discontinuous rock masses.

The extended finite element method (XFEM) was used. Discontinuities (such as joints, faults, and material interfaces) are contained in the elements, thus the mesh can be generated without taking into account the existence of discontinuities. When one element contains no discontinuity, the displacement function is degenerated into that of the conventional finite element. For the element containing discontinuities, the standard displacement‐based approximation is enriched by incorporating level‐set‐based enrichment functions that model the discontinuities, and an element subdivision procedure is used to integrate the domain of the element.

Mesh generation can be simplified considerably and high‐quality meshes can be obtained. A solution with good precision can also be achieved. It is concluded that the XFEM technique is especially suitable in simulating discontinuous rock masses problems.

Crack initiation and propagation should be considered in further studies.

The paper presents a very useful numerical method for a geotechnical engineering problem that has the ability to simulate the failure process of discontinuous rock masses. The method is expected to be used widely in the deformation and stability analysis of complicated rock masses.

The paper provides a new numerical method for discontinuous rock masses that is very convenient for pre‐processing.

This paper aims to present a computational framework to generate numeric enrichment functions for two-dimensional problems dealing with single/multiple local phenomenon/phenomena. The two-scale generalized/extended finite element method (G/XFEM) approach used here is based on the solution decomposition, having global- and local-scale components. This strategy allows the use of a coarse mesh even when the problem produces complex local phenomena. For this purpose, local problems can be defined where these local phenomena are observed and are solved separately by using fine meshes. The results of the local problems are used to enrich the global one improving the approximate solution.

The implementation of the two-scale G/XFEM formulation follows the object-oriented approach presented by the authors in a previous work, where it is possible to combine different kinds of elements and analyses models with the partition of unity enrichment scheme. Beside the extension of the G/XFEM implementation to enclose the global–local strategy, the imposition of different boundary conditions is also generalized.

The generalization done for boundary conditions is very important, as the global–local approach relies on the boundary information transferring process between the two scales of the analysis. The flexibility for the numerical analysis of the proposed framework is illustrated by several examples. Different analysis models, element formulations and enrichment functions are used, and the accuracy, robustness and computational efficiency are demonstrated.

This work shows a generalize imposition of different boundary conditions for global–local G/XFEM analysis through an object-oriented implementation. This generalization is very important, as the global–local approach relies on the boundary information transferring process between the two scales of the analysis. Also, solving multiple local problems simultaneously and solving plate problems using global–local G/XFEM are other contributions of this work.

This paper proposes to extend the combination of Extended Finite Element Method (XFEM) and Level Set Method (LSM) from structural mechanics to electromagnetics. Based on this approach, the actual stage of the research work, dedicated to the investigation, development, implementation and validation of a shape optimization methodology, particularly tailored for 2D electric structures is described.

The proposed numerical approach is based on the efficiency of the XFEM and the flexibility of the LSM, to handle moving material interfaces without remeshing the whole studied domain at each optimization step.

This approach eliminates the conventional use of discrete finite elements and provides efficient, stable, accurate and faster computation schemes in comparison with other methods.

This research is limited to shape optimization of two‐dimensional electric structures, however, the work can be extended to 3D ones too.

The implementation of the proposed numerical approach for the shape optimization of a planar resistor is hereby described.

The main value of the proposed approach is a powerful and robust numerical shape optimization algorithm that demonstrates outstanding suppleness of handling topological changes, fidelity of boundary representation and a high degree of automation in comparison with other methods.

The purpose of this paper is to assess the effect of the statical admissibility of the recovered solution and the ability of the recovered solution to represent the singular solution; also the accuracy, local and global effectivity of recovery‐based error estimators for enriched finite element methods (e.g. the extended finite element method, XFEM).

The authors study the performance of two recovery techniques. The first is a recently developed superconvergent patch recovery procedure with equilibration and enrichment (SPR‐CX). The second is known as the extended moving least squares recovery (XMLS), which enriches the recovered solutions but does not enforce equilibrium constraints. Both are extended recovery techniques as the polynomial basis used in the recovery process is enriched with singular terms for a better description of the singular nature of the solution.

Numerical results comparing the convergence and the effectivity index of both techniques with those obtained without the enrichment enhancement clearly show the need for the use of extended recovery techniques in Zienkiewicz‐Zhu type error estimators for this class of problems. The results also reveal significant improvements in the effectivities yielded by statically admissible recovered solutions.

The paper shows that both extended recovery procedures and statical admissibility are key to an accurate assessment of the quality of enriched finite element approximations.

The purpose of this paper is to evaluate some numerical integration strategies used in generalized (G)/extended finite element method (XFEM) to solve linear elastic fracture mechanics problems. A range of parameters are here analyzed, evidencing how the numerical integration error and the computational efficiency are improved when particularities from these examples are properly considered.

Numerical integration strategies were implemented in an existing computational environment that provides a finite element method and G/XFEM tools. The main parameters of the analysis are considered and the performance using such strategies is compared with standard integration results.

Known numerical integration strategies suitable for fracture mechanics analysis are studied and implemented. Results from different crack configurations are presented and discussed, highlighting the necessity of alternative integration techniques for problems with singularities and/or discontinuities.

This study presents a variety of fracture mechanics examples solved by G/XFEM in which the use of standard numerical integration with Gauss quadratures results in loss of precision. It is discussed the behaviour of subdivision of elements and mapping of integration points strategies for a range of meshes and cracks geometries, also featuring distorted elements and how they affect strain energy and stress intensity factors evaluation for both strategies.

Of particular interest is the ability of the extended finite element method (XFEM) to capture transient solution and motion of phase boundaries without adaptive remeshing or moving-mesh algorithms for a physically nonlinear phase change problem. The paper aims to discuss this issue.

The XFEM is applied to solve nonlinear transient problems with a phase change. Thermal conductivity and volumetric heat capacity are assumed to be dependent on temperature. The nonlinearities in the governing equations make it necessary to employ an effective iterative approach to solve the problem. The Newton-Raphson method is used and the incremental discrete XFEM equations are derived.

The robustness and utility of the method are demonstrated on several one-dimensional benchmark problems.

The novel procedure based on the XFEM is developed to solve physically nonlinear phase change problems.

The purpose of this paper is to achieve numerical simulation of cohesive crack growth in concrete structures.

The extended finite element method (XFEM) using four-node quadrilateral element associated with the fictitious cohesive crack model is used. A mixed-mode traction-separation law is assumed for the cohesive crack in the fracture process zone (FPZ). Enrichments are considered for both partly and fully cracked elements, and it thus makes the evolution of crack to any location inside the element possible. In all. two new solution procedures based on Newton-Raphson method, which differ from the approach suggested by Zi and Belytschko (2003), are presented to solve the nonlinear system of equations. The present formulation results in a symmetric tangent matrix, conveniently in finite element implementation and programming.

The inconvenience in solving the inversion of an unsymmetrical Jacobian matrix encountered in the existing approach is avoided. Numerical results evidently confirm the accuracy of the proposed approach. It is concluded that the developed XFEM approach is especially suitable in simulating cohesive crack growth in concrete structures.

Multiple cracks and crack growth in reinforced concretes should be considered in further studies.

The research paper presents a very useful and accurate numerical method for engineering application problems that has ability to numerically simulate the cohesive crack growth of concrete structures.

The research paper provides a new numerical approach using two new solution procedures in solving nonlinear system of equations for cohesive crack growth in concrete structures that is very convenient in programming and implementation.

The purpose of this paper is to carry out a finite element simulation of a physically non-linear phase change problem in a two-dimensional space without adaptive remeshing or moving-mesh algorithms. The extended finite element method (XFEM) and the level set method (LSM) were used to capture the transient solution and motion of phase boundaries. It was crucial to consider the effects of unequal densities of the solid and liquid phases and the flow in the liquid region.

The XFEM and the LSM are applied to solve non-linear transient problems with a phase change in a two-dimensional space. The model assumes thermo-dependent properties of the material and unequal densities of the phases; it also allows for convection in the liquid phase. A non-linear system of equations is derived and a numerical solution is proposed. The Newton-Raphson method is used to solve the problem and the LSM is applied to track the interface.

The robustness and utility of the method are demonstrated on several two-dimensional benchmark problems.

The novel procedure based on the XFEM and the LSM was developed to solve physically non-linear phase change problems with unequal densities of phases in a two-dimensional space.

This paper aims to present a computationally efficient finite element model for the simulation of isothermal immiscible two-phase flow in a rigid porous media with a particular application to CO₂ sequestration in underground formations. Focus is placed on developing a numerical procedure, which is effectively mesh-independent and suitable to problems at regional scales.

The averaging theory is utilized to describe the governing equations of the involved unsaturated multiphase flow. The level-set (LS) method and the extended finite element method (XFEM) are utilized to simulate flow of the CO₂ plume. The LS is employed to trace the plume front. A streamline upwind Petrov-Galerkin method is adopted to stabilize possible occurrence of spurious oscillations due to advection. The XFEM is utilized to model the high gradient in the saturation field front, where the LS function is used for enhancing the weighting and the shape functions.

The capability of the proposed model and its features are evaluated by numerical examples, demonstrating its accuracy, stability and convergence, as well as its advantages over standard and upwind techniques. The study showed that a good combination between a mathematical model and a numerical model enables the simulation of complicated processes occurring in complicated and large geometry using minimal computational efforts.

A new computational model for two-phase flow in porous media is introduced with basic requirements for accuracy, stability, and convergence, which are met using relatively coarse meshes.

The purpose of this study is to develop an efficient numerical procedure for simulating the effect of printing orientation, as one of the primary sources of anisotropy in 3D-printed components, on their fracture properties.

The extended finite element method and the cohesive zone model (XFEM-CZM) are used to develop subroutines for fracture simulation. The ability of two prevalent models, i.e. the continuous-varying fracture properties (CVF) model and the weak plane model (WPM), and a combination of both models (WPM-CVF) are evaluated to capture fracture behavior of the additively manufactured samples. These models are based on the non-local and local forms of the anisotropic maximum tangential stress criterion. The numerical models are assessed by comparing their results with experimental outcomes of 16 different configurations of polycarbonate samples printed using the material extrusion technique.

The results demonstrate that the CVF exaggerates the level of anisotropy, and the WPM cannot detect the mild anisotropy of 3D-printed parts, while the WPM-CVF produces the best results. Additionally, the non-local scheme outperforms the local approach in terms of finite element analysis performance, such as mesh dependency, robustness, etc.

This paper provides a method for modeling anisotropic fracture in 3D-printed objects. A new damage model based on a combination of two prevalent models is offered. Moreover, the developed subroutines for implementing the non-local anisotropic fracture criterion enable a reliable crack propagation simulation in media with varying degrees of complication, such as anisotropy.