Search results

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.

The purpose of this paper is devoted to present an accurate assessment for determine natural frequencies for uniform and non-uniform Euler-Bernoulli beams and frames by an adaptive generalized finite element method (GFEM). The present paper concentrates on developing the C1 element of the adaptive GFEM for vibration analysis of Euler-Bernoulli beams and frames.

The variational problem of free vibration is formulated and the main aspects of the adaptive GFEM are presented and discussed. The efficiency and convergence of the proposed method in vibration analysis of uniform and non-uniform Euler-Bernoulli beams are checked. The application of this technique in a frame is also presented.

The present paper concentrates on developing the C1 element of the adaptive GFEM for vibration analysis of Euler-Bernoulli beams and frames. The GFEM, which was conceived on the basis of the partition of unity method, allows the inclusion of enrichment functions that contain a priori knowledge about the fundamental solution of the governing differential equation. The proposed enrichment functions are dependent on the geometric and mechanical properties of the element. This approach converges very fast and is able to approximate the frequency related to any vibration mode.

The main contribution of the present study consisted in proposing an adaptive GFEM for vibration analysis of Euler-Bernoulli uniform and non-uniform beams and frames. The GFEM results were compared with those obtained by the h and p-versions of FEM and the c-version of the CEM. The adaptive GFEM has shown to be efficient in the vibration analysis of beams and has indicated that it can be applied even for a coarse discretization scheme in complex practical problems.

This paper aims to propose a generalized finite element technique that can accurately approximate the solution of the flexural-shear cantilever model of wall-frame structures proposed by Heidebrecht and Stafford Smith.

This approach adopts scaled monomials as enrichment functions, and they are highly effective in accurately capturing the solution of the problem, as it consists of smooth functions such as polynomials, hyperbolic and trigonometric functions. Several numerical experiments are performed on the static and modal analyses of the flexural-shear cantilever wall-frame structures using the proposed generalized finite element method (GFEM), and their accuracies are compared with those obtained using the standard finite element method.

The proposed GFEM is able to achieve theoretical convergence rates of the static and modal analyses, which are, in principle, identical to those of the standard FEM, for various polynomial orders of its shape functions such as quadratic, cubic and quartic orders. The proposed GFEM with quartic enrichment functions can provide more accurate solutions than the standard FEM, and thus can be effectively used at the initial design stage of high-rise wall-frame structures.

This work is the first paper where the GFEM is applied to the analysis of high-rise wall-frame structures, and the developed technique can be used as a good analysis tool at the initial design stage.

The purpose of this paper is to check the efficiency of isogeometric analysis (IGA) by comparing its results with classical finite element method (FEM), generalized finite element method (GFEM) and other enriched versions of FEM through numerical examples of free vibration problems.

Since its conception, IGA was widely applied in several problems. In this paper, IGA is applied for free vibration of elastic rods, beams and trusses. The results are compared with FEM, GFEM and the enriched methods, concerning frequency spectra and convergence rates.

The results show advantages of IGA over FEM and GFEM in the frequency error spectra, mostly in the higher frequencies.

Isogeometric analysis shows a feasible tool in structural analysis, with emphasis for problems that requires a high amount of vibration modes.

The purpose of this paper is to evaluate the capabilities of the generalized finite element method (GFEM) under the context of the geometrically nonlinear analysis. The effect of large displacements and deformations, typical of such analysis, induces a significant distortion of the element mesh, penalizing the quality of the standard finite element method approximation. The main concern here is to identify how the enrichment strategy from GFEM, that usually makes this method less susceptible to the mesh distortion, may be used under the total and updated Lagrangian formulations.

An existing computational environment that allows linear and nonlinear analysis, has been used to implement the analysis with geometric nonlinearity by GFEM, using different polynomial enrichments.

The geometrically nonlinear analysis using total and updated Lagrangian formulations are considered in GFEM. Classical problems are numerically simulated and the accuracy and robustness of the GFEM are highlighted.

This study shows a novel study about GFEM analysis using a complete polynomial space to enrich the approximation of the geometrically nonlinear analysis adopting the total and updated Lagrangian formulations. This strategy guarantees the good precision of the analysis for higher level of mesh distortion in the case of the total Lagrangian formulation. On the other hand, in the updated Lagrangian approach, the need of updating the degrees of freedom during the incremental and iterative solution are for the first time identified and discussed here.

This paper presents a novel framework for solving elliptic partial differential equations (PDEs) over irregularly spaced meshes on bounded domains.

Second‐generation wavelet construction gives rise to a powerful generalization of the traditional hierarchical basis (HB) finite element method (FEM). A framework based on piecewise polynomial Lagrangian multiwavelets is used to generate customized multiresolution bases that have not only HB properties but also additional qualities.

For the 1D Poisson problem, we propose – for any given order of approximation – a compact closed‐form wavelet basis that block‐diagonalizes the stiffness matrix. With this wavelet choice, all coupling between the coarse scale and detail scales in the matrix is eliminated. In contrast, traditional higher‐order (n>1) HB do not exhibit this property. We also achieve full scale‐decoupling for the 2D Poisson problem on an irregular mesh. No traditional HB has this quality in 2D.

