JOURNALS

Documents published in Scipedia

• Comparative study on sheet metal forming processes by numerical modelling and experiment

  W. Sosnowski, E. Oñate, C. Agelet de Saracibar

  Journal of Materials Processing Technology (1992). Vol. 34 (1-4), pp. 109-116

  
  

  Abstract

  In this paper some results of a wide experimental program are presented and compared with some finite element solution of sheet metal forming problems using a viscous shell formulation.

  Abstract
  In this paper some results of a wide experimental program are presented and compared with some finite element solution of sheet metal forming problems using a viscous shell [...]

  • 20
  •  read

• High‐density sphere packing for discrete element method simulations

  C. Labra, E. Oñate

  Commun. Numer. Meth. Engng (2009). Vol. 25 (7), pp. 837-849

  
  

  Abstract

  The first step in a discrete element simulation is the discretization of the domain into a set of particles. The cost of generating a good cylindrical or spherical packing has resulted in a great number of approaches during the last years. A new algorithm is proposed for high‐density packing using a scheme that minimizes the distance between each particle. Using the support of a finite element mesh, less time is needed in order to achieve a low porosity configuration. In addition, a boundary constraint is introduced. The application of the same optimization scheme is used as a condition to force a good surface definition. The results obtained show a high efficiency for the generation of low porosity packing, achieving values smaller than 10% in 2D cases and 30% in 3D cases. 

  Abstract
  The first step in a discrete element simulation is the discretization of the domain into a set of particles. The cost of generating a good cylindrical or spherical packing [...]

  • 34
  •  read

• Parallel Delaunay triangulation for particle finite element methods

  Y. Fragakis, E. Oñate

  Commun. Numer. Meth. Engng (2008). Vol. 24 (11), pp. 1009-1017

  
  

  Abstract

  Delaunay triangulation is a geometric problem that is relatively difficult to parallelize. Parallel algorithms are usually characterized by considerable interprocessor communication or important serialized parts. In this paper, we propose a method that achieves high speed‐ups, but needs information regarding locally maximum element circumspheres prior to the beginning of the algorithm. Such information is directly available in iterative methods, like the particle finite element methods. The developed parallel Delaunay triangulation method, has minimum communication requirements, is quite simple and achieves high parallel efficiency.

  Abstract
  Delaunay triangulation is a geometric problem that is relatively difficult to parallelize. Parallel algorithms are usually characterized by considerable interprocessor communication [...]

  • 14
  •  read

• A simple method for automatic update of finite element meshes

  G. Chiandussi, G. Bugeda, E. Oñate

  Commun. Numer. Meth. Engng (2000). Vol. 16 (1), pp. 1-19

  
  

  Abstract

  A simple method to automatically update the finite element mesh of the analysis domain is proposed. The method considers the mesh as a fictitious elastic body subjected to prescribed displacements at selected boundary points. The mechanical properties of each mesh element are appropriately selected in order to minimize the deformation and the distortion of the mesh elements. Different selection strategies have been used and compared in their application to simple examples. The method avoids the use of remeshing in the solution of shape optimization problems and reduces the number of remeshing steps in the solution of coupled fluid–structure interaction problems

  Abstract
  A simple method to automatically update the finite element mesh of the analysis domain is proposed. The method considers the mesh as a fictitious elastic body subjected to [...]

  • 32
  •  read

• An advancing front point generation technique

  R. Löhner, E. Oñate

  Commun. Numer. Meth. Engng (1998). Vol. 14 (12), pp. 1097-1108

  
  

  Abstract

  An algorithm to construct boundary‐conforming, isotropic clouds of points with variable density in space is described. The input required consists of a specified mean point distance and an initial triangulation of the surface. Borrowing a key concept from advancing front grid generators, one point at a time is removed and, if possible, surrounded by admissible new points. This operation is repeated until no active points are left. Timings show that the scheme is about an order of magnitude faster than volume grid generators based on the advancing front technique, making it possible to generate large (>106) yet optimal clouds of points in a matter of minutes on a workstation. Several examples are included that demonstrate the capabilities of the technique. 

  Abstract
  An algorithm to construct boundary‐conforming, isotropic clouds of points with variable density in space is described. The input required consists of a specified mean point [...]

  • 43
  •  read

• Mixing anisotropic formulation for analysis of composites

  S. Oller, E. Oñate, J. Canet

  Commun. Numer. Meth. Engng (1996). Vol. 12 (8), pp. 471-482

  
  

  Abstract

  A general constitutive model adequate for analysis of the thermomechanical response of composite materials is presented. The model is based on the mixture of the basic substances of the composite and allows the evaluation of the interdependence between the constitutive behaviour of different compounding materials. The behaviour of the each compound is modelled by a general anisotropic thermo‐elasto‐plastic model, termed the ‘base model’. The different base models for each compound are combined using mixing theory to simulate the behaviour of the multiphase material.

  Abstract
  A general constitutive model adequate for analysis of the thermomechanical response of composite materials is presented. The model is based on the mixture of the basic substances [...]

  • 15
  •  read

• On the pseudomaterial approach for the analysis of transient forming processes

  M. Cruchaga, E. Oñate, S. Idelsohn

  Commun. Numer. Meth. Engng (1995). Vol. 11 (2), pp. 137-148

  
  

  Abstract

  A fixed‐mesh method for the analysis of transient forming processes is presented. The mesh covers material regions and zones through which the material may flow. These last zones are identified by a pseudomaterial with relatively small physical parameters. During time processing, the interface between both materials is followed by an arbitrary Lagrangian mesh. This technique appears to be suitable for the treatment of moving surfaces with sharp corners. A particular boundary condition for the Navier‐Stokes equations is also introduced in order to model a porous wall.

  Abstract
  A fixed‐mesh method for the analysis of transient forming processes is presented. The mesh covers material regions and zones through which the material may flow. These last [...]

  • 23
  •  read

