Abstract

In this article we present numerical methods for the approximation of incompressible flows. We have addressed three problems: the stationary Stokes’ problem, the transient Stokes’ problem, and the general motion of newtonian fluids. In the three cases a discretization is employed that does not require a mesh of the domain but uses maximum entropy approximation functions. To guarantee the robustness of the solution a stabilization technique is employed. The most general problem, that of the motion of newtonian fluids, is formulated in lagrangian form. The results presented verify that stabilized meshless methods can be a competitive alternative to other approached currently in use.

Abstract

This paper presents the application of stabilized mixed explicit strain/displacement formulation (MEX-FEM) [23, 24] for solving non-linear plasticity problems in solid mechanics with strain localization. In order to use the same linear interpolation order for displacements and strains, the formulation uses the variational subscales method. Compared to the standard irreducible formulation, the proposed formulation yields improved strain and stress fields, and it is capable of addressing nearly incompressible situations. This work investigates the effects of the improved strain and stress fields in problems involving strain softening and localization leading to failure for the Mohr-Coulomb and Drucker Prager plasticity models. Numerical examples validate the ability of the proposed formulation to correctly predict failure mechanisms with localized patterns of strain, virtually free of mesh dependence and without using tracking algorithm.

Abstract

Abstract

A multidimensional extension of the HRPG method using the lowest order block finite elements is presented. First, we design a nondimensional element number that quantifies the characteristic layers which are found only in higher dimensions. This is done by matching the width of the characteristic layers to the width of the parabolic layers found for a fictitious 1D reaction–diffusion problem. The nondimensional element number is then defined using this fictitious reaction coefficient, the diffusion coefficient and an appropriate element size. Next, we introduce anisotropic element length vectors li and the stabilization parameters αi, βi are calculated along these li. Except for the modification to include the new dimensionless number that quantifies the characteristic layers, the definitions of αi, βi are a direct extension of their counterparts in 1D. Using αi, βiand li, objective characteristic tensors associated with the HRPG method are defined. The numerical artifacts across the characteristic layers are manifested as the Gibbs phenomenon. Hence, we treat them just like the artifacts formed across the parabolic layers in the reaction-dominant case. Several 2D examples are presented that support the design objective—stabilization with high-resolution

Abstract

The characteristic‐based split (CBS) stabilization procedure developed originally in fluid mechanics has been adapted successfully to solid mechanics problems. The CBS algorithm has been implemented within a finite element program using an explicit time integration scheme. Volumetric locking of linear triangular and tetrahedral elements has been successfully eliminated. The performance of the numerical algorithm is illustrated with numerical results. Comparisons with an alternative stabilization technique based on the finite calculus method also are given. Copyright © 2006 John Wiley & Sons, Ltd.

Abstract

The use of stabilization methods is becoming an increasingly well-accepted technique due to their success in dealing with numerous numerical pathologies that arise in a variety of applications in computational mechanics.

In this paper a multiscale finite element method technique to deal with pressure stabilization of nearly incompressibility problems in nonlinear solid mechanics at finite deformations is presented. A J2-flow theory plasticity model at finite deformations is considered. A mixed formulation involving pressure and displacement fields is used as starting point. Within the finite element discretization setting, continuous linear interpolation for both fields is considered. To overcome the Babuška–Brezzi stability condition, a multiscale stabilization method based on the orthogonal subgrid scale (OSGS) technique is introduced. A suitable nonlinear expression of the stabilization parameter is proposed. The main advantage of the method is the possibility of using linear triangular or tetrahedral finite elements, which are easy to generate and, therefore, very convenient for practical industrial applications.

Numerical results obtained using the OSGS stabilization technique are compared with results provided by the P1 standard Galerkin displacements linear triangular/tetrahedral element, P1/P1 standard mixed linear displacements/linear pressure triangular/tetrahedral element and Q1/P0 mixed bilinear/trilinear displacements/constant pressure quadrilateral/hexahedral element for 2D/3D nearly incompressible problems in the context of a nonlinear finite deformation J2 plasticity model.

Abstract

This paper deals with the question of strain localization associated with materials which exhibit softening due to tensile straining. A standard local isotropic Rankine damage model with strain-softening is used as exemplary constitutive model. Both the irreducible and mixed forms of the problem are examined and stability and solvability conditions are discussed. Lack of uniqueness and convergence difficulties related to the strong material nonlinearities involved are also treated. From this analysis, the issue of local discretization error in the pre-localization regime is deemed as the main difficulty to be overcome in the discrete problem. Focus is placed on low order finite elements with continuous strain and displacement fields (triangular P1P1 and quadrilateral Q1Q1), although the presented approach is very general. Numerical examples show that the resulting procedure is remarkably robust: it does not require the use of auxiliary tracking techniques and the results obtained do not suffer from spurious mesh-bias dependence.

Abstract

This paper presents the application of a stabilized mixed strain/displacement finite element formulation for the solution of nonlinear solid mechanics problems involving compressible and incompressible plasticity. The variational multiscale stabilization introduced allows the use of equal order interpolations in a consistent way. Such formulation presents two advantages when compared to the standard, displacement based, irreducible formulation: (a) it provides enhanced rate of convergence for the strain (and stress) field and (b) it is able to deal with incompressible situations. The first advantage also applies to the comparison with the mixed pressure/displacement formulation. The paper investigates the effect of the improved strain and stress fields in problems involving strain softening and localization leading to failure, using low order finite elements with continuous strain and displacement fields (P1P1 triangles or tetrahedra and Q1Q1quadrilaterals, hexahedra, and triangular prisms) in conjunction with an associative frictional Drucker–Prager plastic model. The performance of the strain/displacement formulation under compressible and nearly incompressible deformation patterns is assessed and compared to a previously proposed pressure/displacement formulation. Benchmark numerical examples show the capacity of the mixed formulation to predict correctly failure mechanisms with localized patterns of strain, virtually free from any dependence of the mesh directional bias. No auxiliary crack tracking technique is necessary.

Abstract

This paper presents an explicit mixed finite element formulation to address compressible and quasi-incompressible problems in elasticity and plasticity. This implies that the numerical solution only involves diagonal systems of equations. The formulation uses independent and equal interpolation of displacements and strains, stabilized by variational subscales. A displacement sub-scale is introduced in order to stabilize the mean-stress field. Compared to the standard irreducible formulation, the proposed mixed formulation yields improved strain and stress fields. The paper investigates the effect of this enhancement on the accuracy in problems involving strain softening and localization leading to failure, using low order finite elements with linear continuous strain and displacement fields (P1P1 triangles in 2D and tetrahedra in 3D) in conjunction with associative frictional Mohr–Coulomb and Drucker–Prager plastic models. The performance of the strain/displacement formulation under compressible and nearly incompressible deformation patterns is assessed and compared to analytical solutions for plane stress and plane strain situations. Benchmark numerical examples show the capacity of the mixed formulation to predict correctly failure mechanisms with localized patterns of strain, virtually free from any dependence of the mesh directional bias. No auxiliary crack tracking technique is necessary.

Abstract

The expression “finite calculus” refers to the derivation of the governing differential equations in mechanics by invoking balance of fluxes, forces, etc. in a space-time domain of finite size. The governing equations resulting from this approach are different from those of infinitesimal calculus theory and they incorporate new terms which depend on the dimensions of the balance domain. The new governing equations allow to derive naturally stabilized numerical schemes using any discretization procedure. The paper discusses the possibilities of the finite calculus method for the finite element solution of convection-diffusion problems with sharp gradients and incompressible fluid flow.