Abstract

We present a Lagrangian formulation for finite element analysis of quasi-incompressible fluids that has excellent mass preservation features. The success of the formulation lays on a new residual-based stabilized expression of the mass balance equation obtained using the Finite Calculus (FIC) method. The governing equations are discretized with the FEM using simplicial elements with equal linear interpolation for the velocities and the pressure. The merits of the formulation in terms of reduced mass loss and overall accuracy are verified in the solution of 2D and 3D quasi-incompressible free-surface flow problems using the Particle Finite Element Method (PFEM, www.cimne.com/pfem). Examples include the sloshing of water in a tank, the collapse of one and two water columns in rectangular and prismatic tanks and the falling of a water sphere into a cylindrical tank containing water. Copyright c 0000 John Wiley & Sons, Ltd.

Abstract

In this chapter we present recent advances in the Discrete Element Method (DEM) and in the coupling of the DEM with the Finite Element Method (FEM) for solving a variety of problems in non linear solid mechanics involving damage, plasticity and multifracture situations.

Abstract

We present a generalized Lagrangian formulation for analysis of industrial forming processes involving thermally coupled interactions between deformable continua. The governing equations for the deformable bodies are written in a unified manner that holds both for fluids and solids. The success of the formulation lays on a residual-based expression of the mass conservation equation obtained using the Finite Calculus (FIC) method that provides the necessary stability for quasi/fully incompressible situations. The governing equations are discretized with the FEM via a mixed formulation using simplicial elements with equal linear interpolation for the velocities, the pressure and the temperature. The merits of the formulation are demonstrated in the solution of 2D and 3D thermally-coupled forming processes using the Particle Finite Element Method (PFEM).

Abstract

We present a Lagrangian monolithic strategy for solving fluid-structure interaction (FSI) problems. The formulation is called Unified because fluids and solids are solved using the same solution scheme and unknown variables. The method is based on a mixed velocity-pressure formulation. Each time step increment is solved via an iterative partitioned two-step procedure. The Particle Finite Element Method (PFEM) is used for solving the fluid parts of the domain, while for the solid ones the Finite Element Method (FEM) is employed. Both velocity and pressure fields are interpolated using linear shape functions. For quasiincompressible materials, the solution scheme is stabilized via the Finite Calculus (FIC) method. The stabilized elements for quasi-incompressible hypoelastic solids and Newtonian fluids are called VPS/S-element and VPS/F-element, respectively. Other two non-stabilized elements are derived for hypoelastic solids. One is based on a Velocity formulation (V-element) and the other on a mixed Velocity-Pressure scheme (VP-element). The algorithms for coupling the solid elements with the VPS/F fluid element are explained in detail. The Unified formulation is validated by solving benchmark FSI problems and by comparing the numerical solution to the ones published in the literature.

Abstract

The present work introduces a new application of the Particle Finite Element Method (PFEM) for the modeling of excavation problems. Lagrangian descriptions of motion of the continuum medium are the natural way of describing motion in solid mechanics. Particle finite element methods are based in these solid mechanics settings that can treat large material deformations and rapidly changing boundaries. These capabilities are very suitable for modeling fluid motions and moving free surfaces. That is the reason that the most of the research and applications of PFEM can be found in the context of computational fluid dynamics (CFD) instead of solid mechanics. The satisfying results in modeling fluids have been the motivation to use this method in dynamic solid mechanics. The simulation of an excavation process is a non-linear dynamic problem. It contains geometrical, material and contact non-linearities. Modeling the contact process can be classiffied as the main difficulty. The simulation face the problems of: the detection of a changing geometry, the detection of contact between several solid domains, the estimation of correct interacting forces, the computation of the wear related to this contact forces and the removing of the material that has been excavated from the model. The PFEM have its fundamentals in the classical non linear finite element analysis. The formulation provides a foundation in the updated lagrangian formulation for solids. The dynamic problem is integrated using an implicit scheme. Remeshing strategies are employed in order to identify the boundary surfaces. A remeshing of the domain is introduced for the detection of rapidly changing boundaries. The Delaunay tessellation and the alpha-shape concept together, are used as a methodology to de_ne the boundary from a cloud of points. At the same time an interface mesh is created for contact recognition. A particular constitutive contact law has been developed for capturing contact normal and frictional forces using this mesh. By means of an Archard-type law, the excavation and damage caused in the ground is quantified. The erosion and wear parameters of the grounds under study determine the excavability. The behavior of the geomaterials is described using a Damage Model which model fracture in the continuum. The update and storage of the variables in particles have especial schemes for its correct treatment. The method and the meshing process is adapted for the simulation of an excavation problem. PFEM is presented as a very suitable tool for the treatment of excavation problem. The method gives solution for the analysis of all processes that derive from it. The method has a high versatility and a reasonable computational cost. The obtained results are really promising.