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.
ELPLOT program is a passive computer graphics system that could be utilized for display of models and responses of general finite-element analyses. Program includes: Wide range of view-orientation selections, number of alternative data-input formats, extensive family of finite-element types, and capabilities for both static and dynamic-response displays.
In a three-part series of papers, a generalized finite element methodology is formulated to handle traveling load problems involving large deformation fields in structure composed of viscoelastic media. The main thrust of this paper is to develop an overall finite element methodology and associated solution algorithms to handle the transient aspects of moving problems involving contact impact type loading fields. Based on the methodology and algorithms formulated, several numerical experiments are considered. These include the rolling/sliding impact of tires with road obstructions.
The viability of advanced viscoplastic models for nonlinear finite element analyses of structural components is investigated. Several uniaxial and a multiaxial problem are analyzed using the finite element implementation of Freed's viscoplastic model. Good agreement between the experimental and calculated uniaxial results validates the finite element implementation and gives confidence to apply it to more complex multiaxial problems. A comparison of results for a sample structural component (the cowl lip of a hypersonic engine inlet) with the earlier elastic, elastic-plastic, and elastic-plastic-creep analyses available in the literature shows that the elastic-viscoplastic analyses yield more reasonable stress and strain distributions. Finally, the versatility of the finite-element-based solution technology presented herein is demonstrated by applying it to another viscoplastic model.
A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.
The present paper describes a transient translating elastoplastic thermomechanical finite element model to study two-dimensional frictional rolling contact. Frictional two-dimensional contact is simulated by repeatedly translating a nonuniform thermomechanical distribution across the surface of an elastoplastic half space. The half space is represented by a two-dimensional finite element mesh with appropriate boundaries. Calculations are for an elastic-perfectly plastic material and the selected thermophysical properties are assumed to be temperature independent. The paper presents temperature variations, stress and plastic strain distributions and deformations. Residual tensile stresses are observed. The magnitude and depth of these stresses depends on (1) the temperature gradients and (2) the magnitudes of the normal and tangential tractions.
An independent refinement and integration procedure is developed to couple together independently modeled (global and local) regions in a single analysis. The models can have different levels of refinement and along the interface between them the finite element nodes need not coincide with one another. In the local model all the nodes except the nodes at the interface are statically condensed and the reduced stiffness matrix is obtained. For this static condensation a modified frontal solution technique is employed. A spline interpolation function that satisfies the linear isotropic plate bending differential equation is used to relate the local model interface nodal displacements to the global model interface displacements. The proposed independent refinement and integration procedure is evaluated by applying it to two- and three-dimensional cases involving inplane and out-of-plane deformation. The procedure yielded very accurate results for all the examples studied.
Conventional analysis of textile composites is impractical because of the complex microstructure. Global/local methodology combined with special macro elements is proposed herein as a practical alternative. Initial tests showed dramatic reductions in the computational effort with only small loss in accuracy.
Presented is an alternate cost-efficient and accurate elastic-plastic finite element procedure to analyze fatigue crack closure and its effects under general spectrum loading. Both Modes 1 and 2 type cycling loadings are considered. Also presented are the results of an investigation, using the newly developed procedure, of various factors that cause crack growth acceleration or retardation and delay effects under high-to-low, low-to-high, single overload, and constant amplitude type cyclic loading in a Mode 1 situation. Further, the results of an investigation of a centercracked panel under external pure shear (Mode 2) cyclic loading, of constant amplitude, are reported.
The aeroelastic stability of flap bending, lead-lag bending, and torsion of a helicopter rotor blade in hover is examined using a finite element formulation based on the principle of virtual work. Quasi-steady two-dimensional airfoil theory is used to obtain the aerodynamic loads. The rotor blade is discretized into beam elements, each with ten modal degrees of freedom. The resulting nonlinear equations of motion are solved for steady-state blade deflections through an iterative procedure. The flutter solution is calculated assuming blade motion to be a small perturbation about the steady solution. The normal mode method based on the coupled rotating natural modes about the steady deflections is used to reduce the number of equations in the flutter eigenanalysis. Results are presented for hingeless and articulated rotor blade configurations.
A computer program that generates three-dimensional (3D) finite element models for cracked 3D solids was written. This computer program, gensurf, uses minimal input data to generate 3D finite element models for isotropic solids with elliptic or part-elliptic cracks. These models can be used with a 3D finite element program called surf3d. This report documents this mesh generator. In this manual the capabilities, limitations, and organization of gensurf are described. The procedures used to develop 3D finite element models and the input for and the output of gensurf are explained. Several examples are included to illustrate the use of this program. Several input data files are included with this manual so that the users can edit these files to conform to their crack configuration and use them with gensurf.
In this paper a method for simulating electromagnetic fields scattered from complex objects is reviewed; namely, an unstructured finite element code that does not use traditional mesh partitioning algorithms.
A general computational method for analyzing unsteady flow around pitching and plunging airfoils was developed. The finite element method was applied in developing an efficient numerical procedure for the solution of equations describing the flow around airfoils. The numerical results were employed in conjunction with computer graphics techniques to produce visualization of the flow. The investigation involved mathematical model studies of flow in two phases: (1) analysis of a potential flow formulation and (2) analysis of an incompressible, unsteady, viscous flow from Navier-Stokes equations.
In a three-part series of papers, a generalized finite element solution strategy is developed to handle traveling load problems in rolling, moving and rotating structure. The main thrust of this section consists of the development of three-dimensional and shell type moving elements. In conjunction with this work, a compatible three-dimensional contact strategy is also developed. Based on these modeling capabilities, extensive analytical and experimental benchmarking is presented. Such testing includes traveling loads in rotating structure as well as low- and high-speed rolling contact involving standing wave-type response behavior. These point to the excellent modeling capabilities of moving element strategies.
The development of a nonlinear numerical algorithm for the analysis of stresses and displacements in partly wrinkled flat membranes, and its implementation on the SAP VII finite-element code are described. A comparison of numerical results with exact solutions of two benchmark problems reveals excellent agreement, with good convergence of the required iterative procedure. An exact solution of a problem involving axisymmetric deformations of a partly wrinkled shallow curved membrane is also reported.
The development of a nonlinear numerical algorithm for the analysis of stresses and displacements in partly wrinkled flat membranes, and its implementation on the SAP VII finite-element code are described. A comparison of numerical results with exact solutions of two benchmark problems reveals excellent agreement, with good convergence of the required iterative procedure. An exact solution of a problem involving axisymmetric deformations of a partly wrinkled shallow curved membrane is also reported.
Analytical formulations are developed to account for the coupled mechanical, electrical, and thermal response of piezoelectric composite materials. The coupled response is captured at the material level through the thermopiezoelectric constitutive equations and leads to the inherent capability to model both the sensory and active responses of piezoelectric materials. A layerwise laminate theory is incorporated to provide more accurate analysis of the displacements, strains, stresses, electric fields, and thermal fields through-the-thickness. Thermal effects which arise from coefficient of thermal expansion mismatch, pyroelectric effects, and temperature dependent material properties are explicitly accounted for in the formulation. Corresponding finite element formulations are developed for piezoelectric beam, plate, and shell elements to provide a more generalized capability for the analysis of arbitrary piezoelectric composite structures. The accuracy of the current formulation is verified with comparisons from published experimental data and other analytical models. Additional numerical studies are also conducted to demonstrate additional capabilities of the formulation to represent the sensory and active behaviors. A future plan of experimental studies is provided to characterize the high temperature dynamic response of piezoelectric composite materials.
Splitting, ultimate failure load and the damage path in center notched composite specimens subjected to in-plane tension loading are predicted using progressive failure analysis methodology. A 2-D Hashin-Rotem failure criterion is used in determining intra-laminar fiber and matrix failures. This progressive failure methodology has been implemented in the Abaqus/Explicit and Abaqus/Standard finite element codes through user written subroutines "VUMAT" and "USDFLD" respectively. A 2-D finite element model is used for predicting the intra-laminar damages. Analysis results obtained from the Abaqus/Explicit and Abaqus/Standard code show good agreement with experimental results. The importance of modeling delamination in progressive failure analysis methodology is recognized for future studies. The use of an explicit integration dynamics code for simple specimen geometry and static loading establishes a foundation for future analyses where complex loading and nonlinear dynamic interactions of damage and structure will necessitate it.
Assumed-displacement versions of the finite-element method are developed to predict large-deformation elastic-plastic transient deformations of structures. Both the conventional and a new improved finite-element variational formulation are derived. These formulations are then developed in detail for straight-beam and curved-beam elements undergoing (1) Bernoulli-Euler-Kirchhoff or (2) Timoshenko deformation behavior, in one plane. For each of these categories, several types of assumed-displacement finite elements are developed, and transient response predictions are compared with available exact solutions for small-deflection, linear-elastic transient responses. The present finite-element predictions for large-deflection elastic-plastic transient responses are evaluated via several beam and ring examples for which experimental measurements of transient strains and large transient deformations and independent finite-difference predictions are available.