• Comparative study of different discrete element models and evaluation of equivalent micromechanical parameters

  J. Rojek, C. Labra, O. Su, E. Oñate

  Int. J. Solids and Structures (2012). Vol 49 (13), pp. 1497-1517

  
  

  Abstract

  Comparative studies of different discrete element models of a rock-type material are presented. The discrete element formulation employs spherical particles with the cohesive interaction model combining linear elastic behaviour with brittle failure. Numerical studies consisted in simulation of the uniaxial compression test. Two cylindrical specimens with particle size distributions yielding different degree of heterogeneity have been used. Macroscopic response produced by different discrete element models has been compared. The main difference between the compared models consists in the evaluation of micromechanical constitutive parameters. Two approaches are compared. In the first approach, the contact stiffness and strength parameters depend on the local particle size, while in the second approach, global uniform contact parameters are assumed for all the contacting pairs in function of average geometric measures characterizing the particle assembly. The size dependent contact parameters are calculated as functions of geometric parameters characterizing each contacting particle pair. As geometric scaling parameters, the arithmetic and harmonic means, as well as the minimum of the radii of two contacting particles are considered. Two different models with size dependent contact parameters are formulated. The performance of these models is compared with that of the discrete element model with global uniform contact parameters. Equivalence between the models with size dependent and uniform contact parameters has been checked. In search of this equivalence, different methods of evaluation of global uniform parameters have been studied. The contact stiffness has been evaluated in terms of the average radius of the particle assembly or in terms of the averages of the arithmetic and harmonic means of the contact pair radii, the geometric parameters used in the evaluation of the contact stiffness in the size-dependent models. The uniform contact strengths have been determined as functions of the averages of radii squares, squares of arithmetic radii means or squares of minimum radii of the contacting pairs. For the more homogenous specimen, the models with local size dependent parameters and models with global uniform parameters give similar response. The models with uniform parameters evaluated according to the averages of the geometric parameters used in the evaluation of local parameters ensure better agreement with the respective models with size-dependent parameters than the models with uniform parameters evaluated according to the particle radii. Simulations using the more heterogenous specimen reveal differences between the considered models. There are significant differences in stress–strain curves as well as in the failure pattern. The models with local size-dependent parameters are more sensitive to the change of heterogeneity than the model with global uniform parameters.

  Abstract
  Comparative studies of different discrete element models of a rock-type material are presented. The discrete element formulation employs spherical particles with the cohesive [...]

  • 30
  •  read

• Numerical simulation of fiber reinforced composite materials––two procedures

  E. Car, F. Zalamea, S. Oller, J. Canet, E. Oñate

  Int. J. Solids and Structures (2002). Vol. 39 (7), pp. 1967-1986

  
  

  Abstract

  In this work, two methodologies for the analysis of unidirectional fiber reinforced composite materials are presented. The first methodology used is a generalized anisotropic large strains elasto-plastic constitutive model for the analysis of multiphase materials. It is based on the mixing theory of basic substance. It is the manager of the several constitutive laws of the different compounds and it allows to consider the interaction between the compounds of the composite materials. In fiber reinforced composite materials, the constitutive behavior of the matrix is isotropic, whereas the fiber is considered orthotropic. So, one of the constitutive model used in the mixing theory needs to consider this characteristic. The non-linear anisotropic theory showed in this work is a generalization of the classic isotropic plasticity theory (A Continuum Constitutive Model to Simulate the Mechanical Behavior of Composite Materials, PhD Thesis, Universidad Politécnica de Cataluña, 2000). It is based in a one-to-one transformation of the stress and strain spaces by means of a four rank tensor. The second methodology used is based on the homogenization theory. This theory divided the composite material problem into two scales: macroscopic and microscopic scale. In macroscopic level the composite material is assuming as a homogeneous material, whereas in microscopic level a unit volume called cell represents the composite (Tratamiento Numérico de Materiales Compuestos Mediante la teorı́a de Homogeneización, PhD Thesis, Universidad Politécnica, de Cataluña 2001). This formulation presents a new viewpoint of the homogenization theory in which can be found the equations that relate both scales. The solution is obtained using a coupled parallel code based on the finite elements method for each scale problem.

  Abstract
  In this work, two methodologies for the analysis of unidirectional fiber reinforced composite materials are presented. The first methodology used is a generalized anisotropic [...]

  • 37
  •  read

• A coupled thermomechanical model for the solidification of cast metals

  D. Celentano, S. Oller, E. Oñate

  Int. J. Solids and Structures (1996). Vol. 33 (5), pp. 647-673

  
  

  Abstract

  A coupled thermomechanical model to simulate solidification problems in casting is presented. The model is formulated from a phenomenological point of view using a general isotropic thermoelasto-plasticity theory. Generalized phase-change effects accounting for the different thermomechanical behaviour of the solidifying material during its evolution from liquid to solid have been considered. For this purpose, a phase-change function, plastic evolution equations and a temperature-dependent constitutive law have been defined. Full thermomechanical effects as well as variable thermal and mechanical boundary conditions are also taken into account. Particular details concerning the numerical implementation of the model are given, where special emphasis is devoted to the resulting highly non-linear fully coupled finite element equations. The behaviour of this formulation is studied first in a simple quenching problem. Finally, a cylindrical casting test problem including phase-change phenomena, temperature-dependent constitutive properties and contact effects is analysed. Numerical results are compared with laboratory measurements.

  Abstract
  A coupled thermomechanical model to simulate solidification problems in casting is presented. The model is formulated from a phenomenological point of view using a general [...]

  • 14
  •  read