Stability of Non-Axisymmetric Rotor and Bearing Systems Modeled With Three-Dimensional-Solid Finite Elements

2019 ◽  
Vol 142 (1) ◽  
Author(s):  
Joseph Oh ◽  
Alan Palazzolo ◽  
Lingnan Hu
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Computational Efficiency ◽  
Degrees Of Freedom ◽  
Monodromy Matrix ◽  
Three Dimensional ◽  
Element Formulation ◽  
System Matrix ◽  
Rotor Bearing ◽  
Solid Finite Elements

Abstract Although rotors are simplified to be axisymmetric in rotordynamic models, many rotors in the industry are actually non-axisymmetric. Several authors have proposed methods using 3D finite element, rotordynamic models, but more efficient approaches for handling a large number of degrees-of-freedom (DOF) are needed. This task becomes particularly acute when considering parametric excitation that results from asymmetry in the rotating frame. This paper presents an efficient rotordynamic stability approach for non-axisymmetric rotor-bearing systems with complex shapes using three-dimensional solid finite elements. The 10-node quadratic tetrahedron element is used for the finite element formulation of the rotor. A rotor-bearing system, matrix differential equation is derived in the rotor-fixed coordinate system. The system matrices are reduced by using Guyan reduction. The current study utilizes the Floquet theory to determine the stability of solutions for parametrically excited rotor-bearing systems. Computational efficiency is improved by discretization and parallelization, taking advantage of the discretized monodromy matrix of Hsu's method. The method is verified by an analytical model with the Routh–Hurwitz stability criteria, and by direct time-transient, numerical integration for large order models. The proposed and Hill's methods are compared with respect to accuracy and computational efficiency, and the results indicate the limitations of Hill's method when applied to 3D solid rotor-bearing systems. A parametric investigation is performed for an asymmetric Root's blower type shaft, varying bearing asymmetry and bearing damping.

Modelling of rotors with three-dimensional solid finite elements

2001 ◽  
Vol 36 (4) ◽  
pp. 359-371 ◽  
Author(s):  
A Nandi ◽  
S Neogy
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Cross Sections ◽  
Three Dimensional ◽  
Element Model ◽  
Finite Element Formulation ◽  
Element Formulation ◽  
Two Dimensional ◽  
Inertia Effects ◽  
Solid Finite Elements

A shaft is modelled using three-dimensional solid finite elements. The shear-deformation and rotary inertia effects are automatically included through the three-dimensional elasticity formulation. The formulation allows warping of plane cross-sections and takes care of gyroscopic effect. Unlike a beam element model, the present model allows the actual rotor geometry to be modelled. Shafts with complicated geometry can be modelled provided that the shaft cross-section has two axes of symmetry with equal or unequal second moment of areas. The acceleration of a point on the shaft is determined in inertial and rotating frames. It is found that the finite element formulation becomes much simpler in a rotating frame of reference that rotates about the centre-line of the bearings with an angular velocity equal to the shafts spin speed. The finite element formulation in the above frame is ideally suited to non-circular shafts with solid or hollow, prismatic or tapered sections and continuous or abrupt change in cross-sections. The shaft and the disc can be modelled using the same types of element and this makes it possible to take into account the flexibility of the disc. The formulation also allows edge cracks to be modelled. A two-dimensional model of shaft disc systems executing synchronous whirl on isotropic bearings is presented. The application of the two-dimensional formulation is limited but it reduces the number of degrees of freedom. The three-dimensional solid and two-dimensional plane stress finite element models are extensively validated using standard available results.

scholarly journals Polygonal finite elements for three-dimensional Voronoi-cell-based discretisations

2012 ◽  
Author(s):  
Kaliappan Jayabal ◽  
Andreas Menzel
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Three Dimensional ◽  
Specific Property ◽  
Voronoi Cell ◽  
Polycrystalline Materials ◽  
Element Formulation ◽  
Hybrid Finite Element ◽  
Grain Structures ◽  
Polygonal Elements

Hybrid finite element formulations in combination with Voronoi-cell-based discretisation methods can efficiently be used to model the behaviour of polycrystalline materials. Randomly generated three-dimensional Voronoi polygonal elements with varying numbers of surfaces and corners in general better approximate the geometry of polycrystalline microor rather grain-structures than the standard tetrahedral and hexahedral finite elements. In this work, the application of a polygonal finite element formulation to three-dimensional elastomechanical problems is elaborated with special emphasis on the numerical implementation of the method and the construction of the element stiffness matrix. A specific property of Voronoi-based discretisations in combination with a hybrid finite element approach is investigated. The applicability of the framework established is demonstrated by means of representative numerical examples.

Three-Dimensional Solid Finite Element Contact Model for Rotordynamic Analysis: Experiment and Simulation

2020 ◽  
Vol 143 (3) ◽  
Author(s):  
Joseph Oh ◽  
Baik Jin Kim ◽  
Alan Palazzolo
Keyword(s):  
Surface Roughness ◽  
Contact Stiffness ◽  
Three Dimensional ◽  
Contact Model ◽  
Element Formulation ◽  
Interface Surface ◽  
Analysis Experiment ◽  
Rotor Bearing ◽  
Solid Finite Elements ◽  
Bearing System

