Implicit rigid body motion in curved finite elements
Rigid body motion implicit representation in curved finite elements
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Rigid body motion implicit representation in curved finite elements
The coupling of NASTRAN to another finite element program developed for the static analysis of automotive structures is discussed. The two programs were coupled together to use the substructuring capability of the in-house program and the normal mode analysis capability of NASTRAN. Modifications were made to the NASTRAN program in order to make the coupling feasible.
An implicit finite element numerical solution algorithm is derived for the compressible Navier-Stokes equations expressed in generalized coordinates. The theoretical basis utilizes a Galerkin-Weighted Residuals formulation, and extremization of approximation error within the context of a multipole expansion. A von Neumann analysis for a simplified form indicates the algorithm fourth- to sixth-order phase accurate, with third-order dissipation for the elementary linear element construction. Performance is improved for the algorithm constructed using quadratic interpolation. Numerical experiments for shocked duct flows are employed to optimize the several algorithm parameters. Additional numerical solutions validate algorithm accuracy and utility for aerodynamics applications.
A finite element formulation for the analysis of polymer processing is presented and its use in some typical situation including entry flow, transient Couette flow, and the Graetz (forced convection) problem is illustrated. The element formulations are constructed on the premise that momentum convection can be neglected (polymer melt flows typically have very low Reynolds' numbers), but that convective heat transfer may be significant (high Peclet numbers). Nonisothermal effects are considered important in polymer processing, due in part to the significant heating which may occur due to viscous dissipation, and also to the very strong influence of temperature on fluid viscosity. The flow is treated as Newtonian with the flow field being coupled to the heat transfer equation only through the viscous heat generation.
The finite element machine is a prototype computer designed to support parallel solutions to structural analysis problems. The hardware architecture and support software for the machine, initial solution algorithms and test applications, and preliminary results are described.
The Finite Element Machine at the NASA Langley Research Center is a prototype computer designed to support parallel solutions to structural analysis problems. The hardware architecture and support software for the machine, initial solution algorithms and test applications, and preliminary results are described. Directions for future work are presented.
The rate of convergence of the finite element method is a function of the strategy by which the number of degrees-of-freedom are increased. Alternative strategies are examined in the light of recent theoretical results and computational experience.
General nonlinear finite element codes were adapted for the purpose of analyzing the dynamics of gas turbine engines. In particular, this adaptation required the development of a squeeze-film damper element software package and its implantation into a representative current generation code. The ADINA code was selected because of prior use of it and familiarity with its internal structure and logic. This objective was met and the results indicate that such use of general purpose codes is viable alternative to specialized codes for general dynamics analysis of engines.
A finite element method (FEM) is applied to predicting coupled frequencies and pressures within an acoustic cavity in order to characterize sound transmission through a panel. Structural equations of motion are defined and the FEM model is configured with four-noded plate elements, the lateral displacement and two slopes being the unknowns at every node. Each element then has 12 degrees of freedom (DOF) and the displacement variation is expressed by a 12-term nonconforming polynomial. A consistent mass matrix is used to represent the panel mass matrix and a wave equation governs the acoustic volume. Analysis of pressure and displacement over the panel yields a square coupling matrix, and an eigenanalysis leads to a solution of the forced vibration problem.
Early versions of the assumed stress hybrid finite elements were based on the a priori satisifaction of stress equilibrium conditions. In the new version such conditions are relaxed but are introduced through additional internal displacement functions as Lagrange multipliers. A rational procedure is to choose the displacement terms such that the resulting strains are now of complete polynomials up to the same degree as that of the assumed stresses. Several example problems indicate that optimal element properties are resulted by this method.
A finite element analysis of the end notch flexure specimen for mode II interlaminar fracture toughness measurement was conducted. The effect of friction between the crack faces and large deflection on the evaluation of G sub IIc from this specimen were investigated. Results of this study are presented in this paper.
Finite element stiffness matrices are derived using an extended Hellinger-Reissner principle in which internal displacements are added to serve as Lagrange multipliers to introduce the equilibrium constraint in each element. In a consistent formulation the assumed stresses are initially unconstrained and complete polynomials and the total displacements are also complete such that the corresponding strains are complete in the same order as the stresses. Several examples indicate that resulting properties for elements constructed by this consistent formulation are ideal and are less sensitive to distortions of element geometries. The method has been used to find the optimal stress terms for plane elements, 3-D solids, axisymmetric solids, and plate bending elements.
Simple mixed finite element models and a computational procedure are presented for the solution of frictionless contact problems. The analytical formulation is based on a form of Reissner's large-rotation theory of the structure with the effect of transverse shear deformation included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the internal forces (or stress resultants), the generalized displacements, and the Lagrange multipliers associated with the contact conditions. Numerical examples of axisymmetric deformations of a hemispherical shell and planar deformations of a circular ring are presented. Both structures are pressed against a rigid plate. Detailed information about the response of the structures is presented. The numerical studies demonstrate the high accuracy of the mixed models and the effectiveness of the computational procedure based on combining the geometrically nonlinear terms and the contact conditions in one iteration loop.
Standard assumed displacement finite elements with anisotropic material properties perform poorly in complex stress fields such as combined bending and shear and combined bending and torsion. A set of three dimensional hybrid-stress brick elements were developed with fully anisotropic material properties. Both eight-node and twenty-node bricks were developed based on the symmetry group theory of Punch and Atluri. An eight-node brick was also developed using complete polynomials and stress basis functions and reducing the order of the resulting stress parameter matrix by applying equilibrium constraints and stress compatibility constraints. Here the stress compatibility constraints must be formulated assuming anisotropic material properties. The performance of these elements was examined in numerical examples covering a broad range of stress distributions. The stress predictions show significant improvement over the assumed displacement elements but the calculation time is increased.
The computation of engineering data from finite element solutions based on the virtual work formulation is discussed. It is shown that, using specially chosen virtual displacement functions, it is possible to compute certain functionals from the virtual work expression with much greater efficiency than would be possible by direct methods.
A 'brick' solid finite element is presently developed on the basis of the p-version analysis, and used to demonstrate the FEM concept of 'geometric entities'. This method eliminates interelement discontinuities between low- and high-order elements, allowing very fine control over the shape-function order in various parts of the model. Attention is given to the illustrative cases of a one-element model of an elliptic pipe, and a square cross-section cantilevered beam.
A Finite Element Method (FEM) is presented to determine propagation characteristics of deformed inflatable rectangular waveguide. Various deformations that might be present in an inflatable waveguide are analyzed using the FEM. The FEM procedure and the code developed here are so general that they can be used for any other deformations that are not considered in this report. The code is validated by applying the present code to rectangular waveguide without any deformations and comparing the numerical results with earlier published results.
Finite-element modeling (FEM) of 3C-SiC thin-film membranes on Si substrates was used to determine the residual stress and center deflection with applied pressure. The anisotropic, three-dimensional model includes the entire 3C-SiC membrane and Si substrate with appropriate material properties and boundary conditions. Residual stress due to the thermal-expansion-coefficient mismatch between the3C-SiC film and Si substrate was included in the model. Both before-and after-etching, residual stresses were calculated. In-plane membrane stress and normal deflection with applied pressure were also calculated. FEM results predict a tensile residual stress fo 259 MPa in the 3C-SiC membrane before etching. This decreases to 247 MPa after etching the substrate below the membrane. The residual stress experimentally measured on sample made at Case Western Reserve University was 280 MPa on post-etched membranes. This is excellent agreement when an additional 30-40 MPa of residual stress to account for lattice mismatch is added to the FEM results.