skip to main content


Title: Performance Evaluation of ANCF Tetrahedral Elements in the Analysis of Liquid Sloshing
Abstract The performance of the absolute nodal coordinate formulation (ANCF) tetrahedral element in the analysis of liquid sloshing is evaluated in this paper using a total Lagrangian nonincremental solution procedure. In this verification study, the results obtained using the ANCF tetrahedral element are compared with the results of the ANCF solid element which has been previously subjected to numerical verification and experimental validation. The tetrahedral-element model, which allows for arbitrarily large displacements including rotations, can be systematically integrated with computational multibody system (MBS) algorithms that allow for developing complex sloshing/vehicle models. The new fluid formulation allows for systematically increasing the degree of continuity in order to obtain higher degree of smoothness at the element interface, eliminate dependent variables, and reduce the model dimensionality. The effect of the fluid/container interaction is examined using a penalty contact approach. Simple benchmark problems and complex railroad vehicle sloshing scenarios are used to examine the performance of the ANCF tetrahedral element in solving liquid sloshing problems. The simulation results show that, unlike the ANCF solid element, the ANCF tetrahedral element model exhibits nonsmoothness of the free surface. This difference is attributed to the gradient discontinuity at the tetrahedral-element interface, use of different meshing rules for the solid- and tetrahedral-elements, and the interaction between elements. It is shown that applying curvature-continuity conditions leads, in general, to higher degree of smoothness. Nonetheless, a higher degree of continuity does not improve the solution accuracy when using the ANCF tetrahedral elements.  more » « less
Award ID(s):
1632302
NSF-PAR ID:
10290993
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Verification, Validation and Uncertainty Quantification
Volume:
5
Issue:
3
ISSN:
2377-2158
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The continuity of the position-vector gradients at the nodal points of a finite element mesh does not always ensure the continuity of the gradients at the element interfaces. Discontinuity of the gradients at the interface not only adversely affects the quality of the simulation results, but can also lead to computer models that do not properly represent realistic physical system behaviors, particularly in the case of soft and fluid material applications. In this study, the absolute nodal coordinate formulation (ANCF) finite elements are used to define general curvature-continuity conditions that allow for eliminating or minimizing the discontinuity of the position gradients at the element interface. For the ANCF solid element, with four-node surfaces, it is shownthat continuity of the gradients tangent to an arbitrary point on a surface is ensured as the result of the continuity of the gradients at the nodal points. The general ANCF continuity conditions are applicable to both reference-configuration straight and curved geometries. These conditions are formulated without the need for using the computer-aided-design knot vector and knot multiplicity, which do not account properly for the concept of system degrees of freedom. The ANCF curvature-continuity conditions are written in terms of constant geometric coefficients obtained using the matrix of position-vector gradients that defines the reference-configuration geometry. The formulation of these conditions is demonstrated using the ANCF fully parameterized three-dimensional solid and tetrahedral elements, which employ a complete set of position gradients as nodal coordinates. Numerical results are presented in order to examine the effect of applying the curvature-continuity conditions on achieving a higher degree of smoothness at the element interfaces in the case of soft and fluid materials. 
    more » « less
  2. .The objective of this investigation is to develop a new total Lagrangian continuum-based liquid sloshing model that can be systematically integrated with multibody system (MBS) algorithms in order to allow for studying complex motion scenarios. The new approach allows for accurately capturing the effect of the sloshing forces during curve negotiation, rapid lane change, and accelerating and braking scenarios. In these motion scenarios, the liquid experiences large displacements and significant changes in shape that can be captured effectively using the finite element (FE) absolute nodal coordinate formulation (ANCF). ANCF elements are used in this investigation to describe complex mesh geometries, to capture the change in inertia due to the change in the fluid shape, and to accurately calculate the centrifugal forces, which for flexible bodies do not take the simple form used in rigid body dynamics. A penalty formulation is used to define the contact between the rigid tank walls and the fluid. A fully nonlinear MBS truck model that includes a suspension system and Pacejka’s brush tire model is developed. Specified motion trajectories are used to examine the vehicle dynamics in three different scenarios – deceleration during straight-line motion, rapid lane change, and curve negotiation. It is demonstrated that the liquid sloshing changes the contact forces between the tires and the ground – increasing the forces on certain wheels and decreasing the forces on other wheels. In cases of extreme sloshing, this dynamic behavior can negatively impact the vehicle stability by increasing the possibility of wheel lift and vehicle rollover. 
    more » « less
  3. A high-order in space spectral-element methodology for the solution of a strongly coupled fluid-structure interaction (FSI) problem is developed. A methodology is based on a partitioned solution of incompressible fluid equations on body-fitted grids, and nonlinearly-elastic solid deformation equations coupled via a fixed-point iteration approach with Aitken relaxation. A comprehensive verification strategy of the developed methodology is presented, including h-, p-and temporal refinement studies. An expected order of convergence is demonstrated first separately for the corresponding fluid and solid solvers, followed by a self-convergence study on a coupled FSI problem (self-convergence refers to a convergence to a reference solution obtained with the same solver at higher resolution). To this end, a new three-dimensional fluid-structure interaction benchmark is proposed for a verification of the FSI codes, which consists of a fluid flow in a channel with one rigid and one flexible wall. It is shown that, due to a consistent problem formulation, including initial and boundary conditions, a high-order spatial convergence on a fully coupled FSI problem can be demonstrated. Finally, a developed framework is applied successfully to a Direct Numerical Simulation of a turbulent flow in a channel interacting with a compliant wall, where the fluid-structure interface is fully resolved. 
    more » « less
  4. null (Ed.)
    Abstract The convergence characteristics of three geometrically accurate spatial finite elements (FEs) are examined in this study using an eigenvalue analysis. The spatial beam, plate, and solid elements considered in this investigation are suited for both structural and multibody system (MBS) applications. These spatial elements are based on geometry derived from the kinematic description of the absolute nodal coordinate formulation (ANCF). In order to allow for an accurate reference-configuration geometry description, the element shape functions are formulated using constant geometry coefficients defined using the position-vector gradients in the reference configuration. The change in the position-vector gradients is used to define a velocity transformation matrix that leads to constant element inertia and stiffness matrices in the case of infinitesimal rotations. In contrast to conventional structural finite elements, the elements considered in this study can be used to describe the initial geometry with the same degree of accuracy as B-spline and nonuniform rational B-spline (NURBS) representations, widely used in the computer-aided design (CAD). An eigenvalue analysis is performed to evaluate the element convergence characteristics in the case of different geometries, including straight, tapered, and curved configurations. The frequencies obtained are compared with those obtained using a commercial FE software and analytical solutions. The stiffness matrix is obtained using both the general continuum mechanics (GCM) approach and the newly proposed strain split method (SSM) in order to investigate its effectiveness as a locking alleviation technique. 
    more » « less
  5. Origami — the ancient art of paper folding — has been widely adopted as a design and fabrication framework for many engineering applications, including multi-functional structures, deployable spacecraft, and architected materials. These applications typically involve complex and dynamic deformations in the origami facets, necessitating high-fidelity models to better simulate folding-induced mechanics and dynamics. This paper presents the formulation and validation of such a new model based on the Absolute Nodal Coordinate Formulation (ANCF), which exploits the tessellated nature of origami and describes it as an assembly of flexible panels rotating around springy creases. To estimate the crease folding, we mathematically formulate a “torsional spring connector” in the framework of ANCF and apply it to the crease nodes, where the facets meshed by ANCF plate elements are interconnected. We simulate the dynamic folding of a Miura-ori unit cell and compare the results with commercial finite element software (ABAQUS) to validate the modeling accuracy. The ANCF model can converge using significantly fewer elements than ABAQUS without sacrificing accuracy. Therefore, this high-fidelity model can help deepen our knowledge of folding-induced mechanics and dynamics, broadening the appeals of origami in science and engineering. 
    more » « less