0
Research Papers

A Total Lagrangian ANCF Liquid Sloshing Approach for Multibody System Applications

[+] Author and Article Information
Cheng Wei

Department of Aerospace Engineering,
Harbin Institute of Technology,
92 W. Dazhi Street,
359 Mail Box, Heilongjiang,
Harbin 150001, China
e-mail: weicheng@hit.edu.cn

Liang Wang

Department of Mechanical and
Industrial Engineering,
University of Illinois at Chicago,
842 W. Taylor Street,
Chicago, IL 60607

Ahmed A. Shabana

Department of Mechanical and
Industrial Engineering,
University of Illinois at Chicago,
842 W. Taylor Street,
Chicago, IL 60607
e-mail: shabana@uic.edu

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received July 19, 2014; final manuscript received September 27, 2014; published online April 16, 2015. Assoc. Editor: Ahmet S. Yigit.

J. Comput. Nonlinear Dynam 10(5), 051014 (Sep 01, 2015) (10 pages) Paper No: CND-14-1171; doi: 10.1115/1.4028720 History: Received July 19, 2014; Revised September 27, 2014; Online April 16, 2015

The objective of this investigation is to develop a total Lagrangian nonincremental liquid sloshing solution procedure based on the finite element (FE) absolute nodal coordinate formulation (ANCF). The proposed liquid sloshing modeling approach can be used to avoid the difficulties of integrating most of fluid dynamics formulations, which are based on the Eulerian approach, with multibody system (MBS) dynamics formulations, which are based on a total Lagrangian approach. The proposed total Lagrangian FE fluid dynamics formulation, which can be systematically integrated with computational MBS algorithms, differs significantly from the conventional FE or finite volume methods which are based on an Eulerian representation that employs the velocity field of a fixed control volume in the region of interest. The ANCF fluid equations are expressed in terms of displacement and gradient coordinates of material points, allowing for straightforward implementation of kinematic constraint equations and for the systematic modeling of the interaction of the fluid with the external environment or with rigid and flexible bodies. The fluid incompressibility conditions and surface traction forces are considered and derived directly from the Navier–Stokes equations. Two ANCF brick elements, one is obtained using an incomplete polynomial representation and the other is obtained from a B-spline volume representation, are used. The new approach ensures the continuity of the displacement gradients at the nodal points and allows for imposing higher degree of continuity across the element interface by applying algebraic constraint equations that can be used to eliminate dependent variables and reduce the model dimensionality. Regardless of the magnitude of the fluid displacement, the fluid has a constant mass matrix, leading to zero Coriolis and centrifugal forces. The analysis presented in this paper demonstrates the feasibility of developing an efficient nonincremental total Lagrangian approach for modeling sloshing problems in MBS system applications in which the bodies can experience large displacements including finite rotations. Several examples are presented in order to shed light on the potential of using the ANCF liquid sloshing formulation developed in this study.

FIGURES IN THIS ARTICLE
<>
Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Fig. 1

The 8-node brick fluid element

Grahic Jump Location
Fig. 2

Continuity at the element interface

Grahic Jump Location
Fig. 3

C1 Interface discontinuity in two-IPAE mesh

Grahic Jump Location
Fig. 4

C1 Continuity after applying continuity constraints

Grahic Jump Location
Fig. 6

Boundary conditions

Grahic Jump Location
Fig. 7

Fluid element in a container

Grahic Jump Location
Fig. 8

Fluid/ground surface interaction using one element (a) IPAE example and (b) BSAE example

Grahic Jump Location
Fig. 9

Fluid/ground surface interaction using an eight-IPAE mesh: (a) Without continuity constraints and (b) with continuity constraints

Grahic Jump Location
Fig. 10

The Z component of rx of node number 10 in eight-IPAE mesh (—▪— continuity, Δ without continuity)

Grahic Jump Location
Fig. 11

Effect of surface tension using eight-IPAE mesh

Grahic Jump Location
Fig. 12

Sloshing problem solution using one element [23]

Grahic Jump Location
Fig. 13

Sloshing problem using the FFR formulation

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In