Abstract Conventional rotordynamic analyses generally treat the rotor as a continuous body without considering effect of clamped joints. However, in modern rotating machines, rotors are often assembled with multiple complex-shaped parts and joints, which may significantly affect rotordynamic behavior. Several authors have proposed methods for implementing contact effects in rotordynamic analysis, but a more general modeling method for handling arbitrary contact geometries with various levels of surface roughness is needed. The present paper suggests a new contact model for rotordynamic analysis of an assembled rotor-bearing system with multiple parts connected by multiple joints. A contact element formulation is presented using solid finite elements and statistics-based contact theories. A test arrangement was developed to validate the proposed contact model for varying interface surface roughness and preloads. An iterative computation algorithm is introduced to solve the implicit relation between contact stiffness and stress distribution. Prediction results, using the contact model, are compared with measured natural frequencies for multiple configurations of a test rotor assembly. A case study is performed for an overhung type rotor-bearing system to investigate the effect of contact interfaces, between an overhung impeller and a rotor shaft, on critical speeds.

Parallelized Finite Element Analysis of Knitted Textile Mechanical Behavior

2018 ◽  
Vol 141 (2) ◽  
Author(s):  
D. Liu ◽  
S. Koric ◽  
A. Kontsos
Keyword(s):  
Finite Element Analysis ◽  
Finite Element ◽  
Mechanical Behavior ◽  
Computational Efficiency ◽  
High Performance ◽  
Degrees Of Freedom ◽  
Computational Models ◽  
Three Dimensional ◽  
Implicit Methods ◽  
Element Analysis

Direct numerical simulations (DNS) of knitted textile mechanical behavior are for the first time conducted on high performance computing (HPC) using both the explicit and implicit finite element analysis (FEA) to directly assess effective ways to model the behavior of such complex material systems. Yarn-level models including interyarn interactions are used as a benchmark computational problem to enable direct comparison in terms of computational efficiency between explicit and implicit methods. The need for such comparison stems from both a significant increase in the degrees-of-freedom (DOFs) with increasing size of the computational models considered as well as from memory and numerical stability issues due to the highly complex three-dimensional (3D) mechanical behavior of such 3D architectured materials. Mesh and size dependency, as well as parallelization in an HPC environment are investigated. The results demonstrate a satisfying accuracy combined with higher computational efficiency and much less memory requirements for the explicit method, which could be leveraged in modeling and design of such novel materials.

scholarly journals A structure-exploiting numbering algorithm for finite elements on extruded meshes, and its performance evaluation in Firedrake

2016 ◽  
Vol 9 (10) ◽  
pp. 3803-3815 ◽  
Author(s):  
Gheorghe-Teodor Bercea ◽  
Andrew T. T. McRae ◽  
David A. Ham ◽  
Lawrence Mitchell ◽  
Florian Rathgeber ◽  
...  
Keyword(s):  
Finite Element ◽  
Performance Evaluation ◽  
Finite Elements ◽  
Aspect Ratio ◽  
Three Dimensional ◽  
Data Layout ◽  
Element System ◽  
Continuous Galerkin ◽  
Performance Penalty ◽  
The Cost

Abstract. We present a generic algorithm for numbering and then efficiently iterating over the data values attached to an extruded mesh. An extruded mesh is formed by replicating an existing mesh, assumed to be unstructured, to form layers of prismatic cells. Applications of extruded meshes include, but are not limited to, the representation of three-dimensional high aspect ratio domains employed by geophysical finite element simulations. These meshes are structured in the extruded direction. The algorithm presented here exploits this structure to avoid the performance penalty traditionally associated with unstructured meshes. We evaluate the implementation of this algorithm in the Firedrake finite element system on a range of low compute intensity operations which constitute worst cases for data layout performance exploration. The experiments show that having structure along the extruded direction enables the cost of the indirect data accesses to be amortized after 10–20 layers as long as the underlying mesh is well ordered. We characterize the resulting spatial and temporal reuse in a representative set of both continuous-Galerkin and discontinuous-Galerkin discretizations. On meshes with realistic numbers of layers the performance achieved is between 70 and 90 % of a theoretical hardware-specific limit.

Truss and Beam Finite Elements Revisited: A Derivation Based on Displacement-Field Decomposition

1996 ◽  
Vol 11 (4) ◽  
pp. 371-380 ◽  
Author(s):  
Alphose Zingoni
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Symmetry Group ◽  
Displacement Field ◽  
Three Dimensional ◽  
Beam Elements ◽  
Symmetry Properties ◽  
Mass Matrices ◽  
General Displacement

Where a finite element possesses symmetry properties, derivation of fundamental element matrices can be achieved more efficiently by decomposing the general displacement field into subspaces of the symmetry group describing the configuration of the element. In this paper, the procedure is illustrated by reference to the simple truss and beam elements, whose well-known consistent-mass matrices are obtained via the proposed method. However, the procedure is applicable to all one-, two- and three-dimensional finite elements, as long as the shape and node configuration of the element can be described by a specific symmetry group.

Sign in / Sign up

Export Citation Format

Share Document