An official website of the United States government
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
Convergence Characteristics of Geometrically Accurate Spatial Finite Elements
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
- Award ID(s):
- 1852510
- PAR ID:
- 10217891
- Date Published:
- Journal Name:
- Journal of Computational and Nonlinear Dynamics
- Volume:
- 16
- Issue:
- 1
- ISSN:
- 1555-1415
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract We consider finite element approximations of the Maxwell eigenvalue problem in two dimensions. We prove, in certain settings, convergence of the discrete eigenvalues using Lagrange finite elements. In particular, we prove convergence in three scenarios: piecewise linear elements on Powell–Sabin triangulations, piecewise quadratic elements on Clough–Tocher triangulations and piecewise quartics (and higher) elements on general shape-regular triangulations. We provide numerical experiments that support the theoretical results. The computations also show that, on general triangulations, the eigenvalue approximations are very sensitive to nearly singular vertices, i.e., vertices that fall on exactly two ‘almost’ straight lines.more » « less
-
Abstract This article presents a novel derivation for the governing equations of geometrically curved and twisted three-dimensional Timoshenko beams. The kinematic model of the beam was derived rigorously by adopting a parametric description of the axis of the beam, using the local Frenet–Serret reference system, and introducing the constraint of the beam cross ection planarity into the classical, first-order strain versus displacement relations for Cauchy’s continua. The resulting beam kinematic model includes a multiplicative term consisting of the inverse of the Jacobian of the beam axis curve. This term is not included in classical beam formulations available in the literature; its contribution vanishes exactly for straight beams and is negligible only for curved and twisted beams with slender geometry. Furthermore, to simplify the description of complex beam geometries, the governing equations were derived with reference to a generic position of the beam axis within the beam cross section. Finally, this study pursued the numerical implementation of the curved beam formulation within the conceptual framework of isogeometric analysis, which allows the exact description of the beam geometry. This avoids stress locking issues and the corresponding convergence problems encountered when classical straight beam finite elements are used to discretize the geometry of curved and twisted beams. Finally, this article presents the solution of several numerical examples to demonstrate the accuracy and effectiveness of the proposed theoretical formulation and numerical implementation.more » « less
-
The folding motion of an origami structure can be stopped at a non-flat position when two of its facets bind together. Such facet-binding will induce self-locking so that the overall origami structure can stay at a pre-specified configuration without the help of additional locking devices or actuators. This research investigates the designs of self-locking origami structures and the locking-induced kinematical and mechanical properties. We show that incorporating multiple cells of the same type but with different geometry could significantly enrich the self-locking origami pattern design. Meanwhile, it offers remarkable programmability to the kinematical properties of the selflocking origami structures, including the number and position of locking points, and the deformation range. Self-locking will also affect the mechanical characteristics of the origami structures. Experiments and finite element simulations reveal that the structural stiffness will experience a sudden jump with the occurrence of self-locking, inducing a piecewise stiffness profile. The results of this research would provide design guidelines for developing self-locking origami structures and metamaterials with excellent kinematical and stiffness characteristics, with many potential engineering applications.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1115/1.4048731
