Antonio Huerta

Viscous flow with large free surface motion

An arbitrary Lagrangian-Eulerian (ALE) Petrov-Galerkin finite element technique is developed to study nonlinear viscous fluids under large free surface wave motion. A review of the kinematics and field equations from an arbitrary reference is presented and since the major challenge of the ALE description lies in the mesh rezoning algorithm, various methods, including a new mixed formulation, are developed to update the mesh and map the moving domain in a more rational manner. Moreover, the streamline-upwind/Petrov-Galerkin formulation is implemented to accurately describe highly convective free surface flows. The effectiveness of the algorithm is demonstrated on a tsunami problem, the dam-break problem where the Reynolds number is taken as high as 3000, and a large-amplitude sloshing problem.

Enrichment and coupling of the finite element and meshless methods

A mixed hierarchical approximation based on finite elements and meshless methods is presented. Two cases are considered. The first one couples regions where finite elements or meshless methods are used to interpolate: continuity and consistency is preserved. The second one enriches a finite element mesh with particles. Thus, there is no need to remesh in adaptive refinement processes. In both cases the same formulation is used, convergence is studied and examples are shown. Copyright

Efficient unstructured quadrilateral mesh generation

This work is devoted to the description of an algorithm for automatic quadrilateral mesh generation. The technique is based on a recursive decomposition of the domain into quadrilateral elements. This automatically generates meshes composed entirely by quadrilaterals over complex geometries (there is no need for a previous step where triangles are generated). A background mesh with the desired element sizes allows to obtain the preferred sizes anywhere in the domain. The final mesh can be viewed as the optimal one given the objective function is defined. The recursive algorithm induces an efficient data structure which optimizes the computer cost. Several examples are presented to show the efficiency of this algorithm. Copyright

ADAPTIVE FINITE ELEMENT STRATEGIES BASED ON ERROR ASSESSMENT

Two main ingredients are needed for adaptive finite element computations. First, the error of a given solution must be assessed, by means of either error estimators or error indicators. After that, a new spatial discretization must be defined via h-, p- or r-adaptivity. In principle, any of the approaches for error assessment may be combined with any of the procedures for adapting the discretization. However, some combinations are clearly preferable. The advantages and limitations of the various alternatives are discussed. The most adequate strategies are illustrated by means of several applications in solid mechanics. Copyright

Arbitrary Lagrangian-Eulerian formulation for fluid-rigid body interaction

A new formulation for two-dimensional fluid–rigid body interaction problems is developed. In particular, vortex-induced oscillations of a rigid body in viscous incompressible flow are studied. The incompressible Navier–Stokes equations are used to describe the motion of the fluid, while it is assumed that the rigid body is mounted on a system consisting of a spring and a dashpot. An arbitrary Lagrangian–Eulerian formulation (ALE) is used in order to account for large boundary motion. A general formulation for the coupled problem is obtained by uncoupling the translation motion of the body from its rotational motion and developing a specific algorithm to efficiently handle the nonlinear dependence of the rotations. This general formulation can be easily applied to multi-body problems. Two numerical examples involving either translations and rotations are presented as an illustration of the proposed methodologies for fluid–rigid body interaction.

Numerical differentiation for local and global tangent operators in computational plasticity

In this paper, numerical differentiation is applied to integrate plastic constitutive laws and to compute the corresponding consistent tangent operators. The derivatives of the constitutive equations are approximated by means of difference schemes. These derivatives are needed to achieve quadratic convergence in the integration at Gauss-point level and in the solution of the boundary value problem. Numerical differentiation is shown to be a simple, robust and competitive alternative to analytical derivatives. Quadratic convergence is maintained, provided that adequate schemes and stepsizes are chosen. This point is illustrated by means of some numerical examples.

Consistent tangent matrices for substepping schemes

A very simple and general expression of the consistent tangent matrix for substepping time-integration schemes is presented. If needed, the derivatives required for the computation of the consistent tangent moduli can be obtained via numerical differentiation. These two strategies (substepping and numerical differentiation) lead to quadratic convergence in complex nonlinear inelasticity problems.

ALE STRESS UPDATE FOR TRANSIENT AND QUASISTATIC PROCESSES

A key issue in Arbitrary Lagrangian–Eulerian (ALE) non-linear solid mechanics is the correct treatment of the convection terms in the constitutive equation. These convection terms, which reflect the relative motion between the finite element mesh and the material, are found for both transient and quasistatic ALE analyses. It is shown in this paper that the same explicit algorithms can be employed to handle the convection terms of the constitutive equation for both types of analyses. The most attractive consequence of this fact is that a quasistatic simulation can be upgraded from Updated Lagrangian (UL) to ALE without significant extra computational cost. These ideas are illustrated by means of two numerical examples.

Error estimation and adaptivity for nonlocal damage models

Nonlocal damage models are typically used to model failure of quasi-brittle materials. Due to brittleness, the choice of a particular model or set of parameters can have a crucial influence on the structural response. To assess this influence, it is essential to keep finite element discretization errors under control. If not, the effect of these errors on the result of a computation could be erroneously interpreted from a constitutive viewpoint. To ensure the quality of the FE solution, an adaptive strategy based on error estimation is proposed here. It is based on the combination of a residual-type error estimator and quadrilateral h-remeshing. Another important consequence of brittleness is that it leads to structural responses of the snap-through or snap-back type. This requires the use of arc-length control, with a definition of the arc parameter that accounts for the localized nature of quasi-brittle failure. All these aspects are discussed for two particular nonlocal damage models (Mazars and modified von Mises) and for two tests: the Brazilian tensile splitting test and the single-edge notched beam test. For the latter test, the capability of the Mazars model to capture the curved crack pattern observed in experiments – a topic of debate in the literature – is confirmed.

A unified approach to remeshing strategies for finite element h-adaptivity

In h-adaptivity, a remeshing strategy is needed to compute the distribution of required element size using the estimated error distribution. Several authors have introduced different remeshing strategies yielding very different results. In this work these methods are included in a unified framework, emphasizing the role of the underlying hypotheses. Moreover, an objective tool to evaluate the accuracy of the resulting finite element solution is presented. Thus, a new remeshing strategy is introduced to optimize the accuracy of the adapted solutions. The different remeshing strategies are compared with well-known numerical examples.