p=2 Continuous finite elements on tetrahedra with local mass matrix inversion

2012 ◽  
Vol 213-216 ◽  
pp. 289-298 ◽  
Author(s):  
Michael J. Brazell ◽  
Brian T. Helenbrook
Author(s):  
Karin Nachbagauer ◽  
Johannes Gerstmayr

For the modeling of large deformations in multibody dynamics problems, the absolute nodal coordinate formulation (ANCF) is advantageous since in general, the ANCF leads to a constant mass matrix. The proposed ANCF beam finite elements in this approach use the transverse slope vectors for the parameterization of the orientation of the cross section and do not employ an axial nodal slope vector. The geometric description, the degrees of freedom, and a continuum-mechanics-based and a structural-mechanics-based formulation for the elastic forces of the beam finite elements, as well as their usage in several static problems, have been presented in a previous work. A comparison to results provided in the literature to analytical solution and to the solution found by commercial finite element software shows accuracy and high order convergence in statics. The main subject of the present paper is to show the usability of the beam finite elements in dynamic and buckling applications.


Author(s):  
Alexander Olshevskiy ◽  
Oleg Dmitrochenko ◽  
Chang-Wan Kim

The present paper contributes to the field of flexible multibody systems dynamics. Two new solid finite elements employing the absolute nodal coordinate formulation are presented. In this formulation, the equations of motion contain a constant mass matrix and a vector of generalized gravity forces, but the vector of elastic forces is highly nonlinear. The proposed solid eight node brick element with 96 degrees of freedom uses translations of nodes and finite slopes as sets of nodal coordinates. The displacement field is interpolated using incomplete cubic polynomials providing the absence of shear locking effect. The use of finite slopes describes the deformed shape of the finite element more exactly and, therefore, minimizes the number of finite elements required for accurate simulations. Accuracy and convergence of the finite element is demonstrated in nonlinear test problems of statics and dynamics.


Author(s):  
Ahmed A. Shabana

Slope discontinuities and T-sections can be modeled in a straight forward manner using fully parameterized absolute nodal coordinate formulation (ANCF) finite elements that have a complete set of gradient vectors. Linear transformations that define the element connectivity can always be obtained and used to preserve ANCF desirable features that include constant mass matrix and zero Coriolis and centrifugal forces in the case of spinning structures. The objective of this paper is to develop a general method that allows for modeling slope discontinuities and T-sections using gradient deficient ANCF finite elements that do not have a complete set of coordinate lines and gradient vectors. Linear connectivity conditions that preserve all the ANCF desirable features including the constant mass matrix are developed at the nodes of slope discontinuities. At these nodes of discontinuity, one can always define a complete set of independent coordinate lines that lie on the structure. These coordinate lines can be used to define a complete set of independent gradient vectors at these nodes. Since the proposed method is based on linear coordinate transformations, the method can be implemented in a preprocessor computer program. The application of the proposed general method is demonstrated using ANCF gradient deficient beam element example.


1992 ◽  
Vol 114 (1) ◽  
pp. 3-9 ◽  
Author(s):  
R. C. Engels

A method is described to model the dynamics of finite elements. The assumed modes method is used to show how static shape functions approximate the element mass distribution. The deterioration of the modal content of a model can be linked to the neglect of interface restrained assumed modes. Restoration of a few of these modes leads to higher accuracy with fewer generalized coordinates compared to the standard consistent mass matrix approach. Also, no need exists for subdivision of basic elements such as rods and beams. The mass and stiffness matrices for several basic elements are derived and used in demonstration problems.


Author(s):  
Alan S. K. Kwan

The stiffness relationship and the distributed mass matrix for a geometrically nonlinear three dimensional straight axial element is derived for use in prestressed cablenet structures. The justification for the use of a linearised stiffness relationship is provided through a theoretical derivation. Results using this simple element have shown a high degree of correlation with results to those available in the literature obtained with more complex curved finite elements, analogous membrane models and other techniques.


Author(s):  
Karine Reaud ◽  
Claude Vallee ◽  
Danielle Fortune

Abstract Solutions of the vibrations equations are generally obtained, in the linear case, by methods involving either matrix exponential computation or matrix eigendecomposition. However, these methods lead to a loss of symmetry because of the necessary inversion of the mass matrix. In doing so, one can introduce ill-conditioned matrices and, thus, compute eigenvectors with poor accuracy. In order to avoid these inconveniences, our method, which is an extension of Le Verrier-Souriau algorithm, provides solutions to the vibrations equations without inverting the mass matrix or computing eigenvectors. Moreover, we can solve damping systems even when the damping matrix has no specific properties such as the Basile property.


Sign in / Sign up

Export Citation Format

Share Document