skip to main content


Title: Constrained Large-Displacement Thermal Analysis
Abstract Two different cases are encountered in the thermal analysis of solids. In the first case, continua are not subject to boundary and motion constraints and all material points experience same displacement-gradient changes as the result of application of thermal loads. In this case, referred to as unconstrained thermal expansion, the thermal load produces uniform stress-free motion within the continuum. In the second case, point displacements due to boundary and motion constraints are restricted, and therefore, continuum points do not move freely when thermal loads are applied. This second case, referred to as constrained thermal expansion, leads to thermal stresses and its study requires proper identification of the independent coordinates which represent expansion degrees-of-freedom. To have objective evaluation and comparison between the two cases of constrained and unconstrained thermal expansion, the reference-configuration geometry is accurately described using the absolute nodal coordinate formulation (ANCF) finite elements. ANCF position-gradient vectors have unique geometric meanings as tangent to coordinate lines, allowing systematic description of the two different cases of unconstrained and constrained thermal expansions using multiplicative decomposition of the matrix of position-gradient vectors. Furthermore, generality of the approach for large-displacement thermal analysis requires using the Lagrange–D'Alembert principle for proper treatment of algebraic constraint equations. Numerical results are presented to compare two different expansion cases, demonstrate use of the new approach, and verify its results by comparing with conventional finite element (FE) approaches.  more » « less
Award ID(s):
1852510
NSF-PAR ID:
10412705
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Computational and Nonlinear Dynamics
Volume:
18
Issue:
2
ISSN:
1555-1415
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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 simultaneously controlled using the new actuation force definitions. 
    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. 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
  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. null (Ed.)
    Granular-microstructured rods show strong dependence of grain-scale interactions in their mechanical behavior, and therefore, their proper description requires theories beyond the classical theory of continuum mechanics. Recently, the authors have derived a micromorphic continuum theory of degree n based upon the granular micromechanics approach (GMA). Here, the GMA is further specialized for a one-dimensional material with granular microstructure that can be described as a micromorphic medium of degree 1. To this end, the constitutive relationships, governing equations of motion and variationally consistent boundary conditions are derived. Furthermore, the static and dynamic length scales are linked to the second-gradient stiffness and micro-scale mass density distribution, respectively. The behavior of a one-dimensional granular structure for different boundary conditions is studied in both static and dynamic problems. The effects of material constants and the size effects on the response of the material are also investigated through parametric studies. In the static problem, the size-dependency of the system is observed in the width of the emergent boundary layers for certain imposed boundary conditions. In the dynamic problem, microstructural effects are always present and are manifested as deviations in the natural frequencies of the system from their classical counterparts. 
    more » « less