COLLECTIONS

Documents published in Scipedia

• Stabilization techniques for finite element analysis of convection-diffusion problems

  E. Oñate, M. Manzan

  (2000). Research Report, Nº PI183

  
  

  Abstract

  The accurate solution of convection type problems on practical grids has been ever a challenging issue, and invariably some sort of stabilization is needed in order to get a physical solution. This has pushed researchers to develop various stabilization algorithms used in every day practice by numerical analysts. In this chapter some methods are presented along with a new finite increment calculus approach to obtain the different algorithms using higher order conservation equations.

  Abstract
  The accurate solution of convection type problems on practical grids has been ever a challenging issue, and invariably some sort of stabilization is needed in order to get [...]

  • 26
  •  read

• The Monte Carlo method. Application to the stochastic analysis of sheet stamping processes

  J. Hurtado, F. Zárate, E. Oñate

  (2000). Research Report, Nº PI184

  
  

  Abstract

  Stochastic Mechanics is a rapidly growing area of research, whose importance is being recognized not only in academic circles but also in industrial practice. This is no doubt due to the fact that most structural properties and loads are either random or uncertain. The first term refers to a natural chaotic variation of the parameter, while the second is associated to the human lack of knowledge about it. Both kinds of unpredictability work together in rendering doubtful the results of a (usually single) deterministic mechanical analysis. When thinking about the randomness and uncertainty linked to all physical parameters and phenomena a big question mark closes the large list of numbers produced by a finite element calculation. In Stochastic Mechanics there are several techniques to analyse the natural scatter of strains and stresses caused by the dispersion in the given loads and/or the structural parameters. The most general one is the Monte Carlo method. However, it must be recognized that is as well the most costly in computational terms. Nevertheless, this cost has becoming feasible with the advance in Computer Science, specially with the advent of parallel computing, due to the fact that a Monte Carlo calculation is intrinsically a task that can be performed in parallel. The present report is intended to provide the reader an introduction to the Monte Carlo method in the context of Computational Mechanics. The technique is applied to the analysis of the uncertainty spread in a stamping process. The first chapter summarises the Monte Carlo method and its theoretical backgrounds. The second chapter is devoted to the case study, namely, the stochastic analysis of a square cup deep drawing problem. Finally, the basic equations governing the mechanical modelling of the stamping process are summarized in the appendix.

  Abstract
  Stochastic Mechanics is a rapidly growing area of research, whose importance is being recognized not only in academic circles but also in industrial practice. This is no doubt [...]

  • 13
  •  read

• Stabilized finite element approximation of transient incompressible flows using orthogonal subscales

  R. Codina

  (2000). Research Report, Nº PI185

  
  

  Abstract

  In this paper we present a stabilized finite element formulation to solve the Oseen equations as a model problem involving both convection effects and the incompressibility restriction. The need for stabilization techniques to solve this problem arises because of the restriction in the possible choices for the velocity and pressure spaces dictated by the inf-sup condition, as well as the instabilities encountered when convection is dominant. Both can be overcome by resorting from the standard Galerkin method to a stabilized formulation. The one presented here is based on the subgrid scale concept, in which unresolvable scales of the continuous solution are approximately accounted for. In particular, the approach developed herein is based on the assumption that unresolved subscales are orthogonal to the finite element space. The motivation of the method is fully described. It is also shown that this formulation is table and optimally convergent for an adequate choice of the algorithmic on which the method depends.

  Abstract
  In this paper we present a stabilized finite element formulation to solve the Oseen equations as a model problem involving both convection effects and the incompressibility [...]

  • 16
  •  read

• Pressure stability in fractional step finite element methods for incompressible flows

  R. Codina

  (2000). Research Report, Nº PI186

  
  

  Abstract

  The objective of this paper is to analyze the pressure stability of fractional step finite element methods for incompressible flows. For the classical first order projection method, it is shown that there is a pressure control which depends on the time step size, and therefore there is a lower bound for this time step for stability reasons. The situation is much worse for a second order scheme in which part of the extremely weak. To overcome these shortcomings, a stabilized fractional step finite element method is also considered, and its stability is analyzed. Some simple numerical examples are presented to support the theoretical results.

  Abstract
  The objective of this paper is to analyze the pressure stability of fractional step finite element methods for incompressible flows. For the classical first order projection [...]

  • 10
  •  read

• Implementation of a stabilized finite element formulation for the incompressible Navier-Stokes equations based on a pressure gradient projection

  R. Codina, J. Blasco, G. Buscaglia, A. Huerta

  (2000). Research Report, Nº PI187

  
  

  Abstract

  We discuss in this paper some implementation aspects of a finite element formulation for the incompressible Navier-Stokes which allows the use of equal order velocity-pressure interpolations. The method consists in introducing the project of the pressure gradient and adding the difference between the pressure Laplacian and divergence of this new field to the incompressibility equations, both multiplied by suitable algorithmic parameters. The main purpose of this paper is to discuss how to deal with the new variable in the implementation of the algorithm. Obviously, it could be treated as one extra unknown, either explicitly or as a condensed variable. However, we take for granted that the only way or another. Here we discuss some iterative schemes to perform this uncoupling of the pressure gradient projection (PGP) from the calculation of the velocity and the pressure, both for stationary and the transient of the linearization loop and the iterative segregation of the PGP, whereas in the second the main dilemma concerns the explicit or implicit treatment of the PGP.  

  Abstract
  We discuss in this paper some implementation aspects of a finite element formulation for the incompressible Navier-Stokes which allows the use of equal order velocity-pressure [...]

  • 15
  •  read

