skip to main content

Title: 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. 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
Award ID(s):
Publication Date:
Journal Name:
Journal of Computational and Nonlinear Dynamics
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 usingmore »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.« less
  2. The Finite Element Method (FEM) is widely used to solve discrete Partial Differential Equations (PDEs) in engineering and graphics applications. The popularity of FEM led to the development of a large family of variants, most of which require a tetrahedral or hexahedral mesh to construct the basis. While the theoretical properties of FEM basis (such as convergence rate, stability, etc.) are well understood under specific assumptions on the mesh quality, their practical performance, influenced both by the choice of the basis construction and quality of mesh generation, have not been systematically documented for large collections of automatically meshed 3D geometries. We introduce a set of benchmark problems involving most commonly solved elliptic PDEs, starting from simple cases with an analytical solution, moving to commonly used test problem setups, and using manufactured solutions for thousands of real-world, automatically meshed geometries. For all these cases, we use state-of-the-art meshing tools to create both tetrahedral and hexahedral meshes, and compare the performance of different element types for common elliptic PDEs. The goal of this benchmark is to enable comparison of complete FEM pipelines, from mesh generation to algebraic solver, and exploration of relative impact of different factors on the overall system performance. Asmore »a specific application of our geometry and benchmark dataset, we explore the question of relative advantages of unstructured (triangular/ tetrahedral) and structured (quadrilateral/hexahedral) discretizations. We observe that for Lagrange-type elements, while linear tetrahedral elements perform poorly, quadratic tetrahedral elements perform equally well or outperform hexahedral elements for our set of problems and currently available mesh generation algorithms. This observation suggests that for common problems in structural analysis, thermal analysis, and low Reynolds number flows, high-quality results can be obtained with unstructured tetrahedral meshes, which can be created robustly and automatically. We release the description of the benchmark problems, meshes, and reference implementation of our testing infrastructure to enable statistically significant comparisons between different FE methods, which we hope will be helpful in the development of new meshing and FEA techniques.« less
  3. 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 ofmore »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.« less
  4. 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.
  5. Intrinsic residual stresses in woven composites result from the coefficient of thermal expansion mismatch between the fibers and the matrix. Extrinsic residual stresses result from large scale thermal gradients during curing and cooling. Intrinsic residual stresses in 3D woven composites are sometimes severe enough to cause micro-cracking in the matrix. They are also expected to impact the fatigue resistance and the impact resistance. To the best of our knowledge, there have been no spatially resolved measurements of the intrinsic residual stress field as a function of position in the repeating weave pattern. We used digital image correlation (DIC) and electronic speckle pattern interferometry (ESPI) to measure the surface displacement field resulting from drilling a 1 mm diameter hole at four selected locations in two different 3D woven composite architectures that represent low and high through-the-thickness constraint. The two methods are used because the displacements sometimes on the lower end of the resolution for the DIC method and the displacement gradients are sometimes too steep to resolve the fringes for the ESPI method. Finite element models constructed with realistic fiber geometry using Dynamic Fabric Mechanic Analyzer software were utilized to estimate the residual stress field from cooling from the curing temperature.more »Holes were manually inserted by deactivating the elements in the hole region and the resultant displacement fields were compared to the measurements. In general, the measured displacement fields were lower in magnitude than the model predictions. In some cases, the sign of the predicted displacement field is opposite to the observed field which could be attributed to differences between the actual hole location and the hole in the model.« less