skip to main content

Title: Nonlinear finiteelementanalysisofliquidsloshingincomplex vehiclemotionscenarios
.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 more » 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. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of sound and vibration
Page Range or eLocation-ID:
208 - 233
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 formore »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.« less
  2. continuum-based approach for simultaneously controlling the motion and shape of soft robots and materials (SRM) is proposed. This approach allows for systematically computing the actuation forces for arbitrary desired SRM motion and geometry. In order to control both motion and shape the position and position gradients of the absolute nodal coordinate formulation (ANCF) are used to formulate rheonomic specified trajectory and shape constraint equations, used in an inverse dynamics procedure to define the actuation control forces. Unlike control of rigid-body systems which requires a number of independent actuation forces equal to the number of the joint coordinates, the SRM motion/shape control leads to generalized control forces which need to be interpreted differently in order to properly define the actuation forces. While the definition of these motion/shape control forces is demonstrated using air pressure actuation commonly used in the SRM control, the proposed procedure can be applied to other SRM actuation types. The approaches for determining the actuation pressure in the two cases of space-dependent and constant pressures are outlined. Effect of the change in the surface geometry on the actuation pressure is accounted for using Nanson’s formula. The obtained numerical results demonstrate that the motion and shape can be simultaneouslymore »controlled using the new actuation force definitions.« less
  3. The rising global trend to reduce dependence on fossil fuels has provided significant motivation toward the development of alternative energy conversion methods and new technologies to improve their efficiency. Recently, oscillating energy harvesters have shown promise as highly efficient and scalable turbines, which can be implemented in areas where traditional energy extraction and conversion are either unfeasible or cost prohibitive. Although such devices are quickly gaining popularity, there remain a number of hurdles in the understanding of their underlying fluid dynamics phenomena. The ability to achieve high efficiency power output from oscillating airfoil energy harvesters requires exploitation of the complexities of the event of dynamic stall. During dynamic stall, the oncoming flow separates at the leading edge of the airfoil to form leading ledge vortex (LEV) structures. While it is well known that LEVs play a significant role in aerodynamic force generation in unsteady animal flight (e.g. insects and birds), there is still a need to further understand their spatiotemporal evolution in order to design more effective energy harvesting enhancement mechanisms. In this work, we conduct extensive experimental investigations to shed-light on the flow physics of a heaving and pitching airfoil energy harvester operating at reduced frequencies of k =more »fc=U1 = 0.06-0.18, pitching amplitude of 0 = 75 and heaving amplitude of h0 = 0:6c. The experimental work involves the use of two-component particle image velocimetry (PIV) measurements conducted in a wind tunnel facility at Oregon State University. Velocity fields obtained from the PIV measurements are analyzed qualitatively and quantitatively to provide a description of the dynamics of LEVs and other flow structures that may be present during dynamic stall. Due to the difficulties of accurately measuring aerodynamic forces in highly unsteady flows in wind tunnels, a reduced-order model based on the vortex-impulse theory is proposed for estimating the aerodynamic loadings and power output using flow field data. The reduced-order model is shown to be dominated by two terms that have a clear physical interpretation: (i) the time rate of change of the impulse of vortical structures and (ii) the Kutta-Joukowski force which indirectly represents the history effect of vortex shedding in the far wake. Furthermore, the effects of a bio-inspired flow control mechanism based on deforming airfoil surfaces on the flow dynamics and energy harvesting performance are investigated. The results show that the aerodynamic loadings, and hence power output, are highly dependent on the formation, growth rate, trajectory and detachment of the LEV. It is shown that the energy harvesting efficiency increases with increasing reduced frequency, peaking at 25% when k = 0.14, agreeing very well with published numerical results. At this optimal reduced frequency, the time scales of the LEV evolution and airfoil kinematics are matched, resulting in highly correlated aerodynamic load generation and airfoil motion. When operating at k > 0:14, it is shown that the aerodynamic moment and airfoil pitching motion become negatively correlated and as a result, the energy harvesting performance is deteriorated. Furthermore, by using a deforming airfoil surface at the leading and trailing edges, the peak energy harvesting efficiency is shown to increase by approximately 17% and 25% relative to the rigid airfoil, respectively. The performance enhancement is associated with enhanced aerodynamic forces for both the deforming leading and trailing edges. In addition, The deforming trailing edge airfoil is shown to enhance the correlation between the aerodynamic moment and pitching motion at higher reduced frequencies, resulting in a peak efficiency at k = 0:18 as opposed to k = 0:14 for the rigid airfoil.« less
  4. 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.more »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.« less
  5. Origami-inspired structures and material systems have been used in many engineering applications because of their unique kinematic and mechanical properties induced by folding. However, accurately modeling and analyzing origami folding and the associated mechanical properties are challenging, especially when large deformation and dynamic responses need to be considered. In this paper, we formulate a high-fidelity model — based on the iso-parametric Absolute Nodal Coordinate Formulation (ANCF) — for simulating the dynamic folding behaviors of origami involving large deformation. The centerpiece of this new model is the characterization of crease deformation. To this end, we model the crease using rotational spring at the nodes. The corresponding folding angle is calculated based on the local surface normal vectors. Compared to the currently popular analytical methods for analyzing origami, such as the rigid-facet and equivalent bar-hinge approach, this new model is more accurate in that it can describe the large crease and facet deformation without imposing many assumptions. Meanwhile, the ANCF based origami model can be more efficient computationally compared to the traditional finite element simulations. Therefore, this new model can lay down the foundation for high-fidelity origami analysis and design that involves mechanics and dynamics.