• Prediction of damage and failure in civil engineering structures using a finite element model.

  E. Oñate, A. Hanganu, J. Canet

  (2000). Research Report, Nº PI188

  
  

  Abstract

  The paper describes a finite element damage model for non linear analysis of concrete or reinforced structures. It is show how can be effectively used to predict local and global damage up to structural failure. Examples of applications of the model to the analysis of different structures such as a nuclear containment shell, a housing building and the domes of St. Mark Basilica are presented.

  Abstract
  The paper describes a finite element damage model for non linear analysis of concrete or reinforced structures. It is show how can be effectively used to predict local and [...]

  • 24
  •  read

• Rigid-flexible body analysis without rotation parameters

  R. Taylor

  (2000). Research Report, Nº PI189

  
  

  Abstract

  This paper addresses the coupled flexible and rigid body response of solids undergoing large motions and deformations. The formulation is presented in a form which is free of rotation parameters and thus avoids the need of finding compatible integration formulae for translations and rotations. The motion of each rigid body is represented in terms of the nodal parameters of a simplex element subject to constraints which ensure rigid motions. Flexible structural members for rods and shells are then expressed in terms of displacement and relative displacement parameters leading to total system of equations involving only translation degrees of freedom. The motions are integrated using classical energy and momentum conserving schemes, thus leading to systems which are unconditionally stable for Hamiltonian (elastic-rigid) systems. The only requirement for absolute stability is that the non-linear algebraic equations representing the solution at each time step must converge.  The formulation is illustrated in two dimensions by representative numerical problems which involve both rigid and flexible parts or multi-rigid body situations. In all cases the conservation properties are observed.

  Abstract
  This paper addresses the coupled flexible and rigid body response of solids undergoing large motions and deformations. The formulation is presented in a form which is free [...]

  • 22
  •  read

• Locking in the incompressible limit for the element free galerkin method

  A. Huerta, S. Fernández-Méndez

  (2000). Research Report, Nº PI190

  
  

  Abstract

  Volumetric locking (locking in the incompressible limit) for linear elastic isotropic materials is studied in the context of the Element Free Galerkin method. The modal analysis developed here shows here shows that the number of non-physical locking modes is independent of the dilation parameter (support of the interpolation functions). Thus increasing the dilation parameter does not suppress locking. Nevertheless, an increase in the dilation parameter does reduce the energy associated to the non-physical locking behavior. Although more locking modes are present in the Element Free Galerkin method with quadratic consistency than with the standard biquadratic finite element method. Finally, numerical examples are shown.    

  Abstract
  Volumetric locking (locking in the incompressible limit) for linear elastic isotropic materials is studied in the context of the Element Free Galerkin method. The modal analysis [...]

  • 7
  •  read

• Límites de los métodos numéricos

  E. Oñate

  (2000). Research Report, Nº PI191

  
  

  Abstract

  Se plantean en el artículo los límites de los métodos numéricos para resolver cualquier problema que afecte a la vida del hombre. Se discute la capacidad de la razón para expresar todos los problemas del universo en forma matemática, y la posibilidad de encontrar su solución en forma numérica. Finalmente, se especula sobre el alcance y posibilidad de los métodos numéricos en el amplio campo de las ciencias sociales

  Abstract
  Se plantean en el artículo los límites de los métodos numéricos para resolver cualquier problema que afecte a la vida del hombre. Se discute la [...]

  • 85
  •  read

• El bucle de los números

  E. Oñate

  (2000). Research Report, Nº PI192

  
  

  Abstract

  El artículo presenta una sucinta panorámica sobre la evolución de los métodos numéricos dese la antigua Babilonia hasta nuestros días. Se destaca como la máxima Pitagórica de que “todo es número”, adquiere plena actualidad en nuestros días, en que, con la ayuda de los ordenadores, podemos dar respuestas numéricas a prácticamente cualquier problema que afecte a la vida del hombre. Recorriendo ese bucle de los números se relata brevemente en el artículo como la humanidad ha evolucionado paralelamente en su aspiración de cuantificar los fenómenos de la naturaleza, y como el paso de los pueblos ha ido de la mano de los avances en expresar numéricamente la solución de sus problemas más cotidianos.  

  Abstract
  El artículo presenta una sucinta panorámica sobre la evolución de los métodos numéricos dese la antigua Babilonia hasta nuestros días. [...]

  • 10
  •  read