Similar techniques may be applied to scale‐decouple the multiresolution finite element (FE) matrices associated with more general elliptic PDEs.

By decoupling scales in the FE matrix, the wavelet formulation lends itself particularly well to adaptive refinement schemes.

The paper explains second‐generation wavelet construction in a Lagrangian FE context. For 1D higher‐order and 2D first‐order bases, we propose a particular choice of wavelet, customized to the Poisson problem. The approach generalizes to other elliptic PDE problems.

It has been usual to prefer an enrichment pattern independent of the mesh when applying singular functions in the Generalized/eXtended finite element method (G/XFEM). This choice, when modeling crack tip singularities through extrinsic enrichment, has been understood as the only way to surpass the typical poor convergence rate obtained with the finite element method (FEM), on uniform or quasi-uniform meshes conforming to the crack. Then, the purpose of this study is to revisit the topological enrichment strategy in the light of a higher-order continuity obtained with a smooth partition of unity (PoU). Aiming to verify the smoothness' impacts on the blending phenomenon, a series of numerical experiments is conceived to compare the two GFEM versions: the conventional one, based on piecewise continuous PoU's, and another which considers PoU's with high-regularity.

The stress approximations right at the crack tip vicinity are qualified by focusing on crack severity parameters. For this purpose, the material forces method originated from the configurational mechanics is employed. Some attempts to improve solution using different polynomial enrichment schemes, besides the singular one, are discussed aiming to verify the transition/blending effects. A classical two-dimensional problem of the linear elastic fracture mechanics (LEFM) is solved, considering the pure mode I and the mixed-mode loadings.

The results reveal that, in the presence of smooth PoU's, the topological enrichment can still be considered as a suitable strategy for extrinsic enrichment. First, because such an enrichment pattern still can treat the crack independently of the mesh and deliver some advantage in terms of convergence rates, under certain conditions, when compared to the conventional FEM. Second, because the topological pattern demands fewer degrees of freedom and impacts conditioning less than the geometrical strategy.

Several outputs are presented, considering estimations for the J–integral and the angle of probable crack advance, this last computed from two different strategies to monitoring blending/transition effects, besides some comments about conditioning. Both h- and p-behaviors are displayed to allow a discussion from different points of view concerning the topological enrichment in smooth GFEM.

In this paper, the polygonal-FEM technique is presented in modeling large deformation – large sliding contact on non-conformal meshes. The purpose of this paper is to present a new technique in modeling arbitrary interfaces and discontinuities for non-linear contact problems by capturing discontinuous deformations in elements cut by the contact surface in uniform non-conformal meshes.

The geometry of contact surface is used to produce various polygonal elements at the intersection of the interface with the regular FE mesh, in which the extra degrees-of-freedom are defined along the interface. The contact constraints are imposed between polygonal elements produced along the contact surface through the node-to-surface contact algorithm.

Numerical convergence analysis is carried out to study the convergence rate for various polygonal interpolation functions, including the Wachspress interpolation functions, the metric shape functions, the natural neighbor-based shape functions, and the mean value shape functions. Finally, numerical examples are solved to demonstrate the efficiency of proposed technique in modeling contact problems in large deformations.

A new technique is presented based on the polygonal-FEM technique in modeling arbitrary interfaces and discontinuities for non-linear contact problems by capturing discontinuous deformations in elements cut by the contact surface in uniform non-conformal meshes.

This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.

The purpose of this paper is to extend the generalized finite‐difference calculus of flexible local approximation methods (FLAME) to problems where local analytical solutions are unavailable.

FLAME uses accurate local approximations of the solution to generate difference schemes with small consistency errors. When local analytical approximations are too complicated, semi‐analytical or numerical ones can be used instead. In the paper, this strategy is applied to electrostatic multi‐particle simulations and to electromagnetic wave propagation and scattering. The FLAME basis is constructed by solving small local finite‐element problems or, alternatively, by a local multipole‐multicenter expansion. As yet another alternative, adaptive FLAME is applied to problems of wave propagation in electromagnetic (photonic) crystals.

Numerical examples demonstrate the high rate of convergence of new five‐ and nine‐point schemes in 2D and seven‐ and 19‐point schemes in 3D. The accuracy of FLAME is much higher than that of the standard FD scheme. This paves the way for solving problems with a large number of particles on relatively coarse grids. FLAME with numerical bases has particular advantages for the multi‐particle model of a random or quasi‐random medium.

Irregular stencils produced by local refinement may adversely affect the accuracy. This drawback could be rectified by least squares FLAME, where the number of stencil nodes can be much greater than the number of basis functions, making the method more robust and less sensitive to the irregularities of the stencils.

Previous applications of FLAME were limited to purely analytical basis functions. The present paper shows that numerical bases can be successfully used in FLAME when analytical ones are not available.