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 mixed velocity-pressure finite element formulation for solving the updated Lagrangian equations for quasi and fully incompressible fluids. Details of the governing equations for the conservation of momentum and mass are given in both differential and variational form. The finite element interpolation uses an equal order approximation for the velocity and pressure unknowns. The procedure for obtaining stabilized FEM solutions is outlined. The solution in time of the discretized governing conservation equations using an incremental iterative segregated scheme is described. The linearization of these equations and the derivation of the corresponding tangent stiffness matrices is detailed. Other iterative schemes for the direct computation of the nodal velocities and pressures at the updated configuration are presented. The advantages and disadvantages of choosing the current or the updated configuration as the reference configuration in the Lagrangian formulation are discussed.

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 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.

Abstract

This paper presents a computational framework for the numerical analysis of fluid-saturated porous media at large strains. The proposal relies, on one hand, on the particle finite element method (PFEM), known for its capability to tackle large deformations and rapid changing boundaries, and, on the other hand, on constitutive descriptions well established in current geotechnical analyses (Darcy’s law; Modified Cam Clay; Houlsby hyperelasticity). An important feature of this kind of problem is that incompressibility may arise either from undrained conditions or as a consequence of material behaviour; incompressibility may lead to volumetric locking of the low-order elements that are typically used in PFEM. In this work, two different three-field mixed formulations for the coupled hydromechanical problem are presented, in which either the effective pressure or the Jacobian are considered as nodal variables, in addition to the solid skeleton displacement and water pressure. Additionally, several mixed formulations are described for the simplified single-phase problem due to its formal similitude to the poromechanical case and its relevance in geotechnics, since it may approximate the saturated soil behaviour under undrained conditions. In order to use equal-order interpolants in displacements and scalar fields, stabilization techniques are used in the mass conservation equation of the biphasic medium and in the rest of scalar equations. Finally, all mixed formulations are assessed in some benchmark problems and their performances are compared. It is found that mixed formulations that have the Jacobian as a nodal variable perform better.

Abstract

The treatment of (near-)incompressibility is a major concern for the simulation of rubber-like parts, or forming processes. The use of mixed finite element methods is known to prevent the locking of the F.E. approximation in the incompressible limit. However, the stability of these formulations is conditionned by the fullfilment of the inf-sup condition. Recently, finite elements method has evolved with the introduction of the partition of unity. The X-FEM uses it to remove the need to mesh (and remesh) physical surfaces. In this paper, a strategy is proposed for the treatment of holes within X-FEM in the incompressible setting. Numerical examples show that F.E. convergence rate is preserved and that the inf-sup condition is passed.

Abstract

Thermoplastic composites are widely considered in structural parts. In this paper attention is paid to squeeze flow of continuous fiber laminates. In the case of unidirectional prepregs, the ply constitutive equation is modeled as a transversally isotropic fluid, that must satisfy both the fiber inextensibility as well as the fluid incompressibility. When laminate is squeezed the flow kinematics exhibits a complex dependency along the laminate thickness requiring a detailed velocity description through the thickness. In a former work the solution making use of an in-plane-out-of-plane separated representation within the PGD – Poper Generalized Decomposition – framework was successfully accomplished when both kinematic constraints (inextensibility and incompressibility) were introduced using a penalty formulation for circumventing the LBB constraints. However, such a formulation makes difficult the calculation on fiber tractions and compression forces, the last required in rheological characterizations. In this paper the former penalty formulation is substituted by a mixed formulation that makes use of two Lagrange multipliers, while addressing the LBB stability conditions within the separated representation framework, questions never until now addressed.

Abstract

The first superconvergent hybridisable discontinuous Galerkin method for linear elastic problems capable of using the same degree of approximation for both the primal and mixed variables is presented. The key feature of the method is the strong imposition of the symmetry of the stress tensor by means of the well known and extensively used Voigt notation, circumventing the use of complex mathematical concepts to enforce the symmetry of the stress tensor either weakly or strongly. A novel procedure to construct element by element a superconvergent postprocessed displacement is proposed. Contrary to other hybridisable discontinuous Galerkin formulations, the methodology proposed here is able to produce a superconvergent displacement field for low-order approximations. The resulting method is robust and locking-free in the nearly incompressible limit. An extensive set of numerical examples is utilised to provide evidence of the optimality of the method and its superconvergent properties in two and three dimensions and for different element types

Abstract

In this work a stabilized mixed formulation for the solution of non-linear solid mechanics problems in nearly-incompressible conditions is presented. In order to deal with high material deformation, an implicit Material Point Method is chosen. Such choice allows avoiding the classical limitations of the Finite Element Method, e.g., element tangling and extreme mesh distortion. The proposed mixed formulation, with displacement and pressure as primary variables, is tested through classical benchmarks in solid and geo-mechanics where a Neo-Hookean, a J2 and a Mohr-Coulomb plastic law are employed. Further, the stabilized mixed formulation is compared with a displacement-based formulation to demonstrate how the proposed approach gets better results in terms of accuracy, not only when incompressible materials are simulated, but also in the case of compressible ones.

Abstract

Low-order finite elements face inherent limitations related to their poor convergence properties. Such difficulties typically manifest as mesh-dependent or excessively stiff behaviour when dealing with complex problems. A recent proposal to address such limitations is the adoption of mixed displacement-strain technologies which were shown to satisfactorily address both problems. Unfortunately, although appealing, the use of such element technology puts a large burden on the linear algebra, as the solution of larger linear systems is needed. In this paper, the use of an explicit time integration scheme for the solution of the mixed strain-displacement problem is explored as an alternative. An algorithm is devised to allow the effective time integration of the mixed problem. The developed method retains second order accuracy in time and is competitive in terms of computational cost with the standard irreducible